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In this paper a description of the Hubbard model on the square lattice with nearest-neighbor 
transfer integral t, on-site repulsion U, and Na ^ 1 sites consistent with its exact global SO{S) x 
SO{3) X U{1) symmetry is constructed. Our studies profit from the interplay of that recently found 
global symmetry of the model on any bipartite lattice with the transformation laws under a suitable 
electron - rotated-electron unitary transformation of a well-defined set of operators and quantum 
objects. For U /At > the occupancy configurations of these objects generate the energy eigenstates 
that span the one- and two-electron subspace. Such a subspace as defined in this paper contains 
nearly the whole spectral weight of the excitations generated by application onto the zero-spin- 
density ground state of one- and two-electron operators. Our description involves three basic objects: 
charge c fermions, spin-1/2 spinous, and 77-spin-l/2 7)-spinons. Independent spinous and independent 
7)-spinons are invariant under the above unitary transformation. Alike in chromo dynamics the quarks 
have color but all quark-composite physical particles are color-neutral, the r;-spinon (and spinous) 
that are not invariant under that transformation have rj spin 1/2 (and spin 1/2) but are part of 77-spin- 
neutral (and spin-neutral) 2!/-r7-spinon (and 2!/-spinon) composite rjv fermions (and sp fermions) 
where v = 1,2, ... is the number of r;-spinon (and spinon) pairs. The occupancy configurations of the c 
fermions, independent spinous and 2i^-spinon composite sv fermions, and independent 77-spinons and 
2z/-?7-spinon composite rjv fermions correspond to the state representations of the 17(1), spin SU{2), 
and »7-spin SU{2) symmetries, respectively, associated with the model 50(3) x 50(3) x U{1) = 
[SU{2) X SU{2) X U{l)]/Zi global symmetry. The components of the au fermion discrete momentum 
values qj = [Qj ^i, Qj ^2] eigenvalues of the corresponding set of av translation generators in 
the presence of fictitious magnetic fields Bav. Our operator description has been constructed to 
inherently the ai/ translation generators q^^ in the presence of the fictitious magnetic field Bav 
commuting with the momentum operator, consistently with their component operators cjavxi and 
ijai^xi commuting with each other. In turn, unlike for the ID model such generators not commute 
in general with the Hamiltonian, except for the Hubbard model on the square lattice in the one- 
and two-electron subspace. Concerning one- and two-electron excitations, the picture that emerges 
is that of a two-component quantum liquid of charge c fermions and spin-neutral two-spinon si 
fermions. The description introduced here is consistent with a Mott-Hubbard insulating ground 
state with antiferromagnetic long-range order for half filling at a; = hole concentration and a 
ground state with short-range spin order for a well-defined range of finite hole concentrations 2; > 0. 
For < a; ^ 1 the latter short-range spin order has an incomensurate-spiral character. Our results 
are of interest for studies of ultra-cold fermionic atoms on optical lattices and elsewhere evidence 
is provided that upon addition of a small three-dimensional anisotropy plane-coupling perturbation 
to the square-lattice quantum liquid considered here its short-range spin order coexists for low 
temperatures and a well-defined range of hole concentrations with a long-range superconducting 
order so that the use of the general description introduced in this paper contributes to the further 
understanding of the role of electronic correlations in the unusual properties of the hole-doped 
cuprates. 

PACS numbers: 71.10.Fd, 71.10.-l-w, 71.27.-fa, 71.10.Hf, 71.30.-fh 



I. INTRODUCTION 



There is some consensus that the Hubbard model on a square lattice is the simplest toy model for describing the 
effects of electronics correlations in high- Tc superconductorsir— and their Mott-Hubbard insulators parent compounds^. 
In addition, the model can be experimentally realized with unprecedented precision in systems of ultra-cold fermionic 
atoms on an optical lattice and one may expect very detailed experimental results over a wide range of parameters 
to be available^. That includes recent studies on systems of ultra-cold fermionic atoms describing the Mott-Hubbard 
insulating phase of the the Hubbard model on a cubic lattice-. Unfortunately, the model has no exact solution and 
many open questions about its properties remain unsolved. A recent exact result, which applies to the model on any 
bipartite lattice^, is that for on-site repulsion [/ > its global symmetry is 50(3) x SO{3) x C/(l) = [SO{A) x U{l)]/Z2. 
That is an extension of the model well-known 5*0(4) symmetrjsiSi, which becomes explicit provided that one describes 
the problem in terms of the rotated electrons obtained from the electrons by any of the unitary transformations of 
the type considered in Refsi^iii. 
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In this paper a general description is introduced for the Hubbard model on the two-dimensional (2D) square lattice 
with spacing a, = [Na]^ sites, D — 2, Na even, and lattice edge length L = Na a. (Due to editorial decisions, the 
results presented in this paper were distributed over three articles, specifically Refs.— . We thus forward the reader 
to these references, in which the presentation of the results was further improved, yet the scientific message remains 
the same.) 

The representation introduced here refers to suitable quantum objects whose occupancy configurations generate the 
state representations of the group SO{3) x SO{3) x U{1). Addition of chemical-potential and magnetic-field operator 
terms to the Hamiltonian lowers its symmetry. Such terms commute though with it and the momentum and energy 
eigenstates correspond to representations of that group. It is shown in Appendix A that all the physics of the model 
in the whole Hilbert space can be obtained from that of the model in the subspace spanned by the lowest-weight 
states (LWSs) of both the yy-spin SU{2) and spin SU{2) algebras. Hence often our results refer to that large subspace. 
Another subspace that plays a key role in our description is the one- and two-electron subspace defined in Section V. 
For 2D, the states generated by the momentum occupancy configurations of the above objects belonging to such a 
subspace are for J7/4i > energy eigenstates. 

Furthermore, often we discuss the use of our description for the model on an one-dimensional (ID) lattice as well, 
so that the coefficient D appearing in several quantities considered in this paper may have the values D = 1,2. Due 
to the integrability of the ID model, the states generated by the momentum occupancy configurations of the above- 
mentioned objects are exact energy eigenstates for the whole Hilbert space and all values of U/At. Our results refer 
to a very large system whose lattice has ^ 1 sites and for the square (and ID) lattice torus periodic boundary 
conditions (and periodic boundary conditions) are used. The description involves the generalization to all finite values 
of the on-site repulsion [/ > of the exact transformation for separation of spin-1/2 fermions without constraints 
introduced for large U values in Reii^. 

The global 5*0(3) x S0{3) x U{1) symmetry refers to the Hubbard model on any bipartite lattice. Hence it applies 
to the Z?-dimcnsional cubic latticei^ii^. However, some of the results obtained below are valid for the square lattice 
only. In spite of our study focusing onto the model on that lattice, the reason why often our analysis refers to the 
corresponding ID problem as well is that in contrast to real-space dimensions D > 1 there is an exact solution for 
the model on a ID latticei^"— . Our quantum-object description also applies to ID, yet due to the integrability of 
the model often it then leads to a different physics. In turn, one of the procedures used to control the validity of the 
approximations used in the construction of our description profits both from identifying the common proprieties of the 
model on the square and ID lattices and checking whether they are valid for the model on the latter lattice by means 
of its exact solution. Some of these common properties are related to the model commuting with the generators of 
the group SO{3) x S0{3) x U{1) both for the square and ID bipartite lattices. 

For such two lattices the fulfillment of our program reveals the emergence of three basic quantum objects: spin-less 
and 77-spin-less charge c fermions, spin-1/2 spinons, and 7y-spin-l/2 ry-spinons. In principle there are some occupancy 
configurations of such objects that generate the exact momentum and energy eigenstates for [//4i > 0, yet the 
detailed structure of such configurations is a complex problem. It simplifies in the one- and two-electron subspace 
constructed in this paper, where for the model on the square lattice and U /At > the occupancy configurations 
of the above objects generate exact ground states and excited energy eigenstates. Such a subspace contains nearly 
the whole spectral weight generated by application onto the zero-spin-density ground state of one- and two-electron 
operators. Throughout this paper we often call that subspace one- and two- electron subspace. In referencsiS^ the 
spinon occupancy configurations of the states that span that subspace are studied. 

The composite 2z^-77-spinon rjiy fermions (and 2z^-spinon sv fermions) emerge from ry-spin (and spin) neutral rjv 
bond particles (and sv bond particles) through extended Jordan- Wigner transformations'^^—. For the model on the 
square lattice the two-spinon si bond particles are found in Refsiii^ to be closely related to the spin-singlet bonds of 
Ref.^"^. It is confirmed in Appendix B that for ID the discrete momentum values of the c fermions and aiy fermions are 
good quantum numbers for all U /At values and their momentum occupancy configurations generate the exact energy 
eigenstates. Evidence that for the one- and two-electron subspace the states generated by the occupancy configurations 
of such discrete momentum values are for the model on the square lattice and U/At > energy eigenstates is given in 
this paper and in Reffi. 

The paper is organized as follows. The model, the uniquely defined rotated-electron description used in our studies, 
its limitations and usefulness, and the global SO{3) x SO{3) x U{1) symmetry are the subjects of Section II. In Section 
III the c fermion, 77-spin-l/2 77-spinon, and spin-1/2 spinon and corresponding c, 77-spin, and spin effective lattices 
arc introduced. The vacua of the theory, their transformation laws under the electron - rotated-electron unitary 
transformation, and the subspaces they refer to are issues also addressed in that section. The composite ai' bond 
particles and av fermions, corresponding av effective lattices, ground-state occupancy configurations, and a suitable 
complete set of momentum eigenstates are the subjects of Section IV. In Section V the one- and two-electron subspace 
where only the c fermions and si fermions play an active role is introduced. The corresponding c and si effective 
lattices and elementary excitations are discussed. In Section VI it is shown that the description introduced in this 
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paper is consistent with a Mott- Hubbard insulating ground state with long-range antiferromagnetic order for half filling 
and a ground state with short-range spin order for a well-defined range of finite hole concentrations. Strong evidence 
is given that for < x ^ 1 and intermediate and large values of U/it the latter order has an incommensurate-spiral 
character. Section VII contains the concluding remarks. Moreover, that the whole physics can be extracted from the 
Hubbard model in the LWS subspace spanned by the ry-spin and spin LWSs is confirmed in Appendix A. In Appendix 
B it is confirmed that the general operator description introduced in this paper for the Hubbard model on the square 
and ID lattices is consistent with the exact solution of the ID problem. For instance, it is confirmed that the discrete 
momentum values of the c fermions and composite av fermions are good quantum numbers. Furthermore, in that 
Appendix the relation of the creation and annihilation fields of the charge ABCDF algebra^^ and more traditional 
spin ABCD Faddeev-Zamolodchikov algebra^** of the algebraic formulation of the ID exact solution of Ref.^^ to the c 
and av fermion operators is discussed and the consistency between the two corresponding operational representations 
confirmed. 



II. THE MODEL, A SUITABLE ROTATED-ELECTRON DESCRIPTION, AND THE GLOBAL 

SO{3) X SO{3) X (7(1) SYMMETRY 

A. The model, a suitable rotated-electron description, and its limitations and usefulness 

The Hubbard model is given by, 

H^tf+U[N^~ Q]/2 ; f = - T.^4,.a ^r, ,a + h.c\ Q ^ J2 J2 ""r.^a (1 - • (1) 

(rjT-,) cr j = l o-=t,i 

Here D = \ and D = 2 for the ID and square lattices, respectively, t is the nearest-neighbor transfer integral, T the 
kinetic-energy operator in units of Q the operator that counts the number of electron singly occupied sites so that 
the operator t) = [N —Q]/2 counts the number of electron doubly occupied sites, np.^c! — ct. „cp^,<T where —a (and 

—a =1) for (T =1 (and a =t), N = and N^, — "-r,.cr- The kinetic-energy operator T can be expressed 

in terms of the operators, 

(rj-r,.,) o- 

as T = To + T+i + T'„i. These three kinetic operators play an important role in the physics. The operator Tq does 
not change electron double occupancy whereas the operators T+i and T_i do it by -t-1 and —1, respectively. 

The studies of Refi^ consider unitary operators V = V{U /At) and corresponding rotated-electron operators ct ^ = 
1/t cjj^ ^ V such that rotated-electron single and double occupancy are good quantum numbers for U /At > 0. The 
generator Sc of the global C/(l) symmetry given below and in that reference has eigenvalue Sc with 2Sc being the 
number of rotated-electron singly occupied sites. Most choices of the unitary operators V correspond to choices 
of U/At — oo sets {l^'oo)} of 4^° energy eigenstates such that the states |^c//4t) = V^^|*oo) are not energy and 
momentum eigenstates for finite U/At values yet belong to a subspace with constant and well-defined values of Sc, 5,,, 
and Ss- Here S*,, (and Sg) denotes the 77-spin (and spin) value of the energy eigenstates. We call S^j = — [N/;^—N]/2 (and 
51 = — [A^t ~ the corresponding projection. We focus our attention on ground states with hole concentration 

X = [A^jP — N]/N//' > and spin density m = [A'^ — Ni\/N//' = and are particularly interested in the LWS-subspace 
spanned by the LWSs of both the ry-spin and spin algebras. Such energy eigenstates refer to values of Sa and 
such that Sa = —S^ for a = r],s. Indeed, since the off-diagonal generators that generate the non-LWSs from the 
LWSs commute with the operator V"—, the whole physics can be extracted from the model ((T|) in the LWS-subspace, 
as confirmed in Appendix A. 

Let {|4'c//4t)} be a complete set of 4^" energy, momentum, yy-spin, r]-spin projection, spin, and spin-projection 
eigenstates for U/At > 0. In the limit U/At 00 such states correspond to one of the many choices of sets {|^'oo)} 
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of 4^'' C//4t-infinite energy eigenstates. For this choice there exists exactly one unitary operator V = V{U /At) such 
that |^'f7/4t) — V^^l^oo)- Here we consider such a unitary operator and corresponding rotated-electron operators 
ct ^ = ct. ^ V. The states \'^u/4t) = l^^l^oo) (one for each value of U/At > 0) that are generated from the same 
initial state |^oo) belong to the same V tower. 

We associate with any operator O a rotated operator O — V'' O V that has the same expression in terms of rotated- 
electron creation and annihilation operators as O of electron creation and annihilation operators, respectively. Any 
operator O can be written in terms of rotated-electron creation and annihilation operators as, 

6 = vdv^ ^d + [d, s] + ^[[d,s], s] + ...^d + [d, s] + ^[[d,s], s] + ..., 



t [ - 



U 



+ 0{t'/U^)- S=^ r+i-T_i +0{i'/U''). (3) 



The operator S appearing in this equation is related to the unitary operator as V'^ — and V = e^'^. The general 
unknown expression of S involves only the kinetic operators To , T_|_ i , and r_ i of Eq. ([2]) and numerical U/At dependent 
coefficients and for U/At > can be expanded in a series of t/U. The corresponding first-order term has the universal 
form given in Eq. ([3]). To arrive to the expression in terms of the operator S also given in Eq. ([3]), the property 

that the operator V commutes with itself so that V = = V = e~^' and S = S and both the operators V and S 
have the same expression in terms of electron and rotated-electron creation and annihilation operators is used. This 
is behind the expansion S — {t/U) [T+i — T-i] + 0{t^ /U"^) given in that equation for the operator S. Its higher-order 
terms can be written as a product of operator factors whose expressions involve the rotated kinetic operators Tq, T+i, 
and T_i. The full expression of the operator S = S can for U/At > be written as S' = 5'(oo) + 5S = S{oo) + SS 
where S{oo) = S{oo) corresponds to the operator S{1) at Z = oo defined in Eq. (61) of Ref and SS = SS has the 
general form provided in its Eq. (64). The unitary operator V ~ e^^ = considered here corresponds to exactly 
one choice of the coefficients D'^'^^m.) of that equation, where for the operator V the index k = 1,2, ... refers to the 
number of rotated-electron doubly occupied sites. The problem of finding the explicit form of the operators V — V 
and S — S is equivalent to finding all coefficients D^'^^m) associated with the electron - rotated-electron unitary 
transformation as defined in this paper. 

Since for finite U/At values the Hamiltonian H of Eq. ([1]) does not commute with the unitary operator V — , 
when expressed in terms of rotated-electron creation and annihilation operators it has an infinite number of terms 
and according to Eq. (O reads, 

H ^VHV^ =H+[H, S] + ^ [[H, S], S] + ... . (4) 

The commutator [H, S ] does not vanish except for U /At ^ oo so that H ^ H for finite values of U/At. For U /At very 
large the Hamiltonian of Eq. ^ corresponds in terms of rotated-electron creation and annihilation operators to a 
simple rotated-electron t — J model. In turn, the higher-order t/U terms, which become increasingly important upon 
decreasing U /At, generate effective rotated-electron hopping between second, third, and more distant neighboring 
sites. Indeed, the products of the kinetic operators Tq, I+i, and r_i contained in the higher-order terms of = 
{t/U) [r+i — T_i] + 0{t^ /U"^) also appear in the Hamiltonian expression Q in terms of rotated-electron creation and 
annihilation operators. In spite of the operators Tq, T+i, and T_i generating only rotated-electron hoping between 
nearest-neighboring sites, their products generate effective hoping between for instance second and third neighboring 
sites whose real-space distance in units of the lattice spacing a is for the model on the square lattice \/2a and 2 a, 
usually associated with transfer integrals t' and t" , respectively^^. 

It follows that when expressed in terms of the rotated-electron operators emerging from the specific unitary trans- 
formation considered above the simple Hubbard model ([1} as given in Eq. ^ contains Hamiltonian terms associated 
with higher-order contributions which can be effectively described by transfer integrals t' , t" , and of higher order. 
It is found in the following that for hole concentration equal to or larger than zero both the model ground state 
and the excited states that span the one- and two-electron subspace as defined in this paper have vanishing rotated- 
electron double occupancy. It follows from a Renormalization Group analysis based on the results of Refs.^^'^'' that 
for intermediate and large values of U/At obeying approximately the inequality U/At > uq fa 1.302, besides the origi- 
nal nearest-neighboring hoping processes only those involving second and third neighboring sites are relevant for the 
square-lattice quantum liquid described by the Hamiltonian of Eqs. ([ij and (j4]) in the one- and two-electron subspace. 
The value U /At = ~ 1.302 is that where according to the results of Section VI the si fermion spinon pairing energy 
reaches its maximum magnitude for small hole concentration < x 1. Such a pairing energy is found to vanish 
both in the limits U/At ^ and U /At — >■ oo. In contrast, for the ID model it vanishes for all U/At values. 
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Hence for approximately U/At > uq, out of the infinite terms on the right- hand-side of Eq. (U) only the first few 
Hamiltonian terms play an active role in the physics of the Hubbard model on the square lattice in the one- and 
two-electron subspace. Therefore, for intermediate and large values of U/At such a square-lattice quantum liquid 
can be mapped onto an effective t — J model on a square lattice with t, t' = t'{U/At), and t" = t"{U/4£) transfer 
integrals^^ for which the role of the processes associated with t' = t'{U/4:t) and t" — t"{U/At) becomes increasingly 
important upon decreasing the U/4t value. For that U/At range the latter model is equivalent to the Hubbard model 
on the square lattice (j4|) in the subspace under consideration and expressed in terms of rotated-electron creation and 
annihilation operators. Indeed, the t — J model constraint against double occupancy is in that subspace equivalent 
to expressing the Hubbard model in terms of rotated-electron creation and annihilation operators. 

However, the rotated electrons are not the ultimate objects whose occupancy configurations generate the exact 
ground states and excited energy eigenstates. For U/At > suitable occupancy configurations of the charge c 
fermions and spin- neutral two-spinon si fermions simpler than those of the rotated electrons generate such states: 
for the square-lattice model in the one- and two-electron subspace the excited states generated by the occupancy 
configurations of such objects are energy eigenstates. (Actually, for ID they refer to exact energy eigenstates for the 
whole Hilbert space provided that U/At > 0.) Another advantage relative to the rotated electrons is that the processes 
associated with the above effective transfer integrals t' — t'{U/At) and t" — t"{U/At) are implicitly contained in the 
U /At dependence of the energy scales associated with the c and si fermions spectra derived in Ref.^. 

The general operator description introduced in this paper for the Hubbard model on the square lattice has two 
main limitations: 

I) For small and intermediate values of U /At the explicit form of the unitary operator V associated with the rotated- 
electron operators as defined in this paper remains an open problem. It is known that such a unitary operator has for 
U/At > the general form V = where the expression of S involves only the three kinetic operators Tq, T+i, and 
T_i of Eq. However, the finding of the explicit form of the J7/4t-dependent functional of S in terms of the latter 
three operators, valid for the whole range of finite U/At values, is a very involved and unsolved quantum problem, 
beyond the reach of the present status of our scheme. Indeed, the operator description introduced in this paper has 
been constructed to inherently the solution of that problem being equivalent to the solution of the Hubbard model 
on the square lattice. 

II) It turns out that the quantum problem under consideration is non-perturbative in terms of electron operators 
so that, in contrast to a three-dimensional (3D) isotropic Fermi liquid^, rewriting the square-lattice quantum- liquid 
theory emerging from the general description introduced in this paper for the model in the one- and two-electron 
subspace in terms of the standard formalism of many-electron physics is in general an extremely complex problem. 
Fortunately, such a quantum liquid dramatically simplifies when expressed in terms of the c fermion and si fermion 
operators and one can rich limited but valuable information on its physics in spite of the lack of explicit information 
about the matrix elements between energy eigenstates. The point is that as justified below in the above subspace 
the c and si fermion discrete momentum values have been constructed to inherently being good quantum numbers 
so that the interactions of these objects are residua l "'^'^^i^^ . 

Concerning the Hamiltonian, the microscopic processes corresponding to the effective transfer integrals t' = t'{U /At) 
and t" — t"{U/At) are important to characterize the type of order associated with the phases of the square-lattice 
quantum liquid as for instance the short-range incomensurate-spiral spin order considered in Subsection VI-B. Fortu- 
nately, for small finite hole concentrations, it is confirmed in Subsection VI-A that for simple one- and two-electron 
operators O the leading operator term O on the right-hand side of Eq. ^ generates nearly the whole spectral weight. 
Hence in spite of the limitations I and II, our description provides useful information about the physics contained 
in the model on the square lattice. Indeed, there are several reasons why, in spite of both the explicit form of the 
unitary operator V being known only for large values of U/At and the difficulties in rewriting the theory emerging 
from the description introduced here in terms of the standard formalism of many-electron physics, that description is 
rather useful for finite values of U/At and particularly for approximately U/At > uq. This is confirmed in Ref.^ from 
comparison of results obtained by our description and the standard formalism of many-electron physics concerning one 
of the few problems where controlled and reliable approximations exist within the latter formalism for the Hubbard 
model on the square lattice. 

For one- and two-electron operators O the terms of the general expression ([3]) containing commutators involving the 
related operator S = S are found in Section V-B to generate very little spectral weight. Hence one can reach a quite 
faithful representation of such operators by expressing them in terms of the c and si fermion operators. In Subsection 
III-A we provide the expression of the rotated-electron creation and annihilation operators in terms of the operators 
of the c fermions, 77-spinons, and spinous. The construction of the spin- neutral two-spinon composite si bond-particle 
operators and corresponding si fermion operators is a more involved problem addressed in Refs.^'^^. Furthermore, in 
spite of the explicit form of the unitary operator V being imknown except for U/At ^ 1, we are able to access the 
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transformation laws and/or invariance of several operators and quantum objects under the corresponding electron - 
rotated-electron unitary transformation, what provides valuable information about the physics of the Hubbard model 
on the square lattice for C//4t > 0. 

The physical reason why the operator terms of the general expression ([3|) containing commutators involving the 
operator S = S generate very little one- and two-electron spectral weight is that the rotated electrons as defined in 
this paper are directly related to objects whose interactions are residual so that the leading elementary processes in 
terms of them generate nearly the whole spectral weight. Expression of the operator terms of Eq. ^ containing 
commutators involving the operator S* = S* in terms of the operators of such objects reveals that those generate 
higher-order processes, whereas the elementary and leading processes are generated by the operator O. 

In Subsection VI-B we study the U/At dependence of the maximum magnitude 2Aq of the si fermion spinon pairing 
energy of the Hubbard model on the square lattice for small hole concentrations < a; <C 1. Such an energy scale is 
also well defined for a larger range of finite hole concentrations yet then has a different physical meaning and plays an 
important role in square-lattice quantum liquid. For < a; <C 1 it plays the role of short-range incommensurate-spiral 
spin order parameter. That order is absent in the ID model for which 2Ao = 0. For the square-lattice model it occurs 
both for < a; < 1 and temperatures < T < Tg and at half filling a; = for < T < , where Tq « Ao/^b. Our 
study of the J7/4t dependence of the energy parameter 2Aq profits from combination of the description introduced 
here with results of the low-temperature approach of Refi^ to the half-filled model. Such an approach uses a CP^ 
representation and is valid below the temperature Tg called in that reference , which marks the onset of the short- 
range spin order. For the specific value U/At « 1.525 the magnitude of 2Ao provided in Subsection VI-B is obtained 
in Refiii from combination of the description introduced in this paper with the results of Refi^ for the spin-wave 
spectrum obtained by summing up an infinite set of ladder diagrams for the half-filled model in a spin-density-wave- 
broken symmetry ground state. Moreover, in the same subsection the spiral-incommensurate character of such a 
short-range spin order is studied by combining results obtained from the use of our description with the action and 
corresponding CP^ field theory of Refs.^^'?I. 

A realistic program to be carried out here and in following papers includes the fulfillment of two main tasks: (i) 
Introducing the objects whose occupancy configurations generate the energy eigenstates that span the one- and two- 
electron subspace and writing the rotated-electron creation and/or annihilation operators in the expression of the one- 
and two-electron operators O of Eq. ([3]) in terms of the operators of such objects; (ii) Finding the specific occupancy 
configurations of the latter objects that generate such states. 

As discussed above, the expression of the operator O in terms of the c fermion, ry-spinon, and spinon operators 
provides a rather good representation of the corresponding one- or two-electron operator O itself. In this paper we 
fulfill task (i) and start solving the problems needed for the fulfillment of task (ii). The latter task is concluded in 
Rcfs.^'iS,, whose studies profit from the general description introduced here. Furthermore, the further construction of 
the square-lattice quantum liquid is fulfilled in Refsjii^^. The lack for U/4t > of the explicit general form of the 
unitary operator V prevents the derivation of matrix elements between energy eigenstates. It is confirmed elsewhere 
that the information about the Hubbard model on the square lattice in the one- and two-electron subspace reached 
by combining different methods and approximations with the general operator description introduced in this paper 
clarifies important open issues of the physics contained in the model. That includes its usefulness and suitability to 
describe the physics of real materiala^ i^^''^^ . 

B. The model global SO{3) x SO(3) x U{1) symmetry 

The studies of Refj2^ reveal that for U/At oo the Hamiltonian ([T|) has a local SU{2) x SU{2) x [/(I) gauge 
symmetry. In turn, for finite U/At values only its interacting term is invariant under SU{2) x SU{2) x U{1) and 
therefore SU{2) x SU{2) x U{1) becomes a group of permissible unitary transformations. As discussed in that 
reference, the local U{1) canonical transformation is not the ordinary U{1) gauge subgroup of electromagnetism, but 
instead is a "nonlinear" transformation, i.e., a gauge transformation transforming a creation or annihilation operator 
into a sum of nonlinear polynomials of such operators. 

Furthermore, the studies of Ref.^ reveal that for on-site repulsion U > the above local SU{2) x SU{2) x U{1) gauge 
symmetry of the model for J7/4t ^ oo can be lifted to a global [SU{2) x SU{2) x U{l)]/Z^ = SO{3) x S0{3) x U{1) 
symmetry in the presence of the kinetic-energy hopping term of the Hamiltonian ([1]) with t > and U/At > 0. 
The generator of the new found hidden global U{1) symmetry is one half the rotated-electron number of singly- 
occupied sites operator. As mentioned above, addition of chemical-potential and magnetic-field operator terms to 
the Hamiltonian ([T]) lowers its symmetry. However, such terms commute with it and therefore the rotated-electron 
occupancy configurations generate for all densities state representations of its global symmetry. 

Since the expression of the operator S involves only the three kinetic operators Tq, T-|_i, and r_i given in Eq. 
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the electron - rotated-electron unitary operator V — V preserves the occurrence of nearest hopping only for rotated 

electrons alike for electrons and commutes with the momentum operator P and, therefore, according to Eq. ([3]) the 

latter operator is such that P = P and thus has the same expression in terms of electron and rotated-electron creation 
and annihilation operators, 

^= E EH.^^,^= E EH.s,,.- (5) 

(T=t,4. k cr=t,i k 

In addition, according to the studies of Ref.^ the unitary operator V commutes with the three generators of the 
spin SU{2) symmetry and three generators of the 77-spin SU(2) symmetry. Hence, such generators also have the same 
expression in terms of electron and rotated-electron creation and annihilation operators and read, 

si = T.^'-'^^^Aa = T.^''''^.A^' ^. = E^"*''^^.tc^.i = E^"^''''^^.t5^,i, 
^ E 4, 4 ""f^ .t = E 4, ,1 ^ ,t ; Ss^Y. 4, ,t ,4 = E 4, ,t i > (6) 

where for the model on the square (and ID) lattice the vector 7? has Cartesian components 7? = [vr, tt] (and component 

In turn, alike the Hamiltonian of Eq. ([1]), the generator Sc of the charge independent U{1) symmetry does not 
commute with the unitary operator V . On the contrary of the Hamiltonian, that operator has a complicated expression 
in terms of electron creation and annihilation operators and a simple expression in terms of rotated-electron creation 
and annihilation operators^, which reads, 

-Sc = 2 E E ^f^^'^ (1 ~ 'V,,-a) ■ (7) 
j=l <T=t,; 

It follows that Sc = ScV where Sc = ■^Q and the operator Q is given in Eq. (P). The six operators provided in 
Eq. © plus that given in Eq. ^ are the seven generators of the group [S'0(4) x U{l)]/Z2 = 50(3) x SO{3) x U{1) 
associated with the global symmetry of the Hamiltonian (jlj. 

That for U/At > the eigenvalue Sc of the generator ([7]) is a good quantum number implies that the rotated-electron 
doubly occupancy Dc = [N/2 — Sc] is also a good quantum number. Alike for the more general case considered in 
Ref.^, following the unitary character of the operator V = V associated with the rotated electrons as constructed 
in this paper one can either consider that V HV^ = H + [H, S] + ... is the Hubbard model written in terms of 
rotated-electron operators or another Hamiltonian with an involved expression and whose operators ct „ and cp a 

rj<7 J 

correspond to electrons. The studies of Rei^ refer to any unitary operator V associated with the general type of 
unitary transformations considered in Ref.^^. In turn, the operator V used here corresponds to a suitable set of 4^" 
U/At — ^ 00 energy eigenstates {|^oo)} such that for U/4t > the V tower states I'iu/At) — V^^l^oo) are energy and 
momentum eigenstates. In the more general case considered in that reference the 4^° states V'^l^'oo) are eigenstates 
of the generator of the charge U{1) symmetry and of the 77-spin and spin square operators of the SU{2) symmetries so 
that the values of Sc, Sjj, and Ss are well defined yet such states are not in general energy and momentum eigenstates. 

Following the procedures of Ref . ^ ^ , the Hamiltonian V HV^ — H+[H, S* ] -I- . . . is built up by use of the conservation 
of singly occupancy 2Sc — (Q) by eliminating terms in the t > Hubbard Hamiltonian so that Sc is an eigenvalue 
of the operator ([7]). According to the studies of that reference this can be done to all orders of 4:t/U provided that 
U/4t ^ 0. The original lattice remains invariant under the electron - rotated-electron unitary transformation and 
the energy eigenstate description introduced in this paper has been constructed to inherently each site of that lattice 
being unoccupied, doubly occupied, singly occupied by a spin-up rotated electron, or singly occupied by a spin-down 
rotated electron. 
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III. THREE ELEMENTARY QUANTUM OBJECTS AND CORRESPONDING c, ry-SPIN, AND SPIN 

EFFECTIVE LATTICES 

A. Elementary quantum objects and their operators 

The rotated electrons are not the uhimate objects whose occupancy configurations that generate the energy eigen- 
states of the model ([1]) in the one- and two-electron subspace have simplest expressions. There emerge from the 
electron - rotated-electron unitary transformation considered here three basic objects: the spin-less and 77-spin-less c 
fermions, the spin-1/2 spinous, and the 77-spin-l/2 7y-spinons. Such objects are a generalization for U/At > of those 
introduced in RefJ^ for U /At ^ 1, as discussed below. In Appendix B it is confirmed that for the ID Hubbard model 
and all U/At values suitable occupancy configurations of such objects generate a complete set of energy eigenstates, 
which span the whole Hilbert space. In this paper and in Ref.-''' evidence is provided that suitable c and si fermion 
momentum occupancy configurations generate the exact ground states and the excited energy eigenstates that span 
the one- and two-electron subspace defined below in Section V. 

Here we profit from the result of Appendix A that full information about the quantum problem can be achieved by 
defining it in the LWS-subspace. Within the corresponding LWS representation the c fermion creation operator can 
be expressed in terms of the rotated-electron operators and electron operators as follows, 

nr, = 4, ~cr. ; 4, ,.^y^ 4, y-^ ~cr.,^^ ^r. ^ . (8) 

where we have introduced the corresponding c fermion momentum-dependent operators as well, e*'^ ''^ is ±1 depending 
on which sublattice site Vj is on, and is the uniquely defined electron - rotated-electron unitary operator introduced 
in this paper. (For the ID lattice that phase factor can be written as (—1)^ .) The expression in terms of electron 
operators involves the electron - rotated-electron unitary operator V as defined in this paper. The c momentum band 
is studied in RefjA and has the same shape and momentum area as the electronic first-Brillouin zone. 

The three spinon local operators and three ?7-spinon local operators such that ^ = ±, 2; are given by. 



^pj,c q?^ ; = (1 - n-fj.c) qp^, l ^±,z. (9) 



Here qp = q^, ± zgH are the rotated quasi-spin operators whose Cartesian coordinates x,y,z are often denoted in 
this paper by xi,X2,x^, respectively, and, 

is the c fermion local density operator. 

In terms of rotated-electron creation and annihilation operators the rotated quasi-spin operators read, 

<li=i4„t-^'-''^r,^)cr„i; q-^^iqiy-, (11) 

Since the electron - rotated-electron transformation generated by the operator V is unitary, the operators ct ^ and 

Crj,a have the same anticommutation relations as ct, ^ and 0^7^,0-, respectively. Straightforward manipulations based 
on Eqs. (|5|)- ((m) then lead to the following algebra for the c fermion operators, 

{4,c . fry.c} = ; {fl, , 4„ J = {/.-, ,c , fry A - , (12) 

c fermion operators and rotated quasi-spin operators, 

[4,c-4] = [^.,c,4]=0, (13) 

and rotated quasi-spin operators, 

[qp. , qp., ] = i ^j,r ^ ^pp'p" q^. ; p = x,y,z, (14) 
p" 
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{q+,qZ^} = l, {g±,g±}^0, (15) 

Hence the operators anticommute on the same site and commute on different sites. 

_ 

The relations provided in Eqs. ((T^ - P^ confirm that the c fermions associated with the global U{1) symmetry are 
7y-spinless and spinless fermionic objects and the spinons and 7y-spinons are spin-1/2 and 77-spin-l/2 objects whose 
local operators obey the usual corresponding spin and 77-spin SU{2) algebras, respectively. 

We can now fulfill task (i) suggested in Subsection II- A: Writing the rotated-electron creation and/or annihilation 
operators in the expression of any one- or two-electron operator O of Eq. ([3]) in terms of the operators of the objects 
whose simple occupancy configurations generate the energy eigenstates that span the one- and two-electron subspace 
of the model on the square lattice and the whole Hilbert space of the ID model. For the LWS-subspace considered 
here this is simply achieved by inverting the relations given in Eqs. ([5]) and (|lip with the result, 

4,,t = 4,c [I + 4) + e*''^ ^..^ {I - 4) ; 4,4 = (4,c - e*^"^ He) ■ (17) 

For U /At ^ 00 the rotated electrons become electrons and Eqs. ((8))- (IT7|) are equivalent to Eqs. (l)-(3) of Ref^i^ with 
the rotated-clcctron operators replaced by the corresponding electron operators and the c fermion creation operator 
ft ^ replaced by the quasicharge annihilation operator Cr- Therefore, in that limit the c fermions are the "holes" 
of the quasicharge particles of that reference and the spinons and 77-spinons are associated with the local spin and 
pseudospin operators, respectively. 

Since the transformation considered in Ref.^^ does not introduce Hilbert-space constraints, suitable occupancy 
configurations of the objects associated with the local quasicharge, spin, and pseudospin operators generate a complete 
set of states. However, only in the U/At — > 00 limit simple occupancy configurations of the corresponding quasicharge, 
spin, and pseudospin objects generate the energy eigenstates that span the one- and two-electron subspace. In turn, 
for the model on the square lattice the same occupancy configurations of the above corresponding finite-C/ /At objects 
generate the energy eigenstates belonging to the same V tower. 

Our next main goal is the fulfillment of the task (ii) also suggested in Subsection II-A: Finding the suitable occupancy 
configurations of the c fermions, spinons, and ry-spinons that generate energy eigenstates of the one- and two-electron 
subspace. The study of the relationship of the transformation laws of such objects under the electron - rotated-electron 
unitary transformation to the state representations of the group 5*0 (3) x SO {3) x U{1) associated with the model 
global symmetry provides valuable and useful yet partial information about that complex problem. Here we start 
a preliminary analysis of the degrees of freedom of the rotated-electron occupancy configurations that generate a 
complete set of Srj, S^, Ss, S^, and Sc eigenstates and introduce the theory vacua. Such states correspond to 
representations of the model global SO{3) x SO{3) x U{1) symmetry. As discussed in Section IV, they are momentum 
eigenstates of the square-lattice model in the whole Hilbert space. For the ID Hubbard model the states of such a 
complete set are shown in Appendix B to be both momentum and energy eigenstates. 

B. Interplay of the global symmetry with the transformation laws under the operator V: three basic 

eflFective lattices and the theory vacua 

The rotated-electron occupancy configurations associated with the state representations of the global 5*0(3) x 
50(3) X U{1) = [SU{2) X SU{2) x C/(l)]/Z| symmetry are weU defined for U/At > 0. It is found in Subsection IV-D 
that the 77-spin and spin state representations correspond to independent rotated-electron occupancy configurations 
of [N^ — 2Sc] sites and 2Sc sites, respectively, whereas the t/(l) symmetry refers to the relative positions of the sites 
involved in each of these two types of configurations. In the present Na ^ 1 limit it is useful to introduce the numbers 
Na^ and Na, such that iV^^ = [N^ - 25c] and N^^ = 25*^, 

Na,^ = {N^-2S,)'/^; iV,^, = (25,)!/^ ; iVf = < + TV,^ . (18) 

Here N//'^ and N^^ are integer numbers that below are identified with the number of sites of a ?7-spin and spin effective 

lattice, respectively. For the D = 2 square lattice the number is chosen so that the number Na of sites in a row or 
column is an integer. However, the designations N/^^ and N/^^ do not imply that the corresponding numbers Na^ and 
Na, are integers. In general they are not integers but for finite values of x and (1 — x) the closest integer numbers to 
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Na^ and Na^ give the mean value of the number of sites in a row or column of the vy-spin and spin effective lattices, 
respectively. 

The degrees of freedom of the rotated-electron occupancy configurations of each of the N^^ = [A^^f — 2Sc] and 
iV^ = 2Sc sites of the original lattice that generate the Sjj, S^, Sg, S^, and Sc eigenstates studied in Section IV, 
which correspond to state representations of the model global 5*0(3) x S0{3) x U{1) symmetry, naturally separate 
as follows: 

i) The occupancy configurations of the c fermions associated with the operators fl ^ of Eq. where j — 1, 

correspond to the state representations of the global U{1) symmetry found in Ref.^. Such c fermions live on the c 
effective lattice, which is identical to the original lattice. Its occupancies are related to those of the rotated electrons, 
the number of c fermion occupied and unoccupied sites being given by Nc = = 2Sc and Nj^ = = [N^ — 2Sc], 
respectively. Indeed, the c fermions occupy the sites singly occupied by the rotated electrons whereas the rotated- 
electron doubly-occupied and unoccupied sites are those unoccupied by the c fermions. Hence the c fermion occupancy 
configurations describe the relative positions in the original lattice of the N^^ = [N^ — 2Sc] sites of the rj-spin effective 
lattice and — 2Sc sites of the spin effective lattice. 

ii) In turn, the remaining degrees of freedom of rotated-electron occupancies of the iV^ = [N^ — 2S^ and = 2Sc 
original lattice sites correspond to the occupancy configurations associated with the 77-spin SU (2) symmetry and spin 
SU{2) symmetry representations, respectively. The occupancy configurations of the N^^ — [N^ — 25c] sites of the 
77-spin effective lattice and iV^ = 25c sites of the spin effective lattice are independent. They refer to the operators 
p^. of Eq. ([5]), which act only onto the iV^^ = \Nf^ — 25c] sites of the 77-spin effective lattice, and to the operators 
s^. given in the same equation, which act onto the A^^f^ = 25'c sites of the spin effective lattice, respectively. This is 
assured by the operators (1 — np.^c) and n^-^c in their expressions provided in that equation, which play the role of 
projectors onto the 77-spin and spin effective lattice, respectively. 

If follows that for C//4t > the degrees of freedom of each rotated-electron singly occupied site separate into a 
spin-less c fermion carrying the electronic charge and a spin-down or spin-up spinon. Furthermore, the degrees of 
freedom of each rotated-electron doubly-occupied or unoccupied site separate into a 77-spin-less c fermion hole and a 
77-spin-down or 77-spin-up 77-spinon, respectively. The 7;-spin-down or 77-spin-up 77-spinon refers to the 7y-spin degrees 
of freedom of a rotated-electron doubly-occupied or unoccupied site, respectively, of the original lattice. Above and 
in the remaining of this paper we often call c fermion hole an unoccupied site of the c effective lattice. 

A key point of our description is that for [//4i > its quantum objects correspond to rotated-electron configurations 
whose spin-down and spin- up singly occupancy, doubly occupancy, and no occupancy refer to good quantum numbers. 
This is in contrast to descriptions in terms of electronic configurations, whose singly occupancy refers to a good 
quantum number for U/At^ 1 only->^i^>2^. 

The transformation laws of the 77-spinons (and spinous) under the electron - rotated-electron unitary transformation 
associated with the operator V play a major role in the description of the 7/-spin (and spin) state representations in 
terms of occupancy configurations of the iV^ = — 25c] sites of the 77-spin effective lattice (and Nj^ — 25c sites 
of the spin effective lattice). Importantly, a well-defined number of 77-spinons (and spinous) remain invariant under 
the unitary transformation generated by V. Those are called independent ±1/2 77-spinons (and independent ±1/2 
spinous) and as further discussed below play the role of unoccupied sites of the 77-spin (and spin) effective lattice. 
The values of the numbers ±1/2 of independent ±1/2 77-spinons and ±1/2 of independent ±1/2 spinous are fully 
controlled by those of the 77-spin Sn and 77-spin projection 5"^ = —x and spin Ss and spin projection 5| = —771 , 
respectively, and are given by, 

Lq = [Xq _i/2 + La^+1/2] = 25a ; La^±i/2 = [SaT S^]; a = ri,s. (19) 

The invariance of such independent 77-spinons (and spinons) stems from the off diagonal generators of the 77-spin 
(and spin) algebra, which flip their 77-spin (and spin), commuting with the unitary operator V and hence having for 
U/At > the same expressions in terms of electron and rotated-electron operators, as given in Eq. ([5]). 

It follows that the number of sites of the 77-spin (a = 77) and spin (a — s) effective lattice can be written as, 

iVf^ = [25c + 2a] ; a = 77 , s . (20) 

The occupancy configurations of 25*0 sites out of N^^ have for U/At > the same form in terms of electrons and 
rotated electrons. However, that refers to the 77-spin or spin effective lattices and hence only to the 77-spin (a = 77) or 
spin (a = s) degrees of freedom, respectively, of the rotated-electron occupancies of the original lattice. Indeed, for 
finite values of U /At and spin density m < {1 — x) the c fermion occupancy configurations are not invariant under the 
electron - rotated-electron unitary transformation. 
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Nevertheless, the two sets of 25",, and 2C^ (and 2Ss and 2Cs) sites of the 77-spin (and spin) effective lattice correspond 
to two weU-defined sets of 2S',j and 2C,, (and 2Ss and 2Cs) sites of the original lattice. Out of the set of 25,, = 
[L,, _i/2 + (and 2Ss — [^5.-1/2 + is.+i/2]) sites of the original lattice, Ln.-i/2 are doubly occupied and 

Lr). +1/2 unoccupied both by electrons and rotated electrons (and Ls _i/2 and _|_i/2 are singly occupied both by 
spin-down and spin-up, respectively, electrons and rotated electrons.) Out of the 2C,, (and 2Cs) sites of the original 
lattice left over, the corresponding C,, (and Cs) sites are unoccupied by rotated electrons (and singly occupied by 
spin-up rotated electrons) and the remaining C,, (and Cs) sites are doubly occupied by rotated electrons (and singly 
occupied by spin-down rotated electrons). In turn, in terms of electrons these [2C,, 4- 2Cs] sites of the original lattice 
have for finite U /At values very involved occupancies such that singly and doubly occupancy are not good quantum 
numbers and thus are not conserved. 

The site numbers C,, > and C<; > are good quantum numbers given by, 

C, = [<y2-5,] = [iVi^/2-5,-5,]; Cs = [nZ/2-Ss]^[S,-Ss]. (21) 

Hence their values are fully determined by those of the eigenvalue S^. of the global U{1) symmetry generator and 
r/-spin Sr, or spin Ss, respectively, so that Cn and Cs are not independent quantum numbers. 

The physics behind the U{1) symmetry found in Ref.^ includes that brought about by the rotated-electron occupancy 
configurations of the [2C,, -I- 2Cs\ sites of Eq. (PT|) . The use of the corresponding model global 50(3) x 50(3) x U{1) 
symmetry confirms that the integer numbers 2Cc, > are always even and that application onto Sa — ^ states of the 
off-diagonal generators of the 77 spin (a = 77) or spin (a = s) algebra provided in Eq. ^ gives zero. For such states 

On the other hand, application of these generators onto 25q > states flips the 77-spin (a = rf) or spin (a = s) of 
an independent 7y-spinon {a = rf) or independent spinon (a = s) but leaves invariant the rotated-electron occupancy 
configurations of the above considered 2Ca sites. It follows that such 20,, (and 20s) sites refer to Ty-spin-singlet (and 
spin-singlet) configurations involving Cn (and Cs) ~l/2 ?7-spinons (and —1/2 spinous) and an equal number of -1-1/2 
77-spinons (and +1/2 spinous). 

One then concludes that concerning the rotated-electron occupancies of the original lattice, besides the [25,, + 25s] 
sites whose occupancies are the same in terms of electrons and rotated electrons, there are O,, > sites doubly occupied 
by rotated electrons, O,, > sites unoccupied by rotated electrons, Os > sites singly occupied by spin-down rotated 
electrons, and Os > sites singly occupied by spin-up rotated electrons. 

Let Ma = Ma -1/2 + -^a,+i/2 denote the total number of ry-spinons {a = 77) and spinons (a = s). It is given 
by Ma = iVi^ and hence M„ = Nj^^ = [A^f - 25J and Ms = iVj^ = 25^. Out of the M„ = iV^^ = [iVf - 25d 

77-spinons (and Ms — A^^^ — 25c spinons), — 25,, (and Ls — 25s) are invariant under the unitary operator V and 
20,, = [Mri — 25,,] (and 20s = [Mg — 25s]) are not invariant under that unitary operator. Also the c fermions are not 
invariant under V . Indeed, the generators of the occupancy configurations of the latter objects are for finite values of 
U/At and m < [1 — x) different when expressed in terms of rotated-electron and electron operators, respectively. The 
number of 7/-spinons {a — rj) and spinons (a = s) can then be written as. 

Ma. ±1/2 = [La. ±1/2 + Ca] ] Ma = N^^ = [25„ +20^.]; a = 77 , s . (22) 

The 77-spinon and spinon operator algebra refers to well-defined subspaces spanned by states whose number of each 
of these basic objects is constant and given by Mjj = N^^ = [N^ — 25c] and Ms = N/f'^ — 25c, respectively. Hence 
in such subspaces the number 25c of rotated-electron singly occupied sites and the numbers and A"^ of sites of 
the 77-spin and spin effective lattices, respectively, are constant. For hole concentrations < a; < 1 and maximum 
spin density m = {1 — x) reached at a critical magnetic field He parallel to the square-lattice plane for D = 2 and 
pointing along the chain for D = 1 the c fermion operators are invariant under the electron - rotated-electron unitary 
transformation and there is a fully polarized vacuum |0,,s), which remains invariant under such a transformation. It 
reads, 

|0,s) = |0,; N^J X jOs; A^i^) x |05c; 25c) , (23) 

where the 77-spin SU{2) vacuum jO,,; A"^) associated with N^^ independent -fl/2 ry-spinons, the spin SU{2) vacuum 

jOs; A"^) with independent -1-1/2 spinons, and the c U{1) vacuum |G5c;25c) with A'c — 25c c fermions remain 
invariant under the electron - rotated-electron unitary transformation. The explicit expression of the state |05c; 25c) 
used below in Subsection IV-E is Yiqflj^'^c', 0) where the vacuum |G5c; 0) corresponds to the electron and rotated- 
electron vacuum of form (|23l) referring to A"^ — NJ^ and NJ^^ — 25c — 0. Only for a ttt, = (1 — x) fully polarized 
state are the state |05c; 25c) and the corresponding A'c = 25c fermions invariant under the electron - rotated-electron 
unitary transformation for U/At > 0. For the vacuum \Orj;N^) (and |Os;A^i^)) the M,, = independent +1/2 
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Ty-spinons refer to N^^ sites of the original lattice unoccupied by rotated electrons (and the Mg = iV^^ independent 
+1/2 spinous to the spins of spin- up rotated electrons that singly occupy sites of such a lattice). At maximum 
spin density m = [1 — x) the c fermions are the non-interacting spinless fermions that describe the charge degrees 
of freedom of the electrons of the fully polarized ground state. At that spin density there are no electron doubly 
occupied sites and the quantum problem is non-interacting for U /At > 0. 

The studies of Ref.— reveal that the a; = and m = absolute ground state such that A^^f^^ = 0, = A^ — 2Nsi, 
and 2Sc = N/^ where A^^i is the number of si fermions considered below is also invariant under the electron - rotated- 
electron transformation. Such an invariance is related to the Mott-Hubbard insulator physics. This is consistent with 
the Sc — A^/2 vacuum referring to a spin-up fully polarized state, whose numbers are also A"^ = 0, A"^ = A^^f, and 
2Sc = N/^ and correspond to a particular case of the general state (|23p . being invariant under that transformation. 

Following the result of Appendix A that full information about the quantum problem described by the Hamiltonian 
([1]) can be achieved by defining it in the LWS-subspace spanned by states such that L^ _i/2 = and iQ.+i/2 = 2Sa 
for a = rj ,s, in the ensuing section we confirm that within the description introduced in this paper, out of the 
A'jP — [2Sa + 2Cq] sites of the 77-spin [a — rf) and spin (a — s) effective lattice, the 2Sa lattices whose occupancies 

are invariant under V play the role of unoccupied sites, whereas the remaining 2Ca sites play the role of occupied sites. 
This is a natural consequence of the ?7-spin SU{2) vacuum |0^; A^^) (and spin SU{2) vacuum |0s; A^^)) being for all 
U /At and m values invariant under the electron - rotated-electron unitary transformation and such that A'^ = 25*^ 
so that 2C^ = (and A^^ = 25^ so that 2Cs = 0). 



IV. COMPOSITE au FERMIONS, ai^ BOND PARTICLES, AND CORRESPONDING ai^ EFFECTIVE 

LATTICES 



In this section we consider a uniquely chosen complete set of S*^, S^, Ss, S^, 5c, and momentum eigenstates that 
within the present description for the Hubbard model on the ID or square lattices is generated by c fermion, ?7-spinon, 
and spinon occupancy configurations. As discussed in Appendix B, for the ID Hubbard model with U/At > such 
states are a complete set of energy eigenstates. For the model on the square lattice such states are in some particular 
cases energy eigenstates, including ground states and the states that span the one- and two-electron subspace as 
defined in Section V. Here we access partial yet valuable information about the Ty-spinon and spinon occupancy 
configurations that generate the yy-spin and spin degrees of freedom of such Sri, S^, Ss, S^, Sc, and momentum 
eigenstates from the suitable use of the transformations laws of the above objects under the electron - rotated-electron 
unitary transformation. 

Specifically, within the general description introduced in this paper both for the Hubbard model on the square and 
ID lattice, the complete set of Sri, S^, Sg, 5*1, Sc, and momentum eigenstates can in the A^^f — > 00 limit be generated 
by occupancy configurations of c fermions, several branches of composite av fermions labeled by the indices a — t], s 
and v = 1,2,..., and a well-defined number of independent 7]-spinons and independent spinous. The latter objects 
are invariant under the above unitary transformation and thus have a non-interacting character for U/At > 0. The 
composite av fermions and corresponding av bond particles involve 2i/ //-spinons (a = 77) or 21^ spinons {a — s) where 

= 1, 2, ... is the number of ?7-spin-neutral ?7-spinon (and spin-neutral spinon) pairs. 



A. The composite ai' bond particles and ai' fermions 

1. The av fermion operators and au translation generators 

Alike for the si bond particles further studied in RefSfiii^, one can introduce within the A^^P — > 00 limit suitable 
creation operators gt, for the aiy bond particles. Provided that {1 — x) > for ai/ ~ si and Sa/N/^ > for the 
remaining av branches, the av bond-particle operators have been constructed to inherently upon acting onto their 
av effective lattice introduced below in Subsection IV-D anticommuting on the same site and commuting on different 
sites, so that they are hard-core like and can be transformed onto fermionic operators whose general expression reads. 



^j' 2 2 



4. = —L= V e+*«^ -^^.' /t , , (1 - a;) > for a,, = si and 5„/Af > for ai. ^ si . (24) 
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Here (jfj^av is the Jordan- Wigner phase^" operator, the indices j' and j refer to sites of the av effective lattice, and 
/I are the corresponding momentum-dependent av fermion operators. The number N^^^ of discrete momentum 
values of the av momentum band equals that of sites of the effective lattice derived below in Subsection IV-D. 
In turn, if for ai^ ^ si branches one has that Sa — and N^^/N^ <^ 1, provided that Na,y' ~ for all remaining 
ai'' branches with a number of 77-spinon or spinon pairs v' > v one finds below that N^^^^ — Na,y, the phase operator 

4'j,ai' may be replaced by the average phase (ftj^af ~ av — '^j'^j 't'j' ■/j,av, and the momenta of the operators /I. 
are given by qj- « 0. In the latter case all sites of the av effective lattice are occupied and the av momentum band 
is full. If the value of Nav is finite one has = N^v discrete momentum values qj ~ Q distributed around zero 

momentum, whose Cartesian components momentum spacing is 2'k/L. A case of interest is when Nau = 1 so that 
the av effective lattice has a single site and the corresponding av band a single discrete momentum value, q — Q. In 
that case 6„u = <Pi = so that fl = at and ft = ft where q = 0. 

As further discussed in Subsection IV-E, the operators ft act onto subspaces with constant values for the set of 
numbers Sa, N^v, and {N^u'} for v' > v and equivalently of the numbers Sc, N^u^ and {N^v'} for all v' 7^ v. Such 
subspaces are spanned by mutually neutral states, that is states with constant values for the numbers of av fermions 
and av fermion holes, respectively. Hence such states can be transformed into each other by av band particle-hole 
processes. The phases (pj^av given in Eq. ([24|) are associated with an effective vector potential^Sii^, 



Bauirj) = '^r, ><■ ^a,^{rj) = ^o^nf^,,a,^^{rj' - rj)ex^; $0 = 1, (25) 

where e^-i is the unit vector perpendicular to the plane and we use units such that the fictitious magnetic flux quantum 
is given by $0 = 1- It follows from the form of the effective vector potential Aa^^fj) that the present description leads 
to the intriguing situation where the av fermions interact via long-range forces while all interactions in the original 
Hamiltonian are onsite. 

The components of the microscopic momenta of the av fermions are eigenvalues of the two (and one for ID) av 
translation generators qa^xi and qavx2 in the presence of the fictitious magnetic field Bav{rj). That seems to imply 
that for the model on the square lattice the components qxi and qx2 of the microscopic momenta q = [qxi^ qx2] refer 
to operators that do not commute. However, such operators commute in the subspaces where the operators ft act 
onto because those are spanned by neutral states'^^. Since [qav xij<laij X2\ = in such subspaces, for the model on the 
square lattice the av fermions carry a microscopic momentum q = [qxiT(lx2\ where the components q^i and qx2 are 

well-defined simultaneously. Importantly, it follows that the momentum operator P of Eq. ([5]) commutes with the 
av generators q^u whose eigenvalues are the av fermion microscopic momenta q. Hence, within our description the 
momentum operator ([5]) can be written as, 

Cs Cn 

P^t + Y. ^s. + + ^ ^^-n .-1/2 ' (26) 

where the c and av translation generators read, 

(f, = ^giV,(g); 'q,, = Y<i^sA?) ^ 9*,,. = [t? - 9l ^r,. (g) • (27) 

"3 9 9 

Here tt is the momentum carried by a —1/2 ry-spinon, the operator .-1/2 counts the number ,-1/2 = [L,-i ,-1/2 + 
Cn] of such objects, that the rjv generators involve [tt ~ q\ instead of q follows from the anti-bounding character of the 
rjv fermions discussed below, and N^iq) and N^uiq) are the momentum distribution- function operators, 

Nc{q) = ftjq,c; NaM^ fia.ho.., (28) 

respectively. 

For the Hubbard model both on the square and ID lattices the Hamiltonian H of Eq. ([1]) and the momentum 
operator P of Eqs. ([5]) and obey the commutation relations, 

[H, P] - [H, Ma ,±1/2] = [H, t] - ; [P, Af„ ,±1/2] = [P, <fj ^ [P, J = , (29) 



14 



for P = 1, Ca and a = i],s. In turn, the set of commutators [H, q^^] vanish for the ID model whereas for the model 
on the square lattice one has in general that , q^^] ^ 0. It follows from Eq. (j26p that the momentum eigenvalues 
P can be expressed as a sum of the filled c and aiy fermion microscopic momenta. 

In the limit Na — >■ oo one has that X]a=j; s for most subspaces so that the set of aiy translation generators 

q^^ of Eq. (|27|) in the presence of the fictitious magnetic fields Bol- of Eq. psp is infinite. That for ID all such 
operators commute both with the Hamiltonian and momentum operator is behind the integrability of the ID Hubbard 
model in that limit. Indeed, such an infinite set of translation generators is equivalent to the infinite set of conservation 
laws, which are known to be behind that model integrabilitj^i^. Such laws are also equivalent to the conservation 
of the set of av fermion numbers {Na^} where a = ry, s, v = 1, Ca and the gnitude is /2 and 

thus Na/2 for ID^. In turn, the Hubbard model on the square lattice is not integrable and consistently the set of av 
translation generators q^u of Eq. (j27p do not commute in general with the Hamiltonian so that the corresponding av 
fermion discrete momentum values qj = [qxi, qx2] are not good quantum numbers. As further discussed in Subsection 
IV-E, that as given in Eq. (j29p for the square-lattice model such av translation generators commute with momentum 
operator is consistent with it being invariant under the electron - rotated-electron unitary transformation and the 
av fermion operators acting onto subspaces spanned by mutually neutral states. In general such subspaces refer to 
constant values of the number of av fermions. 

It turns out that processes within the one- and two-electron subspace as defined in this paper where two si fermions 
are annihilated and one s2 fermion created are also neutral in the sense that for the model on the square lattice the 
si fermion microscopic momenta q = [qxi, qx2\ are for such excitations associated with si translation generators qaixx 
and qsix2 that commute. Moreover, in such a subspace there are no av fermions other than si fermions and none 
or one zero-momentum s2 fermion. The interest of our operator description lies in that for the Hubbard model on 
the square lattice in the neutral subspaces of the one- and two-electron subspace the si translation generators q^i of 
Eq. (I27p in the presence of the fictitious magnetic field Bsi commuting with both the Hamiltonian and momentum 
operator. As justified below, in spite of the lack of the square-lattice model integrability, for that model in such 
a subspace the si translation generators have been constructed to inherently these commutation relations holding. 
Therefore, for the square-lattice quantum liquid corresponding to the Hubbard model on a square lattice in the one- 
and two-electron subspace as defined in this paper the si fermion discrete momentum values q = [qxi,qx2] are good 
quantum numbers and thus are conserved. We emphasize that according to Eq. (j29p the commutator [iJ , q^ vanishes 
so that for the model on the square lattice the c band momenta are good quantum numbers for the whole Hilbert 
space, alike for the ID model. Hence for the model on the square lattice in the one- and two-electron subspace both 
the c band and si band discrete momentum values are good quantum numbers. This result follows from the way that 
our operator description is constructed. However, the shape of the si momentum band and of its boundary as well 
as the form of the si and c fermion energy dispersions remain unsolved problems, which are addressed in Ref.— . 

The advantage of using momentum-dependent c fermion operators /I. ^ and av fermion operators /I. given in 

Eqs. ([8]) and ((24)) . respectively, is then that for the ID model in the whole Hilbert space and the model on the square 
lattice in the one- and two-electron subspace as defined in this paper such microscopic momenta are good quantum 
numbers. That for the model on the square lattice the av translation generators q^u of Eq. (j27p in the presence of 
the fictitious magnetic field Ba^ do not commute in general with the Hamiltonian is consistent with the set of av 
fermion numbers {Nav} which label the 5^, S*^, Ss^ S^, Sc, and momentum eigenstates considered in Subsection IV-E 
not being in general good quantum numbers. 

The related studies of Refs.^'^^ focus mainly on the square- lattice quantum liquid in the one- and two-electron 
subspace. In turn, here we consider as well general A/'-electron subspaces where Af = 1,2,3,4,.... The av fermion 
operators f^^u ^^'^ fq,av act onto subspaces spanned by neutral states. Nevertheless, creation of one av fermion is 
a well-defined process whose generator is the product of an operator, which fulfills small changes in the av effective 
lattice and corresponding av momentum band and the operator f^^^ appropriate to the excited-state subspace. It 
turns out that for the model on the ID and square lattices creation of one rjv fermion onto a [N^ — N] — 2v ground 
state generates for U/it > a, 2z/-electron charge excited energy eigenstate. Moreover, creation of one sv fermion 
onto a [Nf — N^] = 2v ground state generates a 2i/-electron spin excited energy eigenstate. Furthermore and as 
discussed below in Subsection IV-E, the created rjv fermion or sv fermion is invariant under the electron - rotated- 
electron unitary transformation. However, the detailed internal structure of the composite av fermions is in general 
a unsolved complex problem. 

Since the microscopic momenta of the c fermions are good quantum numbers for the model on the square lattice, 
one can define an energy dispersion ec{q}^. For the ID model the av energy dispersions (.av{,q) are well defined as well. 
In turn, for the model on the square lattice the energies tavisi) have for av ^ si no uncertainty for the momentum 
values for which such objects are invariant under the electron - rotated-electron unitary transformation. This is the 
case of the above r\v fermion created onto a —'N\ —2v ground state and sv fermion onto a [iV-|- — — 2v ground 
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state whose excited-state aiy momentum band corresponds to the zero momentum value only, as confirmed below. For 
general ai' fermion momentum values there is some uncertainty in the magnitude of e^^ (q) yet fortunately it refers to 
a well-defined energy window. There is in general a uncertainty in the magnitude of esi{q) as well yet such an energy 
dispersion is well defined for the one- and two-electron subspace^. 

Within our ^ 1 limit on can access general useful information on the energies Cc and Cav without knowing 
in detail the form of the aiy bond-particle operators gt, and corresponding ai' fermion operators /I which for 
a — Tj (and a = s) involve a set of jy-spinon operators pi = (1 — np.^c)(l\i. (and spinon operators si = nf^ cQp.) 
given in Eqs. (P)- pT|) associated with the 77-spin (and spin) neutral superposition of the corresponding 2^ — 2, 4, 6, ... 
7/-spinons (and spinous). (The si bond-particle operators gt, are studied in Refji^.) Since as confirmed below, 
the ground-state occupancy configurations involve only the charge c fermions and composite two-spinon si fermions, 
that problem is addressed for the one- and two-electron subspace in Refsjii^^ concerning the si fermions. The studies 
of such references profit from the c fermions and si fermions being the only objects playing an active role in that 
subspace. In turn, in the following we access general information about such composite objects including for instance 
the anti-bounding and bounding character of the 2i/-77-spinon composite rjv fermions and 2i^-spinon composite sv 
fermions, respectively, and the energy uncertainty of their energy spectrum, which is smaller than or equals the 
corresponding energy bandwidth. For the Hubbard model on the square lattice in the one- and two-electron subspace 
such objects play the role of "quasiparticles" . Due to integrability, for the ID model they have zero-momentum 
forward-scattering interactions only. In turn, for the model on the square lattice the unusual scattering properties of 
the quantum problem are controlled and determined by the inelastic c - si fermion collisions'^^. 

For the limit 3> 1 considered in this paper the problem of the internal degrees of freedom of the composite ai' 
fermions separates from that of their positions in the corresponding av effective lattice defined below and corresponding 
discrete momentum values in the av band, whose number equals that of sites in such a lattice. In this section it is 
confirmed that for the subspaces spanned by states with constant values of 5*,,, Ss, and Sc the momentum occupancy 
configurations of the composite av fermions and c fermions generate the correct number of state configurations of 
the group 5*0(3) x S0{3) x C/(l). Addition of the dimensions of all such subspaces leads to a total number 4^" of 
independent state representations, which indeed is the dimension of the Hilbert space. Furthermore, we initiate the 
study of the c and si effective lattices for the states that span the one- and two-electron subspace as defined in this 
paper. 

For the particular case of the ID Hubbard model, a preliminary version of the general operator description intro- 
duced here for both the model on the square and ID lattices is presented in Refj^. The studies of that reference profit 
from the exact solution of the ID model. No explicit relations to the rotated-electron operators as those given in 
Eqs. (I5))- ([T7| are derived and no relation to the U{1) symmetry contained in the model S0{3) x 5*0(3) x U{1) global 
symmetry is established. Indeed, the eigenvalue of the generator ([T]) of the U{1) symmetry is one half the number of 
rotated-electron singly occupied sites 25c. Such an eigenvalue plays an important role in the present description. It 
determines the values of the number iVc = 25c of c fermions, Mg — 25c of spinous, N^} = [N^ — 25c] of c fermion 
holes, and = [N^ — 2Sc] of ry-spinons. In reference^ the 77-spinons, c fermions, independent 77-spinons, and 
independent spinous are called holons, c pseudoparticles, Yang holons, and HL spinous, respectively, where HL stands 
for Heilmann and Lieb. Moreover, the av pseudoparticles considered in that reference are the av fermions. 

2. More about the states representation in terms of composite av fermion occupancy configurations 

For the model ([Ij on the square lattice the sites of the 77-spin (and spin) effective lattice are distributed along 
well-defined 0x1 and 02:2 directions, defined by the fundamental vectors = a, Cx-^ and = 1 of the original 
lattice where Cx-^ and ex2 are the usual Cartesian unit vectors and a is the lattice constant. The average distance 
between such sites when going through these directions is Oq = L/Na^ — [Na/NaJ\ a where a = rj,s. For ID this is 
also the average distance between the sites of such effective lattices, a,, is the average distance along the 0x1 and/or 
the 0x2 direction between the sites that are either doubly occupied or unoccupied by rotated electrons, as is the 
average distance along that and/or those directions between the sites that are singly occupied by rotated electrons. 

The concept of a 77-spin (and spin) effective lattice is well defined for finite values of the hole concentration x 
(and electronic density n = (1 — x)) and the present 3> 1 limit. The general representation introduced in this 
paper contains full information about the relative positions of the sites of the ?7-spin and spin effective lattices in 
the original lattice. For each rotated-electron occupancy configuration that generates a given momentum eigenstate, 
such an information is stored in the corresponding occupancy configurations of the c fermions in their c effective 
lattice, which is identical to the original lattice. The latter configurations correspond to the state representations 
of the U{1) symmetry in the subspaces spanned by states with constant values of 5c, 5,,, and Ss- Indeed, the sites 
of the 77-spin (and spin) effective lattice have in the original lattice the same real-space coordinates as the sites of 
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the c effective lattice unoccupied (and occupied) by c fermions. Tliat sucli an information is stored in tfie c fermion 
occupancy configurations is consistent with the 77-spinon and spinon occupancy configurations of the 77-spin and spin 
effective lattices, respectively, being independent of each other. The latter occupancy configurations refer to the state 
representations of the SU{2) ?7-spin and spin symmetries, respectively, in the subspaces with constant values of Sc, 
Srj, and Ss- 

For the model on the ID lattice the ?7-spin (and spin) effective lattice is straightforwardly and uniquely obtained by 
considering only the rotated-electron doubly and unoccupied sites (and singly occupied sites) of the original lattice 
and keeping the same site order as for that lattice. In turn, for the model on the square lattice the 77-spin (and spin) 
effective lattice is in the limit x — > 1 (and x — ^ 0) exactly a square lattice. Indeed, in that limit it becomes identical to 
the original lattice. Moreover, it is an excellent approximation to consider that for small electron density n = {1 — x) 
(and small hole concentration x) the 77-spin (and spin) effective lattice is a square lattice. 

Since the 77-spinon, spinon, and c fermion description contains full information about the relative positions of the 
sites of the 77-spin and spin effective lattices in the original lattice, it turns out that within the 3> 1 limit and for 
finite values of the hole concentration x (and electron density n = (1 — x)) the 77-spin (and spin) effective lattice can 
for 2D be represented by a square lattice with lattice constant (and as) given by, 

a-a = = a; a = t] , s . (30) 

The expression of the operator S such that = and V = involves only the three kinetic operators Tq, T+i, 
and T_i given in Eq. Therefore, when expressed in terms of rotated-electron creation and annihilation operators 
the expression of the electron - rotated-electron unitary operator V ^ V = e^^ involves only the corresponding 
three rotated kinetic operators Tq, T+i, and T_i and thus preserves the occurrence of only nearest hopping for 
rotated electrons, alike for electrons. Indeed, the effective rotated-electron hopping between second and third nearest 
neighboring sites is generated by products of such operators so that the corresponding elementary hopping processes 
correspond only to nearest hopping. This confirms that in the square lattice the rotated electrons hop along the 0x1 
or the 0x2 direction between the sites. Hence the same holds for 77-spinons (and spinous) moving in the corresponding 
ry-spin and spin effective lattices, consistently with those being square lattices. 

The 77-spin (and spin) effective lattice includes only the rotated-electron doubly and unoccupied sites (and singly 
occupied sites) of the original lattice. However, within periodic boundary conditions a pair of rotated electrons or 
rotated holes (and a rotated electron) moving once around the chain in ID and through the whole crystal along the 
0x1 or 0x2 directions in the square lattice and doubly occupying (and singly occupying) a site in each step must pass 
an overall distance L. Therefore, consistently with expression = L/Na^ = [Na/Na^] a introduced below where 
a = r],s, the 77-spin (and spin) effective lattice has both for ID and 2D the same length and edge length L, respectively, 
as the original lattice. Furthermore, the requirement that for the 2D case when going through the whole crystal along 
the 0x1 or 0x2 directions a 77 spinon (and spinon) passes an overall distance L is met by a square effective 77-spin 
(and spin) lattice. However, since the number of sites sum-rule [N^^ + N^J = holds, the 77-spin and spin effective 

lattices have in general a number of sites smaller than that of the original lattice, . It follows that their lattice 
constants given in Eq. (pO)) are larger than that of the original lattice, Oq. > a where a = 77, s. 

The average distance between the sites of the 77-spin and spin effective lattices provided in that equation then plays 
the role of lattice constant of such lattices. The validity for ^ 1 of the square spin effective lattice constructed 
in that way is confirmed by the average value (^'j^dj^') relative to any energy eigenstate \^) belonging to a subspace 
with constant number of rotated-electron singly occupied sites of the distance 5d in real space of any of the sites 
of the spin effective lattice from the rotated-electron singly occupied site of the original lattice closest to it vanishing 
in the limit — ?> 00. We recall that the number of rotated-electron singly occupied sites of the original lattice equals 
that of sites of the spin effective lattice. 

It is useful to consider the occupancy configurations of the 25*^ (and 2Ss) "unoccupied sites" and the 2C^ (and 
2Cs) "occupied sites" of the 77-spin (and spin) effective lattice that generate state representations of the model global 
50(3) X 50(3) X U{1) symmetry. For the LWS-subspace considered here the 25,, (and 25^) "unoccupied sites" of such 
a lattice correspond to 25,, independent 4-1/2 77-spinons (and 25^ independent +1/2 spinons) that remain invariant 
under the unitary transformation associated with the operator V. In turn, the O,, +1/2 r7-spinons and O,, —1/2 
77-spinons (and Cs +1/2 spinons and Cg —1/2 spinons) that refer to the 20,, (and 2Cs) "occupied sites" of such a 
lattice do not remain invariant under that unitary transformation. 

From the relation of the transformation laws under the electron - rotated-electron unitary transformation to the 
number of 77-spin (and spin) singlet configurations of the subspaces with constant values of 5c, 5,,, and Ss one finds 
that such configurations can be expressed in terms of occupancies of 77-spin-neutral 27^-77-spinon composite rjv bond 
particles (and spin-neutral 2;^-spinon composite sv bond particles) in suitable and independent 77;^ effective lattices 
(and sv effective lattices). We recall that here v — 1,2,...,0,, (and v = l,2,...,Os) is the number of ?7-spinon pairs 
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(and spinon pairs) of opposite jy-spin projection (and spin projection) in each composite tji/ bond particle (and su 
bond particle). Among such composite aiy bond particles, the two-spinon si bond particles play a major role since 
they are found below to be the only of the composite objects with finite occupancy in the ground states. Such two- 
spinon objects are found in Refj-^ to be related to the resonating-valence-bond pictures considered long ago^^^"— 
for spin-singlet occupancy configurations of ground states. 

The av fermions have the same internal structure as the corresponding av bond particles and live on the av 
effective lattice as well. Hence a tji' fermion (and si^ fermion) is a 77-spin-neutral 2i/-?7-spinon composite (and spin- 
neutral 2^-spinon composite) object. Since within our description the momentum occupancy configurations of the c 
and av fermions generate momentum eingestates, in the following we refer mostly to aiy fermions, yet many of our 
considerations apply as well to the corresponding av bond particles. Indeed, as given in Eq. (1241) the creation operator 
fl of a local av fermion of real-space coordinate Vj differs from the creation operator gt of the corresponding 
av bond particle with the same real-space coordinate by a phase factor e"^^ "", fl, = e^4>i,au ^,1^ Otherwise, such 
objects have the same internal structure. 

Our description is consistent with the Abelian and non-Abelian character of the U{1) symmetry and two SU{2) 
symmetries, respectively, of the model global [SU{2) x SU{2) x C/(l)]/Z| = SO{3) x 50(3) x U{1) symmetry. 
Indeed, the c fermions occupancy configurations refer to the state representations of a U{1) symmetry associated with 
an Abelian group and hence involve a single c effective lattice. In contrast, the ?7-spinon (and spinon) occupancy 
configurations refer to the state representations of a SU{2) symmetry associated with a non-Abelian group so that 
several types of ry-spin (and spin) singlet configurations emerge, each corresponding to a rjv effective lattice (and 
sv effective lattice) occupied by ry-spin-zero (and spin- zero) composite objects of 2v 77-spinons (and 2v spinous). As 
confirmed below, there exists an independent av effective lattice for each av branch where a = rj, s and v ^ 1, Cq. 
Here a refers to 77-spin (a = rj) and spin (a = s) and v to the number of 77-spin-singlet 77-spinon pairs (Cq, = C,,) and 
spin-singlet spinon pairs (Cq, = C^). 

Within chromodynamics the quarks have color but all quark-composite physical particles are color- neutral"*^. Here 
the r/-spinon (and spinous) that are not invariant under the electron - rotated-electron unitary transformation have rj 
spin 1/2 (and spin 1/2) but are part of ?7-spin-neutral (and spin-neutral) 2i/-?7-spinon (and 2i^-spinon) composite t]v 
fermions (and sv fermions). The occurrence of such 77-spin neutral rjv fermions (and spin neutral sv fermions) leads 
to the following expression for the number 2Ca of "occupied sites" of the 7;-spin (a — rj) and spin {a = s) effective 
lattice of Eq. ([^0]) . which refers to the subspaces spanned by states with constant values of Sc, Sri, and Ss, 



Here Nau denotes the number of composite av fermions. Each subspace spanned by states with constant values of 
Sc, Srj, and Ss and hence also with constant values of Crj and Cs can be divided into smaller subspaces spanned by 
states with constant values for the set of numbers {Nr,^} and {Ns^}, which must obey the sum rules of Eq. ([3T|) . 
The numbers {A^j^jy} and {Ns,y} correspond to operators that commute with the momentum operator and operators 
associated with the numbers S*^, S^, Sg, 5*1, and Sc- In turn, for the model on the square lattice such operators do 
not commute with the Hamiltonian so that the numbers {Nrf^} and {Ng^} are not in general good quantum numbers. 
The same applies to the set of av translation generators of Eq. p7| in the presence of the fictitious magnetic field 
Bcti, of Eq. (j25p . which as given in Eq. (j29p commute with the momentum operator yet in general do not commute 
with the Hamiltonian of the model on the square lattice. For ID both such av translation generators commute with 
the Hamiltonian and the set of numbers {A^ai^} are within our limit A^^f — > 00 associated with the set of conservation 
laws behind the model integrability^ and thus are conserved. Fortunately, we confirm below that they are good 
quantum numbers as well for the model on the square lattice in the one- and two-electron subspace and read N^ii, = 0, 
Nsi = [Sc — Ss — 2Ns2], Ns2 = 0, 1, and A^^^ = for > 3, consistently with and the si translation generator q^i of 
Eq. ([?7|) in the presence of the fictitious magnetic field Bsi of Eq. ([^5)) commuting with both the Hamiltonian and 
momentum operator. 

In the following we consider the complete set of S*,,, S^, Ss, Ss, Sc, and momentum eigenstates also labelled by the 
numbers {N^ii,} and {Nsu}- The expressions of the corresponding av fermion operators involve the spinon (a = s) and 
?7-spinon (a = rf) operators of Eqs. (llOp . and (jlip . Both such spinon and ?7-spinon operators and the c fermion 
operators are expressed in terms of the rotated-electron operators given in Eq. ([5]), which are generated from the 
corresponding electron operators by the electron - rotated-electron unitary operator as defined in this paper. Hence 
and as discussed below in subsection IV-E, for the model on the square lattice both the present Sr,, S^, Ss, S^, Sc, 
and momentum eigenstates and the unknown energy eigenstates are generated by application onto the corresponding 
sets of U/At 00 states by the same electron - rotated-electron unitary operator = V^{U/4:t). Indeed, for that 
model the energy eigenstates are a superposition of a set of the former 5*,,, S^, Ss, Ss, Sc, and momentum eigenstates 




a = r] ,s . 



(31) 
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with the same momentum and S,j, S^, Ss, S^, and Sc values. For the ID model both sets of states are identical, due 
to the infinite conservation laws, which for — > oo are associated with the model integrability. 

For av fermion branches with N^^, > finite occupancy in a given state the maximum v value I'max obeys the 
inequality Vmax < Cq.. The absolute maximum value Vmax = Ca is reached for a subspace with one av fermion and 
vanishing occupancies for the remaining av' branches such that v' 7^ v. For the states spanning such a subspace the 
2Ca Ty-spinons (a = rj) or spinous (a = s) other than the independent 77-spinons or independent spinous, respectively, 
are part of a single av fermion, which moves around in a av effective lattice with 2Sa unoccupied sites. For Sa = such 
a single av bond particle has a "Big-Bang" like character in that the ry-spin (a = r]) or spin (a = s) degrees of freedom 
of a half-filling N = state with an equal number /2 = i'max of doubly and unoccupied sites or = 2i'rnax 
sites singly occupied by an equal number /2 = Vmax of spin- up and spin-down electrons, respectively, are described 
by such an object. Such two Big-Bang states have 77-spin and spin lattice site numbers: (i) N^^^ — 2vrnax — ^'^'^ 
= referring to a single 77-spin-neutral ly — /2 and iV^P-jy-spinon composite ryi/ fermion in a rji' band lattice 
with a single vanishing momentum value and an empty c momentum band with = c fermion holes; (ii) 
= and N^^ = 2vmax — corresponding to a single spin- neutral i/ ~ /2 and A'^-spinon composite si/ 

fermion in a band with a single vanishing momentum value and a full c momentum band filled by Nc — c 
fermions. It turns out that such Big Bang av fermions are invariant under the electron - rotated-electron unitary 
transformation. That invariance implies that they are A'^-electron objects and hence the corresponding Big-Bang 
states have vanishing overlap with one- and two-electron excitations: Their creation requires application onto the 
ground state whose occupancy configurations we study below of a suitable TV-electron operator. As a result in part of 
their invariance under the electron - rotated-electron unitary transformation, the A'^-electron excited state generated 
by creation of such objects is an exact energy eigenstate both for the model on the ID and square lattice. 

B. Processes that conserve and do not conserve the number of sites of the 77-spin and spin effective lattices 

The vacuum of the theory corresponds to independent +1/2 77-spinons and independent -1-1/2 spinous 
so that as mentioned above such objects play the role of "unoccupied sites" of the 77-spin and spin effective lattices, 
respectively. Consistently, the latter objects have vanishing energy and momentum and relative to that vacuum their 
r7-spin and spin fiip processes correspond to "creation" processes of independent —1/2 77-spinons and independent 

— 1/2 spinous, respectively. 

Furthermore, creation of a local 771^ fermion involves the replacement of 27^ independent +1/2 77-spinons corre- 
sponding to 27/ sites of the 77-spin effective lattice by v —1/2 77-spinons and 1/ +1/2 77-spinons in the 77-spin-neutral 
configuration associated with such a 771/ fermion. Also creation of a local fermion involves the replacement of 2;/ 
independent +1/2 spinous referring to 27/ sites of the spin effective lattice by 1/ —1/2 spinous and v +1/2 spinous in 
the spin-neutral configuration associated with such a sv fermion. Finally, creation (and annihilation) of a c fermion 
involves annihilation (and creation) of a c fermion hole. The latter process involves removal (and addition) of one 
site from (and to) the 77-spin effective lattice and addition (and removal) of one site to (and from) the spin effective 
lattice. 

The momentum eigenstates considered in this paper are generated by microscopic momentum occupancy config- 
urations of c and ai/ fermions. Such configurations can be expressed as a superposition of local rotated-electron 
occupancies in the original lattice. The degrees of freedom of the rotated-electron occupancy of a given site of the 
original lattice are always of two types. For the sites singly occupied by a rotated electron of spin projection —1/2 or 
+1/2, those are a spinless c fermion of charge — e associated with the U{1) symmetry and a spinon of spin projection 

— 1/2 or +1/2 associated with the spin SU{2) symmetry, respectively. In turn, for the sites doubly occupied or unoc- 
cupied by rotated electrons, those are a ?7-spinless c fermion hole associated with the U{1) symmetry and a 77-spinon 
of ?7-spin projection —1/2 or +1/2 associated with the 77-spin SU{2) symmetry, respectively. 

It follows that the 2v sites of the spin effective lattice occupied by one local sv fermion correspond to the spin 
degrees of freedom of 2v rotated-electron singly occupied sites whose charge degrees of freedom are described by 2y 
occupied sites of the c fermion lattice. Therefore, these 2v rotated-electron singly occupied sites are described both 
by the local sv fermion and 2v local c fermions. 

Moreover, the 2z/ sites of the 77-spin effective lattice occupied by one local rjv fermion correspond to the 77-spin 
degrees of freedom of v rotated-electron doubly occupied sites and v rotated-electron unoccupied sites whose U{1) 
symmetry degrees of freedom are described by 21/ unoccupied sites of the c fermion lattice. As a result, the v rotated- 
electron doubly occupied sites and v rotated-electron unoccupied sites are described both by the local 771/ fermion and 
2v local c fermion holes (c effective lattice unoccupied sites). 

In summary, local sv fermions share 2v sites of the original lattice with 2v local c fermions and local rjv fermions 
share 2v sites of that lattice with 2v local c fermion holes. Any pair of local av and a'v' fermions always refer to 
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two different sets of 2i/ and 2z/' sites, respectively, of the original lattice. Creation of local rjv and sv fermions are 
processes that conserve the total number of ?7-spinons and spinons, respectively, and hence also conserve the number of 
c fermions and c fermion holes. In contrast, creation (and annihilation) of one local c fermion is a process that involves 
addition (and removal) of one site to (and from) the spin effective lattice and removal (and addition) of one site from 
(and to) the ?7-spin effective lattice. Therefore, the spin and ?7-spin effective lattices are exotic, since the number of 
their sites iV^ = 2Sc and N^^ = [N^ — 2Sc], respectively, changes by ±1 and =Fl upon creation/annihilation of one 
c fermion. Indeed, such processes change the eigenvalue Sc of the generator ([7]) of the global U{1) symmetry. A 
subspace with constant Sc value and hence constant iV^ — [N^ — 2Sc] and iV^ = 2Sc values as well is associated 
with a well-defined vacuum |0,,s) given in Eq. ([23]). An excitation involving a change of such values drives the system 
into a new subspace referring to a different vacuum |0,,s)- In turn, 77-spinon and spinon creation and annihilation 
processes refer to excitations within the same quantum- liquid subspace associated with a well-defined vacuum |0^s)- 

It follows that from the point of view of the ry-spin and spin degrees of freedom, c fermion creation and annihilation 
processes correspond to a change of quantum system. Indeed, the ?7-spin and spin lattices and corresponding number 
of sites change along with the quantum-system vacuum |0,,s) of Eq. (j23p . Therefore, within the T^-spinon and spinon 
representation there is a different quantum system for each eigenvalue Sc of the generator ([7]) of the global U{1) 
symmetry. However, from the point of view of the degrees of freedom associated with the latter symmetry, the 
model ([ij corresponds to a single quantum system and the local c fermions live on a lattice whose number of sites 
= [N^ + N^] is constant, alike that of the electrons and rotated electrons. Consistently, the c effective lattice 
is identical to the original lattice. This property follows from the invariance of the latter lattice under the electron - 
rotated-electron unitary transformation, which is related to the invariance of the momentum operator of Eq. ([5]). 

Creation of local -qv (and sv) fermions always involves virtual processes where 2i/ independent -1-1/2 ?7-spinons (and 
2i> independent +1/2 spinons) are replaced by the ry- spin- singlet (and spin-singlet) 2^'-site occupancy configurations 
of the local rji/ fermions (and sv fermions) in the 77-spin (and spin) effective lattice. For instance, a given process for 
which two local si fermions of the initial state are replaced by one local s2 fermion in the final state is divided into two 
virtual processes: first, two local si fermions are annihilated, i.e. the four sites of the spin effective lattice occupied 
in the initial state by the local si fermions are upon two spin-flip processes occupied in an intermediate virtual state 
by four independent -1-1/2 spinons (annihilation of two si bond particles); second, one local s2 fermion is created on 
such four sites. That involves two opposite spin-flip processes and rearrangement of the spinons associated with the 
creation of the local s2 fermion spin-neutral four-spinon occupancy configurations in the spin effective lattice. 

Concerning creation (and annihilation) of one local c fermion, the corresponding overall excitation always involves 
a virtual process for which the site of the T^-spin effective lattice removed (and added) by such an elementary process 
is occupied in the initial (and final) state by an independent -1-1/2 77-spinon. Also the site of the spin effective lattice 
added (and removed) by such an elementary process is occupied in the final (and initial) state by an independent 
-fl/2 spinon. 

If as occurs for one-electron addition (and removal), creation (and annihilation) of a local c fermion involves creation 
(and annihilation) of a local si fermion, the overall process is divided into two virtual processes. For instance, 
complementarily to creation of the local c fermion in its effective lattice, the virtual processes occurring in the spin 
effective lattice are: First, a site occupied by an independent -1-1/2 spinon is added to that lattice; Second, a local si 
fermion is created, one of the two independent -1-1/2 spinons involved in the final-state two-site si bond configuration 
being that located on the site added to the spin effective lattice. (A corresponding momentum eigenstate involves the 
superposition of many local configurations for which that site has different positions.) In turn, if annihilation of a 
local c fermion at its effective lattice involves annihilation of a local si fermion, one has the following virtual processes 
in the spin effective lattice: First, a local si fermion is annihilated, what involves a rearrangement process, which 
leads to the occupancy of its two sites by two independent -1-1/2 spinons in the intermediate virtual state; Second, 
the site of the spin effective lattice occupied by one of these two independent spinons is removed along with it. 

Creation (and annihilation) of both one local c fermion and one local si fermion corresponds to creation (and 
annihilation) of a spin-down electron. In turn, the corresponding process of creation (and annihilation) of a spin-up 
electron involves addition (and removal) of one local c fermion to (and from) its effective lattice and addition (and 
removal) of one site occupied by an independent -1-1/2 spinon to (and from) the spin effective lattice. Alike creation 
and annihilation of local c fermions leads to addition and removal (and removal and addition) of sites in the spin (and 
77-spin) effective lattice, respectively, it is confirmed below that creation of one local av' fermion gives rise to addition 
of 21/' sites to the ai^ effective lattices of aiy fermion branches such that v < v' . 

It is confirmed below in Section V that for the states that span the one- and two-electron subspace only c fermions 
and si fermions play an active role. Concerning excitations belonging to such a subspace that conserve the number 
A^c = 25c of c fermions and thus the eigenvalues Sc of the generator ([7]) of the global C/(l) symmetry and the 
number Nsi of si fermions, the c fermions can move through their effective lattice independently of the motion of the 
corresponding si fermions. In that case when a given c fermion moves around in the lattice it is not attached always 
to the same spinon and hence to the same two-spinon si fermion. The investigations of Ref<^ reveal that there are 
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correlations between a spin- neutral two-spinon si fermion and a well defined set of c fermion pairs whose centre of 
mass is located at the si fermion position. 

C. Ranges of the av fermion and c fermion energies, transformation laws of such objects, and the 

ground-state occupancies 

The quantum-object occupancy configurations of the ground state are found below. Self-consistency of our descrip- 
tion then confirms that for m > and a; > the elementary energies £3,-1/2 = H and -1/2 — 2/i of Eq. \A7\ of 
Appendix A correspond to creation onto the ground state of an independent —1/2 spinon and an independent —1/2 
?7-spinon, respectively. We recall that /is is the Bohr magneton and H the magnitude of a magnetic field aligned par- 
allel to the plane (2D) or chain (ID). The energy €s^-i/2 = 2/xb H (and e,;,_i/2 = 2/i) refers to an elementary spin-flip 
(and 77-spin-fiip) process, which transforms an independent -1-1/2 spinon (and -1-1/2 ?7-spinon) into an independent 

— 1/2 spinon (and —1/2 77-spinon). Such elementary energies control the range of several physically important energy 
scales. Since within our LWS representation an independent -1-1/2 spinon (and -1-1/2 77-spinon) has vanishing energy 
and an independent —1/2 spinon (and —1/2 vy-spinon) has an energy given by £3,-1/2 = 2/iBi? (and £,,,-1/2 — 2/i), 
the energy of a pair of independent spinous (and 77-spinons) with opposite projections is 2/is H (and 2/^). Indeed, due 
to the invariance of such objects under the electron - rotated-electron unitary transformation, they are not energy 
entangled and the total energy is the sum of their individual energies. 

The magnetization curve is such that the spin density m vanishes at if = and is finite and positive for H > 0. 
(We recall that within the LWS representation the convention that > for x > and H > for m > is used.) 
Consistently with the properties of the x = and rn = absolute ground state discussed below and in Ref.^, for x ^ 
and m = the chemical potential /i belongs to the range /i £ (— Z^*^, Z^*^) where the energy scale = \im.x^o /i equals 
one half the Mott-Hubbard gap and is such that ^° ^ for U/At and ^° > for U/At > 0. For < j: < 1 and 
m — the chemical potential is an increasing function of the hole concentration x such that, 

/ < ^^{x) < n^; 0<x<l,m = 0, (32) 

where /x^ = lima;_).i /j,. fj} reads, 

^'^ = U/2 + 2Dt- L> = 1,2. (33) 
In turn, has the following limiting behaviors, 

6-2^^(17)'^°, C//4t<l; / w [[//2 - 21)^] , L//4t>l, ^ = 1,2. (34) 

Hence indeed — >■ as [//4t — whereas /i" cx) for U/At 00 for both the model on the ID and square lattices. 

For D = 1 the expressions provided in Eq. are obtained by use of the exact Bethe-ansatz solutionis. 

Expression (|33| is exact both for ID and the square lattice. It can be derived explicitly for both lattices, since it 
refers to the non-interacting limit of vanishing electronic density. In turn, the expressions given in Eq. (j34]) for 
the model on the square lattice are obtained by the use of a result discussed in Subsection VI-A: the energy below 
which the long-range antiferromagnetic order survives for x = 0, m = 0, and zero temperature is where 2^q 
is the Mott-Hubbard gap. Indeed, according to the analysis of that subsection the Mott-Hubbard gap 2/i'^, which 
refers to the charge degrees of freedom, affects the spin degrees of freedom as well. Consistently, the D = 2 limiting 

behaviors w 32te'^''V^ and fi° « U/2 of Eq. ^ for U/U < 1 and U/U > 1, respectively, are those of the 
zero-temperature spin gap of Eq. (13) of Ref.^^. 

In the following we confirm that ground states have no rjv fermions and no sv fermions with v > 2 spinon pairs. 
Hence the corresponding energies £,,jy and £3^, respectively, considered below refer to creation onto the ground state 
of one of such objects carrying a given momentum. We start by providing a set of useful properties. We emphasize 
that some of these properties are not valid for descriptions generated by rotated-electron operators associated with 
the general unitary operators V considered in Refi^. Indeed, the following properties refer to the operator description 
associated with the rotated-electron operators ct ^ — V'< cl, of Eq. dS]) as defined in this paper, whose electron 

- rotated-electron unitary operator V = V{U/4£) is such that the set of 4^° states of the form \^u/At) = ^^^l^'oo) 
are energy eigenstates for U/At > 0. Here l^'oo) corresponds to exactly one of the many set of f//4t — >• 00 energy 
eigenstates, chosen according to the well-defined criterion discussed above in Subection II-A. For the model on the 
square lattice our operator description refer to a related complete set of 5,,, S'^, Sg, SI, Sc, and momentum eigenstates 
|'J*c//4t) — ^^^I'l^oo) generated from the corresponding U/it 00 states |$oo) by the same electron - rotated-electron 
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unitary operator V — V{U/4:t) as the energy eigenstates |5'c//4t) = V^^I'I'oo)- Some of the states \^u/it) are identical 
to the energy eigenstates |^c//4t), which are also Srj, Sf^, Sg, S'f , Sc, and momentum eigenstates. (For the ID model 
both sets of states are identical for the whole Hilbert space.) 

Let us profit from the interplay of the transformation laws of the ai^ fermions under the electron - rotated-electron 
unitary transformation with the model global 5*0(3) x SO{3) x U{1) symmetry to reach valuable information about 
the range of the energy for addition onto the ground state of one av fermion. For the model on the square lattice 
there is in general a small uncertainty Se in the energy eav of a fermion of a given momentum smaller than or 
equal to the corresponding energy bandwidth and such that ± Se belongs to that energy range. For that model 
and ai^ fermions that are invariant under the electron - rotated-electron unitary transformation one has Se — 0. That 
holds for instance for the av fermions whose occupancy configurations generate above states \^u/4t) that are identical 
to energy eigenstates \'^u/4t)- For both the model on the ID and square lattices the latter range corresponds to the 
energy bandwidth of the energy dispersion eav For the description associated with rotated-electrons, corresponding 
c fermions, 77-spinons, and spinous whose operators are given in Eqs. l|8|)-(|17p. and composite ry-spinon and spinon 
fermions the following properties hold. 

1. Energy range of a rjv fermion 

A rjiy fermion is a 77-spin-ncutral anti-bounding configuration of v —1/2 77-spinons and v +1/2 77-spinons. Symmetry 
implies that for U/At > there is no energy overlap between the ranges of bandwidth Wjju of the energy erju 
corresponding to different branches = 1, C,; associated with addition onto the ground state of one rjv fermion of 
a given momentum. Since fermions belonging to neighboring branches r/v and 777^ + 1 differ in the number of single 
77-spinon pairs by one, the requirement for the above lack of energy overlap is that Wr/i, < 2/i where 2/i equals the 
energy of a pair of independent ?7-spinons of opposite 77-spin projections. Such properties imply that the energy e^ji, 
obeys the following inequality, 

2j//i < e^i, < 2(i^ + 7^^)/i ; < i^^ < 1 , (35) 

where 2vfi is the energy for creation of v —1/2 independent 77-spinons and v independent -1-1/2 77-spinons onto the 
ground state. Since the latter objects are invariant under the electron - rotated-electron unitary transformation 
associated with the operator V, they are non interacting and thus their energies are additive. For m — the number 
ir/i, decreases continuously for increasing values of U /At, having the limiting behaviors 7^1/ — )■ 1 for U / At and 
iriu — > for U /At ^ 00 so that Wrfv — as J7/4t — 00. The latter behavior is associated with the full degeneracy of 
the 7;-spin configurations reached for U /At — >■ c» when the spectrum of the 27/-77-spinon composite 777^ fermion becomes 
dispersionless. 

2. Energy range of a sv fermion 

A sv fermion of a given momentum is a spin-neutral bounding configuration of —1/2 spinons and v +1/2 spinous. 
Again, symmetry implies that for U/At > there is no energy overlap between the ranges of bandwidth Wsi, of the 
one-sv fermion energy esi, associated with different branches v = 1, ...,Cs. For sv branches with a number of spinon 
pairs V > 2 the bandwidth Wgi, of such an energy range is for the present bounding configurations and for the same 
reasoning as for the r]v fermion such that Wsu ^ H where 2fiB H equals the energy of a pair of independent 
spinons of opposite spin projections. The energy esu for addition onto the ground state of one sv fermion with v > 2 
spinon pairs obeys the inequality, 

'2{v - is^)^J.B H < es^ <2vfiB H ; v>l, < i^^ < 1 , (36) 

where 2vfiB H is the energy for creation of v independent —1/2 spinons and v independent +1/2 spinons onto the 
ground state. Since the latter objects are invariant under the unitary transformation associated with the operator V, 
they are not energy entangled and their energies are additive. For 777 = the number i^i, decreases continuously for 
increasing values of U/At, having the limiting behaviors igi, — ?> 1 for U/At — > and ig^ — > for U/At — > 00 so that 
WsL> — as J7/4t — ^ 00. Such a behavior is associated with the full degeneracy of the spin configurations reached for 
U/At — ^ 00 when the spectrum of the 2i'-spinon composite sv fermion becomes dispersionless. 

In turn, for a tti = and x > ground state all sites of the si effective lattice are occupied so that the corresponding 
si momentum band is full and the energy — e^i for removal from that state of one si fermion carrying momentum 
obeys the inequality, 



0< -eai <max{VF,i,|A|}, 



(37) 
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where the energy scale |A| is introduced below in Subsection VI-B and further studied in Ref.^. For small hole 
concentrations < x ^ 1 it vanishes both in the limits U /At — >■ and U /At — >■ oo and goes through a maximum value 
at U/At — uq 1.302 and at constant U/At decreases for increasing x and vanishes for x > x^, where the critical hole 
concentration x» is studied in Ref.— for approximately U/At > uq. As discussed in that subsection, due to a sharp 
quantum phase transition it has a singular behavior at x = 0, having for U/At > different magnitudes at x ~ 
and a; — > 0, respectively. In turn, for m = the energy bandwidth Wsi refers to the auxiliary si fermion dispersion 
defined in Ref.^ and has it maximum magnitude at U/At = 0. For U /At > it decreases monotonously for increasing 
values of U/At, vanishing for U/At — > oo. That both Wsi — ^ and |A| — > for J7/4t — > oo is associated with the full 
degeneracy of the spin configurations reached in that limit for which the spectrum of the two-spinon composite si 
fermions becomes dispersionless. For instance, we could access the explicit limiting behaviors of the m = energy 
bandwidth W^^ = lim^^^o Wsi = Wsi\x=o, 

W^j^=2Dt, U/At^O; W^^^ ^ 2071 — , [//4i > 19^"^ L> = 1, 2 , (38) 

where that for the D = 2 square lattice the laige-U/At expression is valid for U/At ^ 19 is justified below in Section 
VI. 

In the U /At — > oo limit the si fermion occupancy configurations that generate the spin degrees of freedom of spin- 
density m = ground states considered below become for the ID model those of the spins of the spin-charge factorized 
wave function introduced both by Woynarovich*^ and Ogata and Shiba^'^. In turn, for the model on the square lattice 
such configurations become in that limit and within a mean- field approximation for the fictitious magnetic field Bsi 
of Eq. (pSj) those of a full lowest Landau level with Ngi = N^^_^ = N/2 one-sl-fermion degenerate states of the 2D 
quantum Hall effect Here iV^^^ is the number of both sites of the square si effective lattice and si band discrete 
momentum values. For finite U/At values and a; > the degeneracy of the N^^_^ = N/2 one-sl fermion states of the 
square-lattice quantum liquid is removed by the emergence of a finite-energy-bandwidth si fermion dispersion yet 
the number of si band discrete momentum values remains being given by N^^_^ — Bgi i^/^o and the si effective 
lattice spacing by agi — Isi/\/2tt where Isi is the fictitious- magnetic- field length and we recall that in our units the 
fictitious-magnetic-field flux quantum reads $o = 1~- 

3. The energy range of the c fermtons 

The energy ec for addition onto the ground state of one c fermion of a given momentum and the energy — Ec for 
removal from that state such a c fermion obey the inequalities, 

< Cc < M^i" = - WP] ; < -e^ < W^^ , = 1, 2 , (39) 

respectively. Here = [ADt - W/!] € (0, ADt) increases monotonously for increasing values of hole concentration 
X S (0, 1). The energy bandwidth depends little on U/At and for U/At > has the following limiting behaviors, 

WP =^ADt, x^O; WP = 0; x^l. (40) 

4- Transformation laws of au fermions and c fermions under the electron - rotated- electron unitary transformation 

An useful property is that rjv fermions and sv fermions with v > 2 spinon pairs that remain invariant under the 
electron - rotated-electron unitary transformation have energy given by, 

erju^2v^, = 1, ; ^2v^ibH , v ^2,...,Cs- (41) 

Those are non-interacting objects such that their energy is additive in the individual energies of the corresponding 
2y r?-spinons and spinous, respectively. Therefore, for U/At > they refer to the same occupancy configurations 
in terms of both rotated electrons and electrons. However, note that in contrast to the independent 77-spinons or 
independent spinous here the objects that are invariant under the electron - rotated-electron unitary transformation 
are the composite 2z/-77-spinon r]v fermions or 2z^-spinon sv fermions and not the corresponding individual 2v ry-spinons 
or 2;/-spinons, respectively. In turn, "qu fermions and sv fermions with v > 2 spinon pairs whose energies obey the 
inequalities e^i^ > 2i^/^ and esi^ < 2u^b H are not invariant under the electron - rotated-electron transformation. 
Furthermore, for finite U/At values c fermions and si fermions are not invariant under that transformation. 

Both for the model on the ID and square lattices the initial ground state and the excited states generated by 
creation of one rju fermion or sv fermion whose energy is given by (j4ip are energy eingenstates. 
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5. Ground state occupancies 

Fulfillment of the requirement of self-consistency concerning the set of properties given here, reveals that both for 
the model on the ID and square lattices the subspace spanned by the x > and m > LWS ground states and their 
excited energy eigenstates of energy uj < min {2/x, 2/iB i?} there are no —1/2 r;-spinons, rjv fermions, independent 
— 1/2 spinous, and sv' fermions with v' >2 spinon pairs so that A^,,,y = and Nsu' = for v' > 2. Hence the number 
of c fermions is Nc = 2Sc ^ N = (1 - x) , independent +1/2 77-spinons ^,,,+1/2 = 25,, = [iVf - N] = x , 
independent +1/2 spinous = 25^ — [N^ — N^] = mN^, and si fermions Nsi — Ni- 

From analysis and comparison of the occupancies of the spin LWS ground states {m > 0) and spin highest-weight 
state (HWS) ground states (m < 0), one finds that a m = ground state for which N is even and a; > has no —1/2 
77-spinons, riv fermions, independent ±1/2 spinous, and sv' fermions with v' > 1 spinon pairs so that N^ii, = and 
Ns^' — for i/' > 1. Hence the number of c fermions is Nc — 2Sc = N = {1 — x) , independent +1/2 Ty-spinons 
= 25,, = [iVf - N] =xN^, and si fermions iV^i = N/2 = (l - x) N^/2. 

D. The number of sites of the au effective lattices and discrete momentum values of the av bands 

A local rji' (and sv) fermion refers to a superposition of well-defined ?7-spinon (and spinon) occupancy configurations 
involving 2^ sites of the ry-spin (and spin) effective lattice. The number of sites of the si effective lattice plays an 
important role in the study of the model in the one- and two-electron subspace considered below in Subsection V-A. 
In order to evaluate its expression one needs to solve the same problem for the remaining av effective lattices as well. 
Indeed, the values of the set of numbers {-/V^P } are dependent of each other. Fortunately, in the limit iV^ ^ 1 
one does not need detailed information about the occupancy configurations of the 77-spin (and spin) effective lattice 
whose superposition defines the internal structure of a rjv (and si') fermion to achieve such a goal. The only needed 
information is that the 2^ sites of the ry-spin (or spin) effective lattice involved in such configurations are centered at a 
point of real-space coordinate fj, which defines the real-space coordinate of the local av fermion. Here j = 1, N^^^^ 
of a given subspace where -/V^_^ is the number of sites of the av effective lattice. Our goal is the derivation of the 
expression for that number in terms of the set of numbers Nc and {Na'u'} where a' — i], s and v' = 1, 2, 3, Ca'u' ■ 
The real-space coordinate Vj plays the role of "centre of mass" of the local aiy fermion and has A'^ well-defined 
values associated with the sites of the av effective lattice. 

For A^^ ^ 1 (i) the internal structure of a local av fermion and (ii) its real-space position rj are separated problems. 
The present analysis refers to the problem (ii) only. The problem (i) is addressed for the si bond particle associated 
with the local si fermion in Ref.'^®. Concerning the internal structure of a local fermion, the only issue that 
matters for the present analysis is that the 2i/ sites of the 77-spin (and spin) effective lattice occupied by a given local 
771/ (and sv) fermion correspond to 2i/ sites of the original lattice that are not simultaneously occupied by any other 
such fermions. 

The number of sites N^^^ of the ai' effective lattice is an integer number. We emphasize, however, that for the 
model on the square lattice for which D = 2 the related number Na^^ is not in general integer so that the av effective 
lattice is not a perfect square. However, for Na^^/Na finite and N^ ^ 1 it is nearly a square lattice so that in that 
limit we use the notation A'^^^^ for its number of sites. 

The aiy effective lattice and its A^^^^^ sites are well-defined concepts in a subspace for which the values of the set 
of numbers Nc = 2Sc and {Na,^'} where v' = 1,2,3, ...,Ca remain constant. As justified below, this is equivalent 
to the 7;-spin 5,, {a — 77) or spin 5s (a = s) and values of the set of numbers {Na,y'} where v' = i>,v 1, ■■■,Ga 
remaining constant. For such subspaces the motion of the av fermion through its effective lattice corresponds to a set 
of elementary steps where it hops from a given site of well-defined real-space coordinate to another. For the model 
on the square lattice such elementary steps involve horizontal and vertical virtual stepa^^. 

For a local av fermion the 2v Nav sites of the 7/-spin [a = rj) or spin {a = s) effective lattice occupied by the Nau 
local fermions belonging to the same av branch play the role of the Nai, occupied sites of the av effective lattice. In 
turn, the 25,, (and 25s) sites of the 77-spin (and spin) effective lattice occupied by independent +1/2 77-spinons (and 
spinons) and 2{v' ~ v) sites out of the 2v' sites of that lattice occupied by each local 771^' (and sv') fermion such that 
v' > V play the role of unoccupied sites of the rjv (and sv) effective lattice. 

The number of sites of the av effective lattice is then given by. 



(42) 
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where the number of unoccupied sites reads, 

Co Co 

i^' — iy+l y' —1 

The equivalence of the two expressions given here confirms that the numbers Nc — 2Sc and {Na^r} where v' = 
1,2,3, ...,Cq remaining constant is equivalent to the 7y-spin Srj {a = rf) or spin Ss (a = s) and values of the set of 
numbers {Na^:} where v' = i/, + 1, Ca remaining constant as well. In both cases that implies that N^^^ remains 
constant. 

In the following and in Appendix C it is confirmed that the expressions given in Eqs. ()42p and ()43p . which equal 
the av band corresponding numbers of discrete momentum values, are compatible with the number of representations 
of the group 50(3) x 5*0(3) x U{1) in each subspace with constant values of 5c, 5^, and Ss- From the use of these 
equations one finds that the number of unoccupied sites of the al effective lattices reads, 

Co 

N^i^[NZ~2B^]- B^^Y.N^,- a = v,s. (44) 

This number equals that of al fermion holes in the al band. 

Straightforward manipulations of the above equations lead to the following general expressions for 5^, Ss, and 5c, 

1 A[D 1 1 

5a = [2<i - 0„ + Sc.] , a^rj,s; S, ^ ~ -N^;, - B,, = -N^, + Bs] . (45) 
It follows from the equality of the two 5c expressions given Eq. (|45)) that. 
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Equation (PT|) . the first expression of Eq. P5|) . and Eq. are equivalent to the following sum-rules for the 

numbers of av fermions, 

Co ^ Co ^ 

I/— 1 i^— 1 

respectively, consistently with the expressions Cs = [Sc — Ss] and O,, = [N^ — 5c — 5„] given in Eq. ([2T|) . 

Each subspace with constant values of 5c and hence also with constant values of = [N^ — 25c] and 25c 
that the vacuum of the theory given in Eq. ([23| refers to can be divided into smaller subspaces with constant values 
of 5c, 5^, and Ss and hence also with constant values of O^ and Cs- Furthermore, the latter subspaces can be further 
divided into even smaller subspaces with constant values for the set of numbers {A^i/} and \N sv\, which must obey 
the sum-rules of Eq. (PT|) . 

Let us confirm that expressions (|^^ and (|33]) for the number of discrete momentum values of the av band and 
av fermion holes in that band, repectively, are compatible with the number of state representations of the ?7(1), 
7?-spin 5C/(2) (a = r?), and spin 5C/(2) (a = s) symmetries of the model global 50(3) x 50(3) x U{1) = [5C/(2) x 
SU{2) X t/(l)]/Z| symmetry in the subspaces with constant values of 5c, 5^, and Ss- Those are subspaces of the 
larger subspace that the 5c > vacuum of Eq. (1^51) refers to. Alike in Reff^,, let us divide the Hilbert space of 
the model ([T]) in a set of subspaces spanned by the states with constant values of 5c, 5,,, and 5s and hence also of 
A^c = 25c, Mr, = A^i^ = [N^ - 25c], and Ms = A^i^ = 25c. We recall that for the subspace with constant values of 
5c, Sri, and Ss under consideration Mjj = N^^ (and AIs ~ A'i^) is both the total number of ?7-spinons (and spinous) 
and the number of sites of the ?7-spin (and spin) effective lattice. 

According to the studies of Ref.^ the dimension of each such a subspace is. 



dr- n AA(5c,Af„), (48) 



where dr and Af{Sa, M^) are given by. 
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and are the number of U{1) symmetry state representations and that of 77-spin SU{2) {a = rf) or spin SU{2) {a = s) 
symmetry state representations, respectively. 

The dimension dr given in Eq. ()49|) is here straightforwardly recovered as dr — ) and equals the number of 
occupancy configurations of the Nc = 2Sc c fermions in their c band with discrete momentum values. On the 
other hand, the values of the numbers N^^^^ of discrete momentum values of the aiy band must obey exactly the 
following equality for all subspaces. 

Here (^'"') is the number of occupancy configurations of the Nai, av fermions in their av band with N^^^ discrete 
momentum values and the {Nav} summation runs over all sets of Na^ numbers for v — 1^2, Ca that owing to the 
conservation of Ca exactly obey the subspace sum-rule. 



Ca 

2Ca = 51 = [^^" - 2^«] , a = 77, s ■ (51) 



The general expression of the number N^^^ of av band discrete momentum values that obeys Eq. ([50|) for all 
subspaces corresponds indeed to that given in Eq. (1421) . Consistently, the occupancies of the independent ?7-spinons 
and independent spinous give rise to the usual factors (25*^ + 1) and (25*3 + 1), respectively, appearing in the expressions 
provided in Eqs. (|49l) and (|50p for the dimensions Af{Srj, M,f) and JV{Ss, Mg). In Appendix C it is shown that the 
subspace-dimension summation, 

A^-* - E E E ffc ) n ^^^^0 ^^") - 4""° , (52) 



Sa=0 s„=o s,=o 



leads indeed to the dimension 4^» of the full Hilbert space. It follows that the present description in terms of 
momentum eigenstates is complete. 

In the 1 limit considered here and alike for the ?7-spin and spin effective lattices, provided that N^^^/N^ 

is finite the related ai/ effective lattices can for ID and 2D be represented by ID and square lattices, respectively, of 
lattice constant, 

■a = -r7 — flc,; iVa„„>l, (53) 



Na^. Na^. Na^. 

where f = 1, Ca and a = rj, s. In turn, the corresponding av bands whose number of discrete momentum values 
is also given by N^^^^ are well defined even when N^^^^ is given by a finite small number, = 1, 2, 3, ... 

For the model on the square lattice and states with numbers Nai, = -^a„„ fo^' '^'^ 7^ *1 ^^'^ ^si ~ -^a^i fo'" '^'^ = -^1 
where [N^^_^ — Nsi] vanishes or is of order for -> 00, the lattice constant aau is directly related to the 

fictitious magnetic-field length lai, associated with the field of Eq. ([25)1 . Indeed, in that case one has that {rip. au) ~ 1 

and such a fictitious magnetic field reads Ba^ifj) ~ $0 X^j'^^j ^i^j' ~ ^j) ^^a- It acting on one av fermion differs from 
zero only at the positions of other av fermions. In the mean field approximation one replaces it by the average field 
created by all av fermions at position fj. This gives, Bau{rj) « $0 nau{rj) ~ $0 Wa^^/L"^] e^a = [^J^o/Oai/] ■ 
One then finds that the number N'^^^ of the av band discrete momentum values equals [i3Q.,yL^]/$o and the av 
effective lattice spacing aav is expressed in terms to the fictitious magnetic-field length lav as a^^, — 2^1"^^. This is 
consistent with for such states each av fermion having a flux tube of one flux quantum on average attached to it. For 
the states under consideration the av fermion problem is related to the Chern-Simons theory^, the number of fiux 
quanta being one being consistent with the av fermion and av bond-particle wave functions obeying Fermi and Bose 
statistics, respectively. Hence the composite av fermion consists of 2v 77-spinons (a — rf) or spinous (a — s) plus an 
infinitely thin flux tube attached to it. 

The states that span the one- and two-electron subspace belong to such a class of states so that for the model on the 
square lattice in that subspace the composite si fermion consists of two spinous in a spin-singlet configuration plus 
an infinitely thin fiux tube attached to it^. Thus, each si fermion appears to carry a fictitious magnetic solenoid with 
it as it moves around in the si effective lattice. Finally, since for such a quantum liquid the si fermions play a major 
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role we find from straightforward manipulations of Eqs. (|42p and (|43p for ai' — si that the number N^^_^ = [Nsi + N^j] 
of si effective lattice sites and thus of si band discrete momentum values is given by, 

iVi^, = [5c + 5. + E - 2)^-] ■ (54) 

i/=3 

Hence for a subspace where Ngi, = for > 3 one has that N/;^^ ~ [Sc + Sg] is a good quantum number. 

E. The momentum eigenstates of our description and the good and quasi-good quantum numbers 

1. The momentum eigenstates of our description 

The low-energy physics of a 3D perturbative and isotropic many-electron system for which each energy eigenstate 
|\E'(0)) of the corresponding non-interacting system evolves upon adiabatically switching on the interaction U on an 
energy eigenstate |^([/)) is in general successfully described by Fermi liquid theory^^. (In contrast to the Hubbard 
model, here the interaction U is not necessarily onsite.) Since the two sets of orthogonal and normalized states {|^(0))} 
and {|\I'(?7))} are complete and refer to the same Hilbert space, there is a well-defined unitary transformation such 
that for each U = Q energy eigenstate |5'(0)) there is exactly one C/ > energy eigenstate Y^{U)) = V^|^'(0)) where 
"pt = Vt(C/) is the corresponding unitary operator. For the exact ground state such a statement is equivalent to the 
Gell-Mann and Low theorem**"*. By perturbative system we mean above that (\I'(0)|\1/(C/)) 7^ for all such pairs of 
states of each of the two complete sets. 

If the t7 = energy eigenstates can be written as Slatter determinants of the form, 

i*(o))=nn4,jo)' (55) 
k 

where ct creates an electron of momentum k and spin projection a and |0) is the electronic vacuum, then the 
corresponding states |*([/)) = Vt|^(0)) read, 

i*(^))=nn^Mvi«>=nn[^H.^]^'io)=nn4,jo)= c-^'4,.^= v^io)=io). (56) 

k k " k 

Here ct = ct is the creation operator of a quasiparticle of momentum k and spin projection a. 

Note that since for U/At — > the unitary operator becomes the unit operator, a state |vE'(0)) of Eq. (|55l) equals 
the corresponding state |\['([/)) of Eq. ([55]) for U/At 0. Let us consider the one-electron Green function, 

G{k,a,t-t') = -t{^Gs{U)\Tc^Jt)cljt')\^Gs{U)), (57) 

where for (7 > (and [/ = 0) the A^+1 states ct (t)|*Gs(J7)) and ct {t')\^Gs{U)} and A^-l states ct: (t)\^Gs{U)} 
and ^(t')|\E'Gs(C^)) are not (and are) energy eigenstates. In turn, for [/ > the 'W -I- 1 quasiparticle states" 
dijt)\^Gs{U)) and £ljt')\^GsiU)) and "iV - 1 quasiparticle states" £f:^^{t)\'^GsiU)) and cg ,^(f )|*gs(C/)) are 
energy eigenstates, yet they cannot be proved in real experiments because only the electrons are physical particles. 
Indeed, quasiparticles exist only inside the many-electron system. Therefore, quasiparticles can only be added to or 
removed from the system upon addition to and removal from it of electrons. That for U > the above one-electron 
states are not energy eigenstates implies that the renormalization factor Z{k) is smaller than one and for momenta 
k in the vicinity of the isotropic Fermi surface and small excitation energy uj such states have an inverse lifetime 
proportional to cj^ so that the "quasiparticle lifetime" is finite for w > and becomes infinite as w — > 0^^. 

We emphasize that in spite of such an one-electron lifetime being finite, the quasiparticle occupancy configurations 
()56p generate for U > exact energy eigenstates, whose continuous changes with time occur according to the many- 
electron Schrodinger equation. Concerning one- and two-electron excitations, the concept of a quasiparticle is well- 
defined for low-energy eigenstates whose quasiparticle occupancy configurations differ from those of the initial ground 
state only near the Fermi surface. In turn, in the absence of one- and two-electron quantum measurement processes 
the quasiparticle occupancy configurations as defined here generate all C/ > energy eigenstates. Hence solution 
of the quantum problem is equivalent to derivation of the unitary operator — V^(J7) for U > 0. Unfortunately, 
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this is in general a very involved and unsolved problem. Landau's Fermi liquid theory is a very useful scheme that 
allows solution of it for the subspace spanned by excited energy eigenstatcs generated from the initial U > ground 
state by a finite number of quasiparticle processes near the Fermi surface. 

In turn, for the Hubbard model on the square or ID lattices our choice of 4^" energy eigenstates {|^oo)} for 
U/it — >■ oo reported is, as in Subsection II-A, such that the states |vE'(7/4t) = t^^l^'cjo) are energy eigenstates for finite 
values U/At > 0. Following the result of Appendix A that the general problem of that model in the whole Hilbert 
space can be fully described by the same model in the LWS subspace, for simplicity to start with our discussion refers 
to that subspace. That just means that for states with independent 77-spinon and/or independent spinon occupancy 
the 77-spin and spin projection, respectively, of such objects is given by +1/2. Within our operator description for 
the model on the square lattice each of the J7/4i — 7> 00 LWS energy eigenstates of the corresponding set {|^'lwS;oo)} 
can be expressed as a suitable superposition of Sjj, S^, Ss, S^, Sc, and momentum eigenstates {|<&lws:oo)} with the 
same values for such physical quantities and also for the numbers Cq = Yl^Zi ^ ^au of Eqs. (PU)) . (pijl . ([5T|). and (|T7)) 
where a = 77, s. The latter states can be generated by occupancy configurations of the microscopic momenta q carried 
by the c and av fermions, which as justified below are of the general form, 

ii>w5;oo) = [nnn4,-iO"'^<)Hn4ciG^=;o)] • (s^) 

a v q' q 

Here J^t, and are creation operators of a C//4t — 00 au fermion of momentum q' and c fermion of momentum (f, 
respectively, and the ry-spin SU{2) vacuum |0^; iV^^ ) associated with N^^ independent +1/2 ?7-spinons, the spin SU{2) 
vacuum |0s; A^^) with N^^ independent +1/2 spinous, and the c U{1) vacuum \GSc',Q) such that Y\q'^\ c\^^c';Q) = 
\GSc\'2,Sc) has Nc — 2Sc c fermions are those of the vacuum given in Eq. (1^51) . We recall that for U/4t — > 00 Eqs. 
([5|)- P?)l are equivalent to Eqs. (l)-(3) of RefJ^ with the c fermion creation operator fl ^ replaced by the quasicharge 

annihilation operator Cr- Therefore, the operator J-t is given by, 

q.c 

(59) 




where within the notation of Ref.— , r = fj are the real-space coordinates of the quasicharge particles of that reference. 
The holes of such quasicharge particles are for J7/4i — >■ 00 the spinless fermions that describe the charge degrees of 
freedom of the electrons of the singly occupied sites. Moreover, the spinous and ?7-spinons of the 2z^-spinon operators 
tI, and 2:/- 77-spinon operators Tt, respectively, are associated with the local spin and pseudospin operators, 
respectively, defined in that reference. As discussed in Appendix B, for the ID model all states ([55)) are energy 
eigenstatesi, so that lwS;cx,} = \^lwS:oo)- For the model on the square lattice that holds only for some of these 
states. For instance, it holds for the energy eigenstates that span the one- and two-electron subspace as defined below 
so that \'^lwS;od) = \^LWS;oa) for that subspace. 

Symmetry implies that for the Hubbard model on a square lattice at U/'it > the N^^^ av band discrete momentum 
values are J7/4t independent for states belonging to the same V tower. Importantly, the vacua of the C//4i — > 00 
states (|58p are invariant under the electron - rotated-electron unitary transformation so that for finite values of U /At 
the corresponding 5*,,, 5*,^ Ss, S^, Sc, and momentum eigenstates read, 

\'^Lws;u/,t) = [IT nn 4.0.10-^'^?)] [1141^^-0)] ; = ^^4.-^; 4 = (60) 

a u q' q 

Alike in the U /At — > oo limit, for the one- and two-electron subspace the energy eigenstates \^ lwS;U /it) of the model 
on the square lattice are such that \^ LWS;U/At) = \^ LWS;U/At) , whereas for the ID model the equality LWS;U/At) = 
\^LWS\U/it) holds for all U/At > energy eigenstatesii^. In Appendix B it is confirmed that the discrete momentum 
values of the c and fermion operators appearing on the right-hand side of Eq. (pO)) are indeed the good quantum 
numbers of the exact solution whose occupancy configurations generate the energy eigenstates. Furthermore, in 
that Appendix the relation of the c and fermion operators whose suitable products refer to the operators whose 
application onto the theory vacua generates the states ([50)1 to the creation and annihilation fields of the charge ABCDF 
algebraic and more traditional spin ABCD Faddeev-Zamolodchikov algebra^ of the algebraic formulation of the ID 
exact solution of Ref.^® is discussed and the consistency between the two corresponding operational representations 
confirmed. 

We recall that addition of chemical-potential and magnetic-field operator terms to the Hamiltonian (|T]) lowers its 
symmetry. However, such operator terms commute with that Hamiltonian and the momentum operator so that its 
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energy and momentum eigenstates correspond to state representations of the 50(3) x 5*0(3) x [/(I) group for all 
values of the densities n = (1 — x) and m. The same holds for the 4^= momentum eigenstates considered here, whose 
c fermion and av fermion occupancy configurations and independent T^-spinon and spinon occupancies are of the form, 

\^u/,t)= n ^^"^^'''^ \'^>Lws-M/4t) = [H ^^"^"^'"'^ nn4'.ajO";^ag)][n4,cig^^;o)] ■ (6i) 

a—rj.s vn OL ^ u q' q 

Such states correspond indeed to the state representations of the 50(3) x 50(3) x U{1) group. In this expression 
\^LWS;U/it) are the states of Eq. (pO]) and the normalization constant Ca is given in Eq. (jA2l) of Appendix A. 

The set of 4^° momentum eigenstates of form (|6T|) is complete and corresponds exactly to the c fermion and av 
momentum occupancy configurations associated with the subspace dimensions of Eq. (|48p and the corresponding 
subspace dimensions, 

dr ■ n (^-^a + 1) • AA(5a, M„) , (62) 
a—rj.s 

where dr and J\f{Sa, Ma) are given in Eq. (1491) and here we have accounted for the occupancies of the independent 
ry-spinons and independent spinous associated with the usual factors (25*,, + 1) and (25*5 + 1), respectively. Use of the 
summations of the subspace dimensions of Eq. (|62l) performed in Appendix C confirms that the correct Hilbert space 



dimension 4^" is indeed obtained. (For the ID model one has that such 4^» states {|'&(7/4t)} are both momentum 
and energy eigenstates so that |^[//4t) — \^u/it)-) 

For the the operators ft, there is an independent problem for the model in each subspace with constant values of 
5c, 5s, and number of sites of the ai/ effective lattice and discrete momentum values of the az^ band 7V^^ = [Na^+N^^] 
given in Eq. P^ . Hence the operators fl and fq-j^av act onto such subspaces, where they are the building blocks of 
the generators that change the av fermion occupancy configurations. Hence and as mentioned above, such operators 
act onto subspaces spanned by mutually neutral states and that assures that for the model on the square lattice the 
components qxi and qx2 of the microscopic momenta q — [qxi, qx2\ refer to commuting si translation generators cjsi xi 
and qsix2~- In turn, for the operators ft ^ there are in ID two types of quantum problems depending on the even 
or odd character of the number [B,^ + Bg] — J^a v ^av Indeed, in spite of the c effective lattice being identical to the 
original lattice for the whole Hilbert space, as discussed in Appendix B the c fermions feel the Jordan- Wigner phases 
of the av fermions created or annihilated under subspace transitions that change the number [B,, + Ss]. Hence for 
ID the c band discrete momenta have different values for [i?^ + B^ odd and even, respectively. Also for the model 
on the square lattice the c fermions feel the Jordan- Wigner phases of the av fermions created or annihilated under 
subspace transitions that change the number [B,, + B^. 

The generators that transform into each other the states with constant values of 5c, S's, and A''^^ = \Nav + 
which span such a subspace, can be expressed in terms of creation and annihilation av fermion operators. These 
states can be formally written as given in Eqs. (1601) and (j6ip provided that the av momentum bands are those of the 
state under consideration. The av fermion operators ft^v ^'^"^ ^q.av act onto subspaces spanned by neutral states. 
However, creation of one av fermion is a well-defined process whose generator is the product of an operator that fulfills 
small changes in the av effective lattice and corresponding av momentum band and the operator ft appropriate to 
the excited-state subspace. Since in the case of Eqs. (j60p and (|6ip it is assumed the av momentum bands are those 
of the state under consideration the corresponding generators on the vacua are simple products of ft^.^ operators, as 
given in these equations. 

As discussed in Appendix B, the use of the exact solution confirms that for ID the states ([60|) and ([61]) are momentum 
and energy eigenstates. Following the way that our operator description is constructed they are momentum eigenstates 
of the model on the square lattice as well. This is consistent with general properties that play a key role in the model 
physics: Both for the ID and square lattices and in contrast to the Hamiltonian ([1}, the momentum operator commutes 
with the set of av translation generators q^^ of El- (E3) and unitary operator V associated with the electron - rotated- 
electron unitary transformation. Since the Hamiltonian does not commute with the av translation generators q^vi 
the av fermion operators labelled by the av band discrete momentum values q^ where j = 1, ^1,^^ act onto and are 
defined in subspaces with constant values of A^^ , and the values of their discrete momenta are subspace dependent, 
the microscopic momenta qj of such objects are not in general conserved. Consistently, the states of Eqs. (j60p and 
(j6ip generated by momentum occupancy configurations of such av fermions are not in general energy eigenstates. 

The unitary operator V'^ considered in Subsection II-A has been constructed to inherently generating ex- 
act ?7/4t > energy and momentum eigenstates \^ lwS;U j^t) — V'^|^'nyS;oo) and |5'(7/4t) = V'^l^'oo) = 
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Ylaii^a)^''' ^^^^ / ^/^all"^ LWS;U/it) ■ For the model on the square lattice the energy and momentum eigenstates |5'c//4t) 
are within our description a well-defined superposition of momentum eigenstates \^u/4t) oi Eq. (|61l) of equal 5,,, S^, 
Ss, 5*1, and Sc values and momentum eigenvalue. The energy eigenvalues of the former states are U /At dependent. In 
turn, the invariance of the momentum operator under V imposes that the momentum eigenvalues of both the states 
\^u/it) and |$;7/4t) are independent of U/At. Therefore, provided that the corresponding set {|^'oo)} of U/At >• 1 
energy and momentum eigenstates and thus that {|$oo)} of momentum eigenstates associated with the same unitary 
operator are suitably chosen according to the recipe of Subsection II-A, the momentum eigenvalues of such states 
are the same as those of the corresponding U/At > momentum and energy eigenstates |^'c//4t) = V^^l^'oo) and 

momentum eigenstates |$[//4t) = ^^^|*&oo), respectively. 

Hence the rotated electrons used as starting building blocks of the description of the quantum problem in terms of 
c fermions, av fermions, and independent 77-spinons and spinous have been constructed to inherently the momentum 
eigenvalues being for U/At > fully determined by those of the corresponding U/At — >■ 00 system. For both the 
model on the square and ID lattices it follows from Eqs. I^E^ and (|27p that the momentum eigenvalues are indeed 
independent of U/At and given by, 

Cs Crj 

P = Y.'il ^-(^^ + QNsAQ) + mil'? - 9] N^A'l) + ^ A^r,,-l/2 ■ (63) 

q v=l q l/=l q 

Here Nc{q) and Nav{q) are the expectation values of the momentum distribution- function operators ([28]) . For the 
model on the ID lattice the c and av fermion discrete momentum values are good quantum numbers so that the 
momentum distribution functions Nc{q) and Na^iq) are eigenvalues of the operators p8)) and have values 1 and for 
occupied and unoccupied momentum values, respectivelyi^i^. As justified below, for the model on the square lattice 
that remains true for Nc{q) and for such a model in the one- and two-electron subspace defined in Section V this is 
true for Nsi{q) as well. 

For D = 1 the validity of expression ((63|) is confirmed by the exact solutionii^. Let us provide further information 
on why for the Hubbard model on the square lattice the momentum operator can for U/At > be expressed as given 
in Eqs. ([26| and ([27| so that the momentum eigenvalues are of the general form (|63l) . alike for the ID model. For 
U/t — > 00 electron single and doubly occupancy become good quantum numbers so that the electrons that singly 
occupy sites do not feel the on-site repulsion. That implies that the discrete momentum values q associated with 
their occupancy configurations separate into (i) those associated with electron hopping, which are good quantum 
numbers and involve the charge degrees of freedom and finite kinetic energy and are carried by the c fermions and (ii) 
those associated with the 77-spinon and spinous, including the composite rji/ and sv fermion occupancy configurations, 
respectively, whose momentum bands are for U/t — >■ 00 dispersionless. 

That for the model on the square lattice the momentum eigenvalues have for U /At > the general form, 

P = J2 + Pr,-sp,n + Psp^n , (64) 

q 

where Pr,-spin and Pspm are the momentum contributions associated with the ?7-spin and spin configurations, re- 
spectively, which for J7/i — >■ oo do not contribute to the kinetic energy, is an exact result. Indeed, in that limit 
the rotated electrons are the electrons, which for the finite-energy physics can only singly occupy sites and due to 
the exclusion principle behave as spinless and yy-spinless fermions. Since the only effect of the interactions on their 
charge degrees of freedom is to impose such a single occupancy constrain, those lead to a momentum contribu- 
tion '^„^^qNc{q) = '^qQ^io) where N{q) = N^{q) is the J7/t — ^ oo electron momentum distribution and 

Naicf) = (*oo|ct Cg,<T|*oo)- In that limit one then has Nc{q) = N{q). For U/At finite the U/t oo distribution N{q) 
is replaced by the rotated-electron momentum distribution, 

NroM = E(*C//4t|5L^?vl*C//4t> = Y.('^-o\n^^cl,VV^Cq^,aV]VM'^oo) = J] (*oo [ct .C,v I *oo) , (65) 

and the c fermion momentum distribution reads Nc{q) — Nrot{q)- Indeed, for rotated electrons singly occupancy 
remains a good quantum number for U /At > and the c effective lattice occupied sites have exactly the same spatial 
variables as the sites singly occupied by rotated electrons in the original lattice. 

For the model on the square lattice an energy eigenstate |$c//4t) is a superposition of 5,,, S*,^, Ss, S^, Sc, and 
momentum eigenstates \^u/4t) of Eq. (1611) with the same momentum eigenvalue, S'^, S*,^, Ss, S^, Sc values, and 
c fermion momentum distribution function Nc{q) = Nrotiq)- Hence all states |'&;7/4t) of such a set have the same 
distribution function Nc{q) — Y^^i^u/itl^'q cr^q,t7\^u /4t) ■ As discussed below, consistently with the commutator of 
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the Hamiltonian with the av translation generators q^jj of Eq. (|27p not vanishing, such states difFer in general in 
the av fermion momentum occupancies. Since for momentum eigenstates |$;7/4t) = V^^|*&oo) of Eq. (|6ip (one for 
each value of [//4i > 0) belonging to the same V tower the momentum eigenvalues are independent of U j^t, the 
contribution to the momentum of the c fermions is '^^qNcictj where alike in Eq. (|65p the distribution function 

Nc(,q) = Nrot{q) = Y.a'\'^u /it\c\^„c^j^a\'^u/it) has for U/At > the same magnitude N{q) = E<t(*oo|4,'t'^'?''^I^°°^ 
for the corresponding U/At — > oo initial momentum eigenstate |$oo)- 

This holds both for the model on the ID and square lattices. For the model on the former lattice the electron 
momentum distribution N{q) is for the x = and m — ground state plotted for C//t — >■ oo in Fig. 3 (a) of Rcf."*^. 
Note that the energy spectrum depends on U/At so that for the model on the square lattice the a; > and m ~ 
ground states may belong to different V towers for different values of U /At. The exception is the a; = and to = 
ground state, which belongs to the same V tower for all U/At > values, as confirmed in Ref."'^. 

The above analysis confirms that consistently with as given in Eq. (|29l) the Hamiltoninan and momentum operator 
commute with the c translation generator of Eq. (P7|) . for U/At > the distribution Nc{q) is an eigenvalue of the 
operator Nc{q) of Eq. ([SS)) and has values 1 and for filled and unfilled momentum values, respectively. That result 
combined with the validity of the momentum spectrum (I64[) for all occupancy configurations of the c fermions that 
generate the momentum eigenstates of form (j6ip confirms that the discrete momentum values q of the c fermion band 
are good quantum numbers for the Hubbard model on the square lattice. 

In turn, for U /At — > oo all 77-spin and spin configurations are degenerate and do not contribute to the kinetic energy. 
Nevertheless they lead to overall momentum contributions P^-spm and Pspim respectively. The main point is that 
the av fermion operators are defined in and act onto subspaces spanned by mutually neutral states for which for the 
square-lattice model the av translation generators cja^xi and qau x2 commute. That implies that alike in ID the av 
translation generators q^^ of Eq. (P7)) in the presence of the fictitious magnetic field B^i^ of Eq. (P5|) commute with 
the momentum operator both in the U/At ^ 00 limit and for U/At finite, as given in Eq. (I29[) . Consistently with the 
momentum operator expression provided in Eqs. (|26p and (j27p and as given in Eq. (j63p . that then implies that both 
in the U/At — > 00 limit and for U/At finite the momentum contributions Pr,-spin and Pspin of Eq. (j64p have the form 
Pjj-spin = J2u J2q[^ ~ <?] ^v^io) + ttM^ -1/2 and Pspin = J2u J2qQ^si^{q), respectively. 

For the Hubbard model on the square lattice the momentum spectrum of Eq. is found iiii to lead to the correct 
spin spectrum in some limiting cases. Note however that the momentum area N^^^ [^'''/i]^ and discrete momentum 
values number of the av bands are known. In contrast, the shape of their boundary remains an open issue. The 
shape of the c band is the same as that of the first Brillouin zone, yet the discrete momentum values may have small 
overall shifts under transitions between subspaces with different [Brj + Bg] = J^au ^<^'^ values. The studies of Ref- 
use several approximations to obtain information on the c and si bands momentum values and the corresponding c 
and si fermion energy dispersions of the model on the square lattice in the one- and two-electron subspace, for which 
as justified below the si fermion discrete momentum values are good quantum numbers. 

The square-lattice quantum liquid studied in this paper corresponds to the Hubbard model on the square lattice 
in such a subspace. For that quantum problem the c and si fermions play the role of "quasiparticles" . There are 
three main differences relative to an isotropic Fermi liquid^: i) Concerning the charge degrees of freedom, the non- 
interacting limit of the theory refers to At^ /U — )■ rather than to the limit of zero interaction J7 — ?► 0; ii) In the 
At^ /U — )■ limit the c fermions and si fermions become the holes of the "quasicharges" of RefJ^ and spin-singlet two- 
spin configurations of the spins of such a reference rather than electrons and only the charge dynamical factor becomes 
that of non-interacting spinless fermions, the one-electron and spin spectral distributions remaining non-trivial; iii) 
For U/At > the si band is full and displays a single hole for initial to = ground states and their one-electron excited 
states, respectively, and as found in Ref.^ its boundary line is anisotropic, what is behind anomalous one-electron 
scattering properties^. 

2. The good and quasi-good quantum numbers 

That Sri, S^, Ss, and are good quantum numbers for both the Hubbard model on the ID and square lattices 
implies that the numbers of independent 77-spinons (a — rj) and independent spinons (a — s) La,±i/2 — [Sa T S^] of 
Eq. are good quantum numbers as well. That the eigenvalue Sc of the generator of the new U (1) symmetry of the 
model global SO{3) x 50(3) x U{1) = [SU{2) x S';7(2) x [/(1)]/Z| symmetry found in Ref.^ is a good quantum number 
implies that the numbers N//*^ — Nc = Mg = 2Sc and = N^} = — [N^ — 2Sc] are also good quantum numbers. 
It then follows that the numbers Ca of Eqs. ([201), ^T^, dST]), and (H?]), which can be expressed as Ca = [N^^/2 - Sa], 
are good quantum numbers as well. The same then applies to the total numbers Ma^±i/2 = ^q,±i/2 + C'a of ±1/2 
7y-spinons (a — rj) and ±1/2 spinons (a = s). 
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The integrability of the ID Hubbard model is for Na oo associated with an infinite number of conservation laws 
such that the Hamiltonian commutes with the infinite av translation generators q^^^ of Eq. (j27p in the presence of 
the fictitious magnetic fields B^^. According to the results of Ref4§. such laws are equivalent to the independent 
conservation of the corresponding set of numbers {Na^} of av fermions, which are good quantum numbers. The lack 
of integrability of the Hubbard model on the square lattice is behind the av translation generators not commuting 
with the Hamiltonian and the set of numbers {Nai,} not being in general good quantum numbers. It follows that for 
the ID model the numbers Ba of Eqs. (|44|) and ([47|) and related numbers N^^-^ and iV^i are good quantum numbers as 
well, whereas for the model on the square lattice are quasi-good quantum numbers. By that we mean that in general 
they are not good quantum numbers for that model, yet they are good quantum numbers for it defined in the following 
subspaces: 

A) B^ and N^i (and Bg, iV^i, and iV^i) are good quantum numbers for the model on the square lattice in subspaces 
spanned by S",, = (and Ss = 0) momentum eigenstates provided that N^i^ > is finite for a single r]v branch (and 
Nsi, > is finite for the si branch and a single siy branch other than it) and iV^t/' = (and A^si^' = 0) for the remaining 
ijiy' (and siy') branches and B^ = N^^ (and [B^ - iV^i] = N^^) is such that 2vNr,^/N^ (and 2vNs,y/Nl/Nl 
and N^^_^/N^ is finite) for — > oo. In such subspaces the Hamiltonian commutes with the rjiy translation generators 
qrjiyxi and qriiyx2 in the presence of the corresponding fictitious magnetic field Bjji, (and si translation generators qsixi 
and qsix2 plus su translation generators qsi^xi and qsiyx2 in the presence of the corresponding fictitious magnetic 
fields Bsi and Bgu, respectively). In these subspaces the Nap av ^ si fermions correspond to a zero-momentum and 
vanishing-energy occupancy and are invariant under the electron - rotated-electron unitary transformation. Indeed, 
they obey the criterion of Eq. (|4ip . When all au 7^ si fermions of a given branch have that invariance their number 
Nai, is conserved so that in the present case the numbers B^ and N^^ where a = 77, s and Nsi if a s can be 
expressed solely in terms of good quantum numbers so that they are good quantum numbers as well. 

B) In subspaces where -/Va,,„ = "^^n (and N"^^^ = 28 s) for all branches with v > 1 (and v > 2) then N^i, = and 
NI^^j = 2Sr) (and Ngu = and A^^i^ — 2S's) are good quantum numbers for all such branches. In addition, B^^ = (and 
Bs — [Sc — Ss], N^^_^ = [Sc + Ss], Nsi — [Sc — Ss], and N^^ = 2Ss) arc good quantum numbers as well. This follows 
straightforwardly from Sr^ Ss, and Sc being conserved. In such subspaces the Hamiltonian commutes with the si 
translation generators qsixi and qsix2 in the presence of the corresponding fictitious magnetic field Bsi- 

The Hubbard model on a square-lattice in a subspace (A) and or in a subspace (B) provided that Ss = has in the 
limit U /At — ^ 00 an interesting physics. For such subspaces the av fermion occupancy configurations are in that limit 
and within a mean-field approximation for the fictitious magnetic field Ba,y of Eq. (j25p closely related to the physics 
of a full lowest Landau level of the 2D quantum Hall effect (QHE). Consider for instance a subspace (A) with Ss = 0, 
Nsi and Nsu finite where v =^ 1, Nsv' = for all remaining su' branches, and 2vNsv/N'^/N'^ — ^ and N^^_^/N^ 
finite for N^ 00. From the use of Eqs. dH]) and (gH]) one then finds that N^^^ = Nsi[l + 0(l/iVf )] = Nsi and 
^Osi, = for N^ — 00. That corresponds to the case when {np.^si) ~ 1 and {np-^su) ~ 1 and if in the mean field 
approximation one replaces the corresponding fictitious magnetic fields by the average fields created by all si and sv 
fermions, respectively, at a given position one finds as above Bsi{fj) w [$o/asi] f^xs and Bsu{rj) ~ ['i'o/aL] ^^s- The 
numbers of si and si^ band discrete momentum values can then be written as, 

A similar relation N"^^^ — Bjji, _L^/<f>o is found for subspaces (A) when there is an occupied rjiy band and for the si 

fermions alone for subspaces (B) and Ss — 0. It is found in Ref.^ that the aiy fermions have a momentum dependent 
energy dispersion and only in the limit U/it — 00 their energy bandwidth vanishes. Hence only in that limit are 
the A^a^i one-sl-fermion states and N^^^ one-si/-fermion states of Eq. (|66p degenerate. It follows that in such a limit 
N^^_^ = Bsi L^/^o and N^_^^ = Bs^L^/^o play the role of the number of degenerate states in each Landau level of 
the 2D QHE. In the subspaces considered here the ai^ fermion occupancies correspond to a full lowest Landau level 
with filling factor one, as given in Eq. (j66p . Only for the U/At — > 00 limit there is fully equivalence between the av 
fermion occupancy configurations and, within a mean-field approximation for the fictitious magnetic field B^u of Eq. 

the 2D QHE with a = 1 full lowest Landau level. In spite of the lack of state degeneracy emerging upon 
decreasing the value of U/4t, there remains for finite U/At values some relation to the 2D QHE. The occurrence of 
QHE-type behavior in the square-lattice quantum liquid shows that a magnetic field is not essential to the 2D QHE 
physics. Indeed, here the fictitious magnetic fields arise from expressing the effects of the electronic correlations in 
terms of the av fermion interactions. 
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The mean-field analysis associated with Eq. (pS)) is consistent with in the subspaces (A) and (B) the square- 
lattice model Hamiltonian commuting with the av translation generators g^.^, in the U/At — > oo limit so that their 
eigenvalues are good quantum numbers. Indeed, in that limit such an analysis refers to an effective description where 
the Hamiltonian is the sum of a c fermion kinetic-energy term and a QHE like Hamiltonian for each av branch. The 
av translation generators (f^^ then commute with all such Hamiltonian terms so that in the U/At — ^ oo limit their 
eigenvalues are indeed good quantum numbers for the model on the square lattice in a subspace (A) or (B). Since 
the electron - rotated-electron transformation is unitary such commutation relations also hold for U/At > so that in 
such subspaces the av fermion discrete momentum values are indeed good quantum numbers for U/At finite as well. 

In summary, for the Hubbard model on the square lattice each C//4t > energy eigenstate \'^u/4t) — V^^l^oo) where 
V is the uniquely defined electron - rotated-electron unitary operator considered in this paper can be expressed as a 
superposition of a set of S,j, S^, Ss, S^, Sc, and momentum eigenstates |<I';7/4t) = V^^|$oo) of general form given in 
Eq. (j6ip with equal values for such physical quantities, c fermion momentum distribution function Nc{q) = Nrotio) 
of Eq. dSSD, and numbers Ca = Y^^ZivNa^ of Eqs. dH]), (EH), and (gT]) where a = 77, s but different values 

for the av fermion numbers {Nav}- However, in some limiting cases such as for the energy eigenstates that span the 
one- and two-electron subspace that set reduces to a single state \^u/it) so that \^u/it) = \^u/it)- bi turn, that 
for the ID model the whole set of av fermion numbers {Nav} are good quantum numbers and the corresponding av 
translation generators q^^, commute with the Hamiltonian is behind the equality Y^jj/it) = \^u/it) always holding. 



F. Number of si fermion holes of the elementary excitations of a; > and m = Q ground states 

The number N'^-^ of unoccupied sites of the si effective lattice equals that of si fermion holes in the si momentum 
band and plays an important role in the general energy spectrum of the square-lattice quantum liquid introduced 
in this paper and further studied in Refsjii^^. Indeed, for the one- and two-electron subspace it is a good quantum 
number. Here we address the problem of the N^^ expression in terms of the number N of electrons, which contains 
important physical information. For the index choice a = s, Eq. (|44p reads N^-^ = [N^ — 2Bs\ where = 2Sc is the 
number of sites of the spin effective lattice and Bs the number defined in Eqs. (|44p and (|T7|) . From the use of expression 
Oq, = L/Na^ = [Na/Na^] tt given in Eq. (pO| where a = r],s, the number of unoccupied sites of the si effective lattice 
can be rewritten as N^i — [N^ — N^^ — 2Bs\ where N^^ = [N^ — 25c] is the number of sites of the jy-spin effective 
lattice. From combination of the expressions given in Eqs. and one finds N^^ = [Nj^ — N + 2C^ + 2L,, _i/2] 
where is given in Eq. (j3ip and L^ _i/2 is the number of independent —1/2 77-spinons. The use of that result in 
the above general 7Vj\ expression leads to, 

N':^ = [N~ 2Bs - 2Cri - 2L^,_i/2] . (67) 

For the model on the square lattice this number labels the momentum eigenstates of Eq. (j6ip but in general is not a 
good quantum number because the number Bg on the right-hand side of Eq. ([67)) given in Eqs. (|44|) and (|47)) is not 
in general conserved. (All remaining numbers on its right-hand side are.) 

Since all sites of the si effective lattice of x > and m — ground states are occupied and hence there are no 
unoccupied sites, from the use of such a N^^ expression one finds that for the corresponding excited states the number 
of unoccupied sites of the si effective lattice and thus of si fermion holes in the si momentum band reads, 

7Vi\ = [(5iV-2(5iV,i-2^Ar,, -2C^-2L„^_i/2] (general), 

1^=2 

N^^ ^ [6N - 26Nsi - 2Ns2] (one- and two - electron subspace) . (68) 

Here SN and SNsi are the electron and si fermion number deviations, respectively, generated in the transitions from 
the above ground states to the excited states. Since a; > and m = ground states have no 2z^-spinon composite sv 
fermions with ly > 2 spinon pairs, no 2i/-r/-spinon composite 771^ fermions, and no independent —1/2 ?7-spinons, the 
deviations in the corresponding numbers Ng^, C,i — ^^lii^ N^ju, and Lt),-i/2 equal the corresponding excited-state 
numbers. The expressions given in Eq. (j68|) refer to the latter numbers. 

The first expression is valid for the whole Hilbert space of the Hubbard model and hence takes into account all 
available excited states. It reveals that the number of holes in the si momentum band can be written as N!,\ = 
d[N — N^{^] where the absolute value of the number iVf/^ = [6N — N^^] given in Eq. (l68l) is always an even number. 
Therefore, for SN = ±1 one-electron excited states Nj}^ must be always an odd integer. In turn, for SN = and 
SN = ±2 two-electron excited states N^i must be always an even integer. Hence one-electron (and two-electron) 
excitations do not couple to excited states with two holes (and a single hole) in the si momentum band. 



33 



The second expression provided in Eq. ([BS]) refers to the one- and two-electron subspace defined in the ensuing 
section. The states that span such a subspace are of type (A) or (B) so that for them all numbers involved in that 
expression are good quantum numbers for the model on the square lattice, as confirmed below. 

V. THE SQUARE-LATTICE QUANTUM LIQUID: A TWO-COMPONENT FLUID OF CHARGE c 
FERMIONS AND SPIN-NEUTRAL TWO-SPINON si FERMIONS 

In this section we introduce and define the one- and two-electron subspace and corresponding spin and si effective 
lattices, whereas the c and si momentum bands associated with the c and si effective lattices, respectively, are studied 
in Ref.- . Moreover, here we address the problem of the generation of the one- and two-electron spectral weight in 
terms of processes of c fermions and spinous. The picture which emerges is that of a two-component quantum liquid 
of charge c fermions and spin neutral two-spinon si fermions. For the model on the square lattice it refers to the 
square-lattice quantum liquid further studied in RefSfii^^. 

A. The one- and two-electron subspace as defined in tiiis paper 

Let I^Gs) be the exact ground state for x > and m = and O denote an one- or two-electron operator. Then 
the state, 

6|*G5) = E^.l*^); C, = {^,\d\^Gs). (69) 
j 

generated by application of O onto that ground state is contained in the general one- and two-electron subspace. 
This is the subspace spanned by the set of energy eigenstates such that the corresponding coefficients Cj are not 
vanishing. 

Our goal is finding a set of excited states that have nearly the whole spectral weight J2j P ~ 1 of the above one- 
and two-electron excitations. For a a; > (and a; = 0) and m = ground state I^'gs) it is found here that such states 
have excitation energy a; < 2/i (and uj < /i*'). The inequality w < 2/i applies to hole concentrations < x < 1 and 
m — 0. In turn, for the ^ = and m — absolute ground state with the chemical-potential zero level at the middle 
of the Mott-Hubbard gap the smallest energy required for creation of both one rotated-electron and one rotated-hole 
doubly occupied site is instead given by the energy scale This justifies why for such initial ground state the above 
inequality w < 2/i is replaced by a; < fjP. 

From the use of the invariance under the electron - rotated-electron unitary transformation of the independent ±1/2 
spinous, one finds that the number of such objects generated by application of A/'-electron operators onto a ground 
state is exactly restricted to the following range, 

Ls - [Ls, -1/2 + is. +1/2] = 2Ss = 0, 1, 2, . (70) 

It follows that for the one- and two-electron subspace only the values Lg = [Lg^ _i/2 + Lg, +1/2] — 2S's ~ 0, 1, 2 are 
allowed. Such a restriction is exact for both the model on the square and ID lattices. 

An important property is that for a number ^ > 2 of spinon pairs the siy fermions created onto a a; > and 
771 = ground state have vanishing energy and momentum. Since at 777 = one has that H = 0, such objects obey 
the criterion egi, = 2vpLB H = oi Eq. (|4ip . Hence they are invariant under the electron - rotated-electron unitary 
transformation. Therefore, for U /At > they correspond to the same occupancy configurations in terms of both 
rotated electrons and electrons. That reveals that such si/ fermions describe the spin degrees of freedom of a number 
27/ of electrons. It follows that nearly the whole spectral weight generated by application onto the ground state of 
A/'-electron operators refers to a subspace spanned by energy eigenstates with numbers in the following range, 

[is,-i/2 + is,+i/2 + 2C,-2i?,] =0,1,2,...,AA. (71) 

Note that owing to the above invariance of the sv fermions with ly > 2 spinon pairs under consideration provided 
that A/'/A^jP — > and Bg/N//' — > for N/^ — > 00 the number Bg = Ngi, is a good quantum number. (This is a 
generalization of the subspace (A) considered in Subsection IV-D.) Consistently, the x > (and a: = and fJ, = 0) and 
777 = ground state and the set of excited states of energy cj < 2/i (and uj < ijP) that span the one- and two-electron 
subspace introduced in this paper have no —1/2 7/-spinons, rjiy fermions, and siy' fermions with v' > S spinon pairs so 
that Njjiy = and Ng^i = for i/' > 3. As further discussed below, summation over the states that span that subspace 
as defined here gives J2j IQP ~ 1 ^o^' the coefficients of the one- or two-electron excitation 0\^gs} of Eq. and 
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TABLE I: The deviations 5N1^ = ~25Sc and numbers iVi\ = [2Ss + 2Ns2] of Eq. ((TS]) for the fourteen classes of one- and 
two-electron excited states of the x > and m — ground state that span the one- and two-electron subspace as defined 
in this paper and Ref.— , corresponding electron number deviations 5Nf and 5N_i, and independent-spinon numbers I;s,+i/2 
and Ls,-i/2 and s2 fermion numbers Ns2 restricted to the values provided in Eq. (|72|) . The spin Ss and deviations SSc, 
SNsi = [SSc — 5*8 — 2A'^s2], and SNa^i ~ [SSc + Ss] of each excitation are also provided. 



there is both for the model on the ID and square lattices an extremely small weight corresponding mostly to states 
with Ns3 = 1, which is neglected within our definition of the one- and two-electron subspace. 

Thus, according to Eq. ((7T|) the numbers of independent ±1/2 spinous and that of s2 fermions of the excited states 
that span such a subspace are restricted to the following ranges, 

Ls.±i/2 =0,1; A^,2 = 0, for AA-1, 
2Ss + 2Ns2 = [i,, _i/2+i,, +1/2 + 2A^,2] =0,1,2, for = 2 . (72) 

Here = 1,2 refers to the A/'-electron operators O whose application onto the ground state |^gs) generates the 
above excited states, as given in Eq. ([M)) . Furthermore, Nc = N = {1 — x) and A^si = [N/2 — 2Ns2 — Sg] = 
(1 — x) /2 — [2Ns2 + Ss]. If in addition we restrict our considerations to the LWS-subspace of the one- and two- 
electron subspace, then _i/2 = in Eq. (|72|l . whereas the values Ls,+i/2 = 0,1 for Ns2 = and J\f ~ 1 remain 
vahd and in [2Ss + 2Ns2] = 0, 1, 2 one has that 2Ss = is,+i/2 for TV = 2. 

We recall that the numbers of independent ?7-spinons {a = rj) and independent spinous {a = s) ^q,±i/2j total 
numbers of 77-spinons [a = 77) and spinous {a = s) Ma,±i/2 — [La, ±1/2 + C'a], number of sites of the spin effective 
lattice iV^ — 2Sc, number of sites of the ry-spin effective lattice iV^ — [iV^ — 2iS'c], number of c fermions Nc — 2Sc, 
and number of c fermion holes — [iV^ — 2Sc] are good quantum numbers of both the Hubbard model on the ID 
and square lattices. The good news is that in the one- and two-electron subspace as defined here the numbers N^^_^ , 
Nsi, N^i, and Ns2 are good quantum numbers of the Hubbard model on the square lattice as well. Indeed, for spin 
Ss = that subspace is a subspace (A) as defined in Subsection IV-D so that Nsi = [Sc — 2Ns2], iV^i = 2Ns2, and Ns2 
and hence N^^^ = [Ngi + N^-^] = Sc are good quantum numbers. Furthermore, for the remaining spin values Ss = 1/2 
and Ss = I such a subspace is a subspace (B) as defined in that subsection so that Nsi = [Sc — Ss], N^^ — 2Ss, 
and N^^_^ = [Sc + Sg] are good quantum numbers. That implies that the microscopic momenta q carried by the si 
fermions are good quantum numbers. 

The number of sites, unoccupied sites, and occupied sites of the c and si effective lattices equal those of discrete 
momentum values, unfilled momentum values, and filled momentum values of the c and si bands, respectively. Use 
of the general expressions of si band discrete momentum values and number N^i of si band unfilled momentum 
values given in Eq. (|44|) for a = s together with the restrictions in the values of the numbers of Eqs. ((TO)) . (|7T|) . and 
([7^ and the result proved above in Subection IV-E that one-electron (and two-electron) excitations have no overlap 
with excited states with none and two (and one) si band holes (and hole) reveals that nearly the whole one- and 
two-electron spectral weight is contained in the subspace spanned by states whose deviation SNc in the number of c 
band holes and number N^i of si band holes are given by, 

dN'c' = -2SSc = -SN = 0,Tl,T'2, 

N^^ = 2Ss + 2Ns2 = ±{SN^ - 6N^) + 2Ls,^i/2 + 2iV,2 = Ls, _i/2 + i.,+1/2 + 2iV,2 = 0, 1, 2 . (73) 

Here SN is the deviation in the number of electrons, SNf and SN^ those in the number of spin-projection f and J, 
electrons, respectively, Ns2 the number of the excited-state s2 fermions, and ^^,±1/2 that of independent spinous 
of spin projection ±1/2. The deviations SN^} and numbers N^^ of Eq. ([75]) for the fourteen classes of one- and 
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two-electron excited states of the a; > and m = ground state that span the one- and two-electron subspace, 
corresponding electron number deviations SN-^ and SN^, and independent-spinon numbers and is^_i/2 and 

s2 fermion numbers A^s2 restricted to the values provided in Eq. (|72|) are given in Table HI For excited states with 
Ns2 = 1 the s2 effective lattice has a single site and the corresponding s2 band a single vanishing discrete momentum 
value, q = 0, occupied by the s2 fermion. Such a s2 fermion is invariant under the electron - rotated-electron unitary 
transformation and has zero energy, consistently with the invariance condition of Eq. (j4ip for ai' — s2 and zero 
magnetic field H = 0. 

The initial x > and m — ground states of the one- and two-electron subspace have zero holes in the si band so 
that SN^i — NI^i for the excited states. The one- and two-electron subspace as defined here is spanned by the states of 
Table U generated by creation or annihilation of \SN^\ = 0, 1, 2 holes in the c momentum band and N^^ = 0, 1, 2 holes 
in the si band plus small momentum and low energy particle-hole processes in the c band. The charge excitations of 
a; > and m = initial ground states consist of a single particle-hole process in the c band of arbitrary momentum 
and energy compatible with its momentum and energy bandwidths, plus small-momentum and low-energy c fermion 
particle-hole processes. Such charge excitations correspond to state representations of the global U{1) symmetry and 
refer to the type of states denoted by "charge" in the table. The one-electron spin-doublet excitations correspond to the 
four types of states denoted by "±lael." in Table [T] where -1-1 and —1 denotes creation and annihilation, respectively, 
and a =t,i- The spin-singlet and spin-triplet excitations refer to the four types of states denoted by "singl.spin" 
and "tripl.spin" in the table. The two-electron excitations whose electrons are in a spin-singlet configuration and 
those whose two created or annihilated electrons are in a spin-triplet configuration correspond to the five types of 
states "±2 fiel." and "±2(Tel." of that table where +2 and —2 denotes creation and annihilation, respectively, of two 
electrons. 

Alike for the ID model, for the model on the square lattice such fourteen types of states are energy eigenstates and 
contain nearly the whole one- and two-electron spectral weight. Excited states of classes other than those of the table 
contain nearly no one- and two-electron spectral weight. Such a weight analysis applies to the ID Hubbard model as 
well. Moreover, for the quantum liquid corresponding to the Hamiltonian ([T]) in the one- and two-electron subspace 
as defined in this paper, the numbers 2Sc, 25^, 2Ss, and —28^ associated with the global SO{3) x SO{3) x U{1) 
symmetry are given by, 

2Sc ^ (l-x) ; 2S,-, = -2S'^ ^ X ■ 2Ss ^ {l~x) ~2[Nsi+2Ns2]; -25,^ = m iVf = 25,-2L,,_i/2 • (74) 
For such a quantum liquid the number of sites and lattice spacing of the spin effective lattice read, 

< = (l-x)iVi^; a,^ ' (1 - a;) > 1/iVf . (75) 

We emphasize that for the model in the one- and two-electron subspace the concept of a 77-spin lattice has no 
physical significance because the 77-spin degrees of freedom refer to the Cr, = {N^ /2 ~ Sc — S,-,) =0 vacuum |0^; Nj^^) 
of Eq. (1^51) . which for S'^ > corresponds to a single occupancy configuration of the A^^P -1-1/2 ?7-spinons. Only when 

both N^JN^ = [1 - '^.Sc/N^] > and = (iVf /2 - Sc - S^) > is the concept of a ?7-spin effective lattice well 
defined and useful. Indeed, for that subspace the numbers and of sites of the original lattice whose electron 
occupancy configurations that generate the energy and momentum eigenstates are not invariant under the electron 
- rotated-electron unitary transformation are given by = and Cs = [Ngi + 2Ns2], respectively. For x = (and 
X = 1) one finds iV^ = and N° = (and N^^ = and N^^ =0), so that there is no jy-spin (and no spin) 
effective lattice and the spin (and ?7-spin) effective lattice equals the original lattice. 

As further confirmed below in Subsection V-D, for the one- and two-electron subspace defined here only the c and si 
fermions have an active role. Straightforward manipulations of the above general expressions given in Eqs. (|42l) - ((53)) 
reveal that for that subspace the number A^^^^ of si band discrete momentum values, A^^i of si fermions, and N!^^ of 
si fermion holes are given by, 

nF^, = [^^.1 + Ni\] = [Sc + Ss] ; A^.i = [Sc -Ss- 2N,2] ; N^i = [25, + 2N,2] = 0, 1, 2 , (76) 
respectively. In turn, the corresponding c effective lattice, c momentum band, and c fermion numbers read, 

N^^^[N, + N'^^]^N^; N, = 2S, = {l-x)N^; ^ x . (77) 

As mentioned above, for states belonging to the one- and two-electron subspace the 77-spin degrees of freedom 
correspond to a single occupancy configuration referring to the ?7-spin vacuum \Orj.N^^) of Eq. ([23]) (and to no 
configuration since iVf = in that vacuum). Therefore, when defined in that subspace the Hamiltonian (jlj has an 
effective global U{2)/Z2 = 50(3) x U{1) symmetry. This is anyway the global symmetry relevant to more general 
models belonging to the same universality class as that considered here. Indeed, the 77-spin SU{2) symmetry disappears 
if one adds to the Hamiltonian ([l} kinetic-energy terms involving electron hopping beyond nearest neighbors. In 
contrast, the global 50(3) x U{\) symmetry is robust under the addition of such terms. 
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B. Generation of the one- and two-electron spectral weight in terms of processes of c fermions and spinons 

Expression by use of the operator relations given in Eq. (|17|) of the operator O corresponding to the one- and 
two-electron operators of the excitations 0\^gs) of Eq. (|69|) in terms of the of c fermion and rotated quasi-spin 
operators reveals that the elementary processes associated with the number value ranges of Eqs. (fTTj) . ([72l) . and (|73|) 
are for the one- and two-electron subspace fully generated by the operator O. 

In turn, expressing by means of the operator relations provided in Eq. (jl7[) the operator terms of Eq. ([3]) containing 
commutators involving the operator S — [t/U) [T+i — T-i\+ 0{t^ /U'^) in terms of the of c fermion and rotated quasi- 
spin operators and taking into account that independently of their form, the additional operator terms 0{t^ /U^) of 
higher order are products of the kinetic operators Tq, T+i, and T_i of Eq. one finds that the processes generated 
by such operator terms refer only to excitations whose number value ranges are different from those of Eqs. (j7ip . 
(17^ . and ((75)) for the one- and two-electron subspace. That confirms that such processes generate very little one- and 
two-electron spectral weight, consistently with the discussions and analysis of Subsections II-A and V-A. 

For the model on the ID lattice the spectral-weight distributions can be calculated explicitly by the pseudofermion 
dynamical theory associated with the model exact solution, exact diagonalization of small chains, and other methods. 
The relative one-electron spectral weight generated by different types of microscopic processes is studied by means 
of such methods in YieiM-. The results of that reference confirm the dominance of the processes associated with the 
number value ranges provided in Eqs. ([7^ and (175)) . In this case the general operators O and O are electron and 
rotated-electron, respectively, creation or annihilation operators. The operator O generates all processes associated 
with the number value ranges of Eqs. (|7^ and (|75)) and number values of Table |I| In addition, it also generates some 
of the non-dominant processes. That is confirmed by the weights given in Table 1 of Ref*^, which correspond to the 
dominant processes associated with only these ranges. The small missing weight refers to excitations whose number 
value ranges are not those of Eqs. (|72)) and (|73| but whose weight is also generated by the operator O. Indeed, that 
table refers to U /4t — > cx) so that O = O and the operator terms of Eq. ([3]) containing commutators involving the 
operator 5* vanish because all such commutators vanish in that limit. 

On the other hand, for finite values of U/it all dominant processes associated with the number value ranges of Eqs. 
(l72t and (|73)) and number values provided in TableUare generated by the operator O whereas the small spectral weight 
associated with excitations whose number value ranges are different from those are generated both by that operator 
and the operator terms of Eq. ^ containing commutators involving the operator S. For the model on the ID lattice 
the small one-electron spectral weight generated by the non-dominant processes is largest at half filling and f7/4i « 1. 
For the range of hole concentrations a: > considered in this paper the one- and two-electron subspace is spanned by 
states with vanishing rotated-electron double occupancy. Its generalization for the range a; < reveals that then it 
is spanned by states with vanishing rotated-hole double occupancy. That property combined with the particle-hole 
symmetry explicit at x = and ^ = 0, implies that the relative spectral-weight contributions from different types of 
excitations given in Fig. 2 of Ref.— for half-filling one-electron addition leads to similar corresponding relative weights 
for half-filling one-electron removal. Analysis of that figure confirms that for the corresponding one-electron removal 
spectrum the dominant processes associated with the number value ranges of Eqs. (|72p and (1731) and number values 
given in Table U refer to the states called in the figure 1 holon - 1 si hole states, whose minimum relative weight of 
about 0.95 is reached at U /At ~ 1. For other hole concentrations a; > and values of U/At the relative weight of such 
states is always larger than 0.95, as confirmed from analysis of Figs. 1 and 2 and Table 1 of Ref4i. 

For the Hubbard model on the square lattice the explicit derivation of one- and two-electron spectral weights is 
a more involved problem. The number value ranges of Eqs. (I71|) . ()72[) . and ()73p and number values provided in 
Table U for the one- and two-electron subspace also apply, implying similar results for the relative spectral weights of 
one- and two-electron excitations. Furthermore and as discussed above, analysis of expression (1551) reveals that for 
6N = ±1 (and 6N — ±0 and 6N — ±2) excited states the number Nl\ of holes in the si band must be always an odd 
(and even) integer so that one-electron (and two-clcctron) excitations do not couple to excited states with two holes 
(and one hole) in the si band. Both such selection rules hold also for the one- and two-electron subspace for which 
iVj\ = [SN - 26N,i - 2Ns2l as given in Eq. 

C. The spin and si effective lattices for the one- and two-electron subspace 

Within the operator description introduced in this paper, the subspace with relevance for the one- and two-electron 
physics is the one- and two-electron subspace as defined above, which is spanned by states with no r]i/ fermions and 
no sf fermions with a number v > 3 oi spinon pairs. According to the restrictions and numbers values of Eqs. ()72p 
and ()73l) and TablelH such states may involve none or one s2 fermion. As confirmed in Refs.^'^^''^^, it is convenient to 
express the one- and two-electron excitation spectrum relative to initial x > and to = ground states in terms of the 
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deviations in the numbers of c effective lattice and si effective lattice unoccupied sites given explicitly in Eq. (|73p and 
TablelH Note that the number of si fcrmions provided in Eq. (|76p can be written as A^^i = [(1 — x)N^ /2 — Sg — 2Ns2] 
where Sg ~ for Ns2 = 1 and Sg = 0, 1/2, 1 for Ns2 ~ 0. The square-lattice quantum liquid associated with the 
Hubbard model on the square lattice in the one- and two-electron subspace is further studied in Refsfii^^. The studies 
of Ref.^'^ refer to the same quantum liquid perturbed by small 3D anisotropy associated with weak plane coupling. 

Since for Ns2 = 1 states the s2 fermion has vanishing energy and momentum and in addition consistently with Eq. 
()4ip for vanishing magnetic field _ff = is invariant under the unitary transformation associated with the operator V, 
the only effect of its creation and annihilation is on the numbers of occupied and unoccupied sites of the si effective 
lattice and thus in the occupancies of the discrete momentum values of the si band, as discussed below. Therefore, 
for the study of the Hamiltonian ([Ij in the one- and two-electron subspace the only composite object whose internal 
occupancy configurations in the spin effective lattice are important for the physics is the two-spinon si fermion and 
associated si bond particle studied in Ref.'^®. 

It is found below in Subsection V-D that for the model in the one- and two-electron subspace as defined in Subsection 
V-A the presence of independent -1-1/2 spinous or a composite s2 fermion is felt through the numbers of occupied 
and unoccupied sites of the si effective lattice, whereas the number of independent -1-1/2 77-spinons equals that of the 
unoccupied sites of the c effective lattice. Therefore, when acting onto that subspace, the Hubbard model refers to a 
two-component quantum liquid described in terms of c fcrmions and si fermions. 

For ^ 1 the description introduced in Ref.^S of both the spinon occupancy configurations of a si bond-particle 
associated with a local si fermion and of the underlying si effective lattice is a good approximation provided that 
the ratio N^jNa is finite. Since the number of sites of the spin effective lattice is given by = (1 — x)N^ , that 
requirement is met provided that the electronic density n = (1 — x) is finite. Such a description is exact for hole 
concentrations Q < x 1 when the spin-effective lattice becomes the original lattice. In that limit all sites of the 
spin effective lattice become sites of the original square or ID lattice. In turn, for the model on the square lattice 
it does not apply when x ~ 1 and there is a finite number of sites in the spin effective lattice only. So it is a good 
approximation for the range (1 — a;) > and the limit 3> 1 considered in our studies. Particularly, it is a 

good approximation for hole concentrations x € (0,a;*) considered in Refsii^^, where a;* G (0.23,0.28) for the range 
U/At G {uo,ui) of interest for the studies of these references where uq « 1.3 and ui « 1.6, respectively. 

For a; < 1 the spin effective lattice has a number of sites = {1 — x) smaller than that of the original 
lattice. Nevertheless for ^ 1 a site of the spin effective lattice with real-space coordinate fj in that lattice has 
a well-defined position r ~ fj in the original real space contained within the system area or length where D — 2 
and D — 1, respectively. This is consistent with the lattice constant Os defined in Eq. ([30]) for a = s, which for the 
one- and two-electron subspace reads, 

a;< 1, I? - 1,2, (78) 



' (l-x)i/c 

being such that the area {D = 2) or length {D — 1) of the system is preserved. Obviously, any real-space point 
within the spin effective lattice corresponds to the same real-space point in the system original lattice. So the real- 
space coordinates rji of the j' — 1, sites of the spin effective lattice are for a; > different from the real-space 
coordinates rj of the j — 1,...,N^ sites of the original lattice. The vectors fj and rj' refer though to well-defined 
positions in the same ID or 2D real space contained within such a lattice. 

The suitability of the description introduced in this paper in terms of occupancies of effective lattices is confirmed 
for ID by the results of Appendix B: The discrete momentum values of the c fermions and au fermions correspond to 
quantum numbers of the model exact solution. In turn, no exact solution exists for the model on the square lattice. 
However, the discrete momentum values obtained for such objects refer to good quantum numbers for that model in 
the one- and two-electron subspace. 

For such a ^ 1 subspace there is commcnsurability between the real-space distributions of the A^^^ « Ngi 
sites of the si effective lattice and the N^^ w 2Nsi sites of the spin effective lattice. For {1 — x) > l/N^ and ^ 1 
the spin effective lattice has A^^ = (1^2;) sites and from the use of the expression given in Eq. ([75)1 for A^^^ and 
Eq. ([5^ for Usi we find, 

- a. « a. (l - J^^^^] - a. , 23. = 0,1,2. (79) 

For N!}i = states the bipartite ID and square spin effective lattices have two well-defined sub-lattices. As 
discussed in RefJ^, for Nj}^ — 1,2 states the spin effective lattice has two bipartite lattices as well, with one or 
two extra sites corresponding to suitable boundary conditions. For the model on the square lattice the two spin 
effective sub-lattices have lattice constant Ugi = v^fls- In turn, for ID the sites of each spin effective sub-lattice are 
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distributed alternately along the chain, the corresponding nearest-neighboring sites being separated by Ugi = 2as. 
The fundamental translation vectors of such sub-lattices read, 

Osi = flsi exi for ID ; a^i = -^i^^i + ^^2) , b^i = -^{e^^ - e^^) for 2D , (80) 

where and are the unit vectors and xi and X2 Cartesian coordinates. As confirmed in RefJ^, the vectors given 
in Eq. ([50]) are the fundamental translation vectors of the si effective lattice. 



D. The two-component quantum liquid of charge c fermions and spin-neutral two-spinon si fermions 

In the case of a; > 0, to = 0, and Nl^i = ground states of the model on the square lattice whose si momentum 
band is full and all N^^_^ sites of the si effective lattice are occupied we consider that the square root Na^ of the 
number N^^ of sites of the spin effective lattice is an integer so that the spin effective lattice is a square lattice with 
Na^ X Na^ sites. It follows that the square root Na^^ of the number iV^_^^ of sites of the si effective lattice cannot in 
general be an integer number yet N^^_^ is. However, within the ^ 1 limit considered here and as mention above 
for general av branches, we use the notation N^^_^ for the number of sites of the si effective lattice where D = 1 and 
£) = 2 for the model on the ID and square lattice, respectively. 

For the one- and two-electron subspace as defined in this paper the number A^^^ of sites of the si effective lattice 
and si band discrete momentum values, Nsi of si fermions, and N^^ of si fermion holes have expressions given in Eq. 
([75)1 . The corresponding numbers for the c effective lattice and c band are provided in Eq. ([77|. For that subspace 
the si band is either full by Ngi — N^_^_^ — 2Sc — A^/2 si fermions or has one or two holes. Furthermore, one- 
electron and two-electron excitations have no overlap with excited states with two holes and one hole in the si band, 
respectively. Hence as given in Table |T] excited states with a single hole in the si band correspond to one-electron 
excitations and those with N^^ = 0,2 holes in that band refer to two-electron excitations. Excited states with Nj}^ = 3 
(and NI^i = 4) holes in the si momentum band correspond to very little one-electron (and two-electron) spectral 
weight and are ignored within our definition of one- and two-electron subspace. 

According to the number value ranges of Eqs. (j72p and (j73p and number values provided in Table HI the one- and 
two-electron subspace is spanned by excited states having either none Ns2 = or one A^s2 = 1 spin-neutral four-spinon 
s2 fermion. When A^s2 = 1 there are no independent spinous and one finds for the Ns2 = 1 states spanning that 
subspace that N^^^ = 1 for A^s2 = 1 so that the s2 fermion occupies a s2 band with a single vanishing momentum value 
and hence A'^^2 = holes. Indeed, when Ns2 — 1 the s2 fermion has vanishing energy and momentum and consistently 
with Eq. (j41[) is for H = invariant under the electron - rotated-electron unitary transformation. Therefore, the only 
effect of its creation and annihilation is on the numbers of occupied and unoccupied sites of the si effective lattice 
and corresponding numbers of si band si fermions and si fermion holes. Specifically, according to the expressions 
provided in Eq. (|76|) and number values of Table U the deviations SSc = SSg = and 6Ns2 = 1 generated by a state 
transition involving creation of one s2 fermion lead to deviations in the number of si fermions and si fermion holes 
given by SN^i = ~26Ns2 = -2 and SN^^ = 25Ns2 = 2, respectively. 

It follows that for A^s2 — 1 excited states one has from the ranges of Eqs. ([72]) and ((73)) and number values of Table 
Uthat 5*5 = so that according to Eq. (1751) the number of holes of the si band is given by Afj\ = 2Ns2 = 2, in contrast 
to A^^i = for the initial ground state. This is consistent with what was mentioned previously: Upon creation onto 
the ground state of a s2 fermion the number A^i^^ of sites of the si effective lattice remains unaltered and following 
the annihilation of two si fermions and creation of one s2 fermion two unoccupied sites appear in it and hence two 
holes emerge in the si band. The emergence of these unoccupied sites and holes involves two virtual processes where 
(i) two si fermions are annihilated and four independent spinous are created and (ii) the latter independent spinous 
are annihilated and the s2 fermion is created. 

The only net effect of creation of the vanishing-energy and zero- momentum s2 fermion is the annihilation of two si 
fermions and corresponding emergence of two holes in the si band and two unoccupied sites in the si effective lattice. 
Therefore, in the case of the one- and two-electron subspace one can ignore that object in the theory provided that the 
corresponding changes in the si band and si effective lattice occupancies are accounted for. For excited states with a 
single s2 fermion two of the four spinous of such an object are used within neutral si fermion particle-hole excitations 
in the motion of si fermions around in the si effective lattice as unoccupied sites of that latticei^, consistently with 
the expression N^^ — 2Ns2 given in Eq. (|76p . 

A similar situation occurs for the Ls — 2Ss independent spinous, whose number Ls ~ [Ls.-i/2 + ^^,-1-1/2] — '^Ss 
belongs in the present subspace to the ranges given in Eqs. d72]) and (173)) and Table U such that 25*8 = 0,1,2. For 
2Ss = 1,2 one has that A^s2 = so that following Eq. (|75| the number of unoccupied sites of the si effective lattice 
and of si fermion holes is for the excited states under consideration given by A^^^ = 25'^ = 1, 2. However, in contrast 
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to the s2 fermion, a deviation S2Ss ~ 1,2 generated by a state transition may lead to deviations in the numbers 
of occupied and unoccupied sites of the si effective lattice and corresponding si fermion and si fermion holes that 
do not obey the usual equality SN^i = -5N^^. Indeed, now 26 Sc = ±1 for SN^^ = 26 Ss = 1 and 26 Sc = 0,±2 for 
6N!}i = 26Ss = 2 and according to the expressions provided in Eq. ([76|) such deviations lead to deviations in the 
numbers of occupied and unoccupied sites of the si effective lattice and corresponding numbers of si fermions and 
si fermion holes given by ^A^si = [6Sc — 6Ss] and 6N!}i = 62Ss, respectively. It follows that the total number of 
sites and thus of discrete momentum values of the si band may also change. That leads to an additional deviation 
6N^^_^ = [6Sc + 6Ss]- As given in TableHl for one-electron excited states one has that 6Nl\ — 26Ss — 1 and 26Sc = ±1 
so that 6Nsi ±1/2 - 1/2 = -1,0 and 6N^^^ =±1/2+1/2 = 0, -1. For Ns2 = two-electron excited states one 
has (5A^i\ = 2,55'^ = 2 and 26Sc = 0, ±2 so that'^iV^i = -1, ±1 - 1 = -2, -1, and JTV^^ = 1, ±1 + 1 = 0, 1, 2. 

Excitations that involve changes 6N^^_^ = [6Sc + 6Ss] in the number of sites and discrete momentum values of the si 
effective lattice and si band, respectively, correspond for the si fermion operators and f^^si to transitions between 
different quantum problems. Indeed, such operators act onto subspaces spanned by neutral states, which conserve 
Sc, Ss, and N^^_^. In turn, the generator of non- neutral excitations is the product of an operator that makes small 
changes associated with the small variations 6N^_^_^ — [6Sc + 6Ss] both in the si effective lattice and corresponding si 
momentum band with a si fermion operator or a product of such operators appropriate to the excited-state subspace. 

Also the vanishing momentum and energy -/j,j.+i/2 — xN^ independent +1/2 77-spinons are invariant under the 
electron - rotated-electron unitary transformation and have no direct effect on the physics. For a; > and the one- 
and two-electron subspace considered here they correspond to a single occupancy configuration associated with the 
ry-spin vacuum |0,;;7V^) of Eq. (^5]) . In turn, the degrees of freedom of the rotated-electron occupancies of such 

xNj^ sites of the original lattice associated with the U{1) symmetry refer to the N'^ = x unoccupied sites of the 
c effective lattice of Eq. (l77l) and corresponding c band holes. Therefore, alike for the one- and two-electron subspace 
the effects of creation and annihilation of a s2 fermion or independent spinous are taken into account through the 
corresponding changes in the numbers of occupied and unoccupied sites of the si effective lattice and corresponding 
si fermion and si fermion hole numbers, the creation and annihilation of independent +1/2 ?7-spinons has no effects 
on the physics except that their number equals that of unoccupied sites of the c effective lattice and c band holes. 

In summary, when acting onto the one- and two-electron subspace as defined in Subsection V-A, the Hubbard model 
refers to a two-component quantum liquid described in terms of two types of objects on the corresponding effective 
lattices and momentum bands: The charge c fermions and spin-neutral two-spinon si fermions. For the model on 
the square lattice that is the square-lattice quantum liquid further studied in Refsiii^^. The presence of independent 
+1/2 spinous or a composite s2 fermion is taken into account through the numbers of occupied and unoccupied sites 
of the si effective lattice and corresponding si fermion and si fermion holes, whereas the number of independent 
+1/2 77-spinons equals that of the unoccupied sites of the c effective lattice and c band holes. Otherwise, the presence 
of vanishing momentum and energy independent +1/2 spinous, spin-neutral four-spinon s2 fermion, and independent 
+1/2 ?7-spinons has no effects on the physics. This property follows from all such objects being invariant under the 
electron - rotated-electron unitary transformation. 

Spin-singlet excitations generated by application onto a m = and x > initial ground state of the operator 
fo s2 fq,si fq'.si where q and q' are the momenta of the two emerging si fermion holes are neutral states which 
conserve Sc, Ss, and N^^_^. (See Table HI) For the model on the square lattice the role of the s2 fermion creation 
operator /g ^2 is exactly canceling the contributions of the annihilation of the two si fermions of momenta q and q' 

to the commutator [qsixi,qsix2] of the si translation generators in the presence of the fictitious magnetic field Bsi 
of Eq. (j25p . so that the overall excitation is neutral. Since the s2 fermion has vanishing energy and momentum and 
the si band and its number N^^_^ of discrete momentum values remain unaltered, one can effectively consider that 
the generator of such an excitation is f^^si fq'.si and omit the s2 fermion creation, whose only role is assuring that 
the overall excitation is neutral and the two components of the si fermion microscopic momenta can be specified. It 
follows that for the one- and two-electron subspace the operators fq',si, /^/^^i fq',si, and /g-,^! /I, ..^ 

generate neutral excitations. 

The quantum-liquid c fermions are 77-spinless and spinless objects without internal degrees of freedom and structure 
whose effective lattice is identical to the original lattice. For the complete set of energy eigenstates that span the 
Hilbert space the occupied sites (and unoccupied sites) of the c effective lattice correspond to those singly occupied 
(and doubly occupied and unoccupied) by the rotated electrons. The corresponding c band has the same shape and 
momentum area as the first Brillouin zone. 

In contrast, the quantum-liquid composite spin-neutral two-spinon si fermions have internal structure and the 
definition of the si effective lattice in terms of both the original lattice and the spin effective lattice as well as the 
spinon occupancy configurations that describe such objects is for the one- and two-electron subspace a more complex 
problem, which deserves and requires further studies^*. It is simplified by the property that the si fermion occupancy 
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configurations of the states that span the one- and two-electron subspace refer to a si effective lattice with none, one, 
or two unoccupied sites. Also the shape of the corresponding si momentum band and corresponding boundary and 
the form of the c and si energy dispersions are complex unsolved problems investigated in Refj^. 

VI. LONG-RANGE ANTIFERROMAGNETIC ORDER AND SHORT-RANGE SPIN ORDER OF THE 
MODEL ON THE SQUARE LATTICE FOR a; = AND < a; < a;., RESPECTIVELY 

Here we profit from the rotated-electron description introduced in this paper to study the occurrence at zero 
temperature of a long-range antiferromagnetic order and a short-range spin order in the Hubbard model on the square 
lattice at X = and for a well defined range of finite hole concentrations, respectively. For small hole concentrations 
< X <C 1 and intermediate and large U /At values the latter order has a spiral-incommensurate character—. 

A. Extension of the Mermin and Wagner Theorem to the half- filling Hubbard model for (7/4t > 

We recall that for U /At ^ 1 the spin degrees of freedom of the half-filling Hubbard model on a square lattice are 
described by a spin-1/2 isotropic Heinsenberg model on a square lattice. It follows that the Mermin and Wagner 
Theorem^® is valid for the former model at half filling and U/At ^ 1. The theorem states that then there is no 
long-range antiferromagnetic order for finite temperatures T > 0. 

Let us provide evidence that the Mermin and Wagner Theorem applies to the half-filling Hubbard model on a 
square lattice for all values U/At > 0. The possibility of such an extension to C//4t > is strongly suggested by 
evidence involving the transformation laws of the spin configurations under the electron - rotated-electron unitary 
transformation. We recall that in terms of rotated electrons, the occupancy configurations that generate the energy 
eigenstates are the same for all finite values U/At > 0. Moreover, such rotated-electron occupancy configurations 
equal those that generate such states in terms of electrons in the U/At ^ oo limit. 

However, the electronic occupancy configurations that in the U/At ^ cjo limit generate the energy eigenstates that 
span the one- and two-electron subspace associated with our description correspond to an involved problem. Indeed, 
the suitable objects whose occupancy configurations that generate such states have a simple form, rather than the 
electrons, are the "quasicharges" , spins, and pseudospins of Ref ji^,, which are our c fermion holes, spinous, and r]- 
spinons, respectively, in that limit. For the model on the square lattice in the one- and two-electron subspace the 
energy bandwidths of the av fermion dispersions vanish for a; > in the U/At — > oo limit and the c fermion dispersion 

has for < a; < 1 the simple form edq) ~ ^'^i'l2i=i i'^^^ilxi) ~ cos{qFcxi)] hi terms of the c band momentum 
components- . 

The use of the following two properties provides below strong evidence that the Mermin and Wagner Theorem 
holds for the half-filling Hubbard model on the square lattice for U /At > 0: I) The x = and m — absolute ground 
state is invariant under the electron - rotated-electron unitary transformation^ so that the occurrence of long-range 
antiferromagnetic order for U/At — ^ oo, associated with that of the spin-1/2 isotropic Heisenberg model, implies the 
occurrence of that long-range order for U/At > as well; II) Since in terms of rotated electrons for U/At > single 
and double occupancy are good quantum numbers, the rotated-electron occupancy configurations that generate the 
energy eigenstates are more ordered than the corresponding electron occupancy configurations. It follows that the 
lack of long-range antiferromagnetic order in terms of the spins of the rotated electrons implies as well a lack of such 
an order in terms of the spins of the electrons whose occupancy configurations generate the same states. 

We recall that the rotated-electron occupancy configurations that for U /At > generate the energy eigenstates 
are identical to the electron occupancy configurations that for U/At oo generate the energy eigenstates belonging 
to the same V tower. Hence concerning the original-lattice rotated-electron occupancies, the Mermin and Wagner 
Theorem applies to all finite values of U/At > 0: for the occupancy configurations of the rotated-electron spins there 
is no long-range antiferromagnetic order for temperatures T > 0. That the lack of long-range antiferromagnetic 
order for U /At > and T > of the rotated-electron spins implies a similar lack of such an order for the spins of the 
original electrons follows from the above property II, being consistent with for finite values of U /At the spin occupancy 
configurations being more ordered for the rotated electrons than for the electrons. 

This is also consistent with the existence of long-range antiferromagnetic order for the Hubbard model on a square 
lattice at hole concentration x = 0, temperature T = 0, and U/At > only, associated with that of the a; = 
and m = absolute ground state, which the above property I refers to, in agreement with the numerical results of 
Refs3-'-P. As discussed below, there is strong evidence that for x = 0, T > 0, and U/At > such an order is replaced 
by a short-range spin order. 

Since the electron - rotated-electrom unitary operator V does not commute with the Hamiltonian, the rotated- 
electron occupancy configurations that generate the energy eingenstates are for C//4i > the same as for U/At ^ oo 
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in terms of electrons. However the corresponding energies depend on the value of U /At. It turns out that the rotated- 
electron occupancy configurations associated with the long-range antiferromagnetic order refer to a subspace that for 
U/At finite exists below an energy given by the energy scale Eq. (|34l) . Such a subspace is spanned by energy 
eigenstates with iV^^ = [7V^ — 2Sc\ = and N'^^ = 2Sc = sites so that there is no Ty-spin effective lattice and the 
spin effective lattice has as many sites as the original lattice and is identical to it. In such a subspace there are no 
rotated-electron doubly occupied sites and no rotated-electron unoccupied sites. 

The point is that the number of sites of the spin effective lattice being given by — is a necessary condition 
for the occurrence of long-range antiferromagnetic order: Concerning the spins of the rotated-electron occupancy 
configurations such an order occurs when the ?7-spin effective lattice does not exist and the quantum-system vacuum 
(Hg becomes \0s) = |0^; N^) x IGS^; iVf ) and the system global symmetry is U{2)/Z2 = SO{2,) x U{1). For U/At oo 
when the spin degrees of freedom of the half-filling Hubbard model are described by the Heisenberg model, the subspace 
for which N'^ = describes the whole finite-energy physics. This follows from according to Eq. the energy 
scale being given by /i" = C//2 — > oo. However, for finite U/At > one has that the spin occupancy configurations 
in a spin effective lattice with N^^ = sites, which are expected to be behind the long-range antiferromagnetic 
order, exist only for energies below The Mott- Hubbard gap 2fjP, which refers to the charge degrees of freedom, 
affects the spin degrees of freedom as well: /z*' is the energy above which the ry-spin effective lattice emerges and hence 
the number of sites of the spin effective lattice decreases so that < N^. It follows that for the Hubbard model 
on a square lattice fjP is the energy below which the long-range antiferromagnetic order of the rotated-electron spins 
survives for a; = 0, m = 0, and temperature T = 0. 

Finally, it follows from the above property II that since the ground-state rotated-electron occupancy configurations 
are for finite values of U/At more ordered than those of the electrons for the same state, a lack of long-range anti- 
ferromagnetic order of the rotated-electron spins for N^^ < implies a similar lack of long-range antiferromagnetic 
order for the spins of the original electrons. In contrast to the x = and m = absolute ground state of property I, 
< X < 1 and m = ground states are not invariant under the electron - rotated-electron unitary transformation. 



B. Short-range spin order of the model on the square lattice for m = and < x <^ 1 

The requirement for the occurrence of long-range antiferromagnetic order at zero temperature T = in the Hubbard 
model on the square lattice that the number of sites of the 77-spin and spin effective lattices are given by N^^ = 

and = N"^, respectively, implies that there is no such a long-range order for m = 0, a; > 0, and T — Q. Indeed, 
then 2Sc < Nl so that Nl^ = [Nl - 2Sc] > and N^^ = 2Sc < N^- For the one- and two-electron subspace such 
expressions read N^^ = x > and N^^ = {1 — x) < N^, respectively. 

One may use as zero-temperature order parameter of the m — and a; = quantum phase the excitation energy /i*^ 
below which the long-range antiferromagnetic order survives, whose limiting behaviors for U/At <^ 1 and U/At^ 1 
are given in Eq. (j34l) for D = 2. In addition to those, here we provide an approximate formula valid for intermediate 
U/At values uq < U/At < ui, 

« 32te-'^^ , C//4t < 1, 

2e^ t 

« + (U/At - uo) , uo<U/At<ui, 

TT 

w [U/2 - At] ; U/At > 1 , (81) 

where the behavior ^q{U /At) « ^0(^0)^/1 + {U/At — uq) is expected to be a good approximation for the small range 
uo ^ U/At < ui, uo ~ 1.302, and ui sa 1.600. The approximate magnitude /io(uo) ~ [2e^/7r] t is that consistent with 
the relation /lo(uo) ~ fJ-oiu*) / y/ 1 + (m* — uq) where uq = 1.302 is the U /At value at which the energy parameter Aq is 
found below to reach its maximum magnitude and = 1.525 is that for which the studies of Ref.— lead by a completely 
different method to fJ,o{u^) « 566 meV for t — 295 meV and U/At k, u^, — 1.525. The use of /io('Uo) ~ [2e^ /tt] t in the 
formula ^oiU/At) ~ /xo(mo)\/1 + (U/At — uq) leads for t ~ 295 meV, U/At « = 1.525, and uq — 1.302 to nearly the 
same magnitude, /zo(m*) ~ 565 meV. 

There is strong evidence of the occurrence in the half-filling Hubbard model on the square lattice of strong short- 
range antiferromagnetic correlations for T > and below a crossover temperature called in Ref.^^, which here 
we denote by Tg*. This is consistent with then the system being driven into a phase with short-range spin order. 
Furthermore, that the occurrence of long-range antiferromagnetic order requires T — 0, = 0, and = is 
consistent with the short-range spin order occurring for m = 0, 0<a;<?;l, and < T < Tq being similar to that 
occurring for m = 0, x = 0, and < T < Tq* and was studied previously in Ref.— for < T ^ Tq . 
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One may choose the excitation energy 2Ao below which the short-range spin order with strong antiferromagnetic 
correlations survives to play the role of zero-temperature order parameter of the quantum phase displaying such an 
order occurring for m — and < a; ^ 1. That energy parameter has a U/At dependence qualitatively similar to 
that of the energy scale 2kB Tq, with the equality 2Ao ~ 2kB Tq holding approximately. For both vanishing and finite 
temperatures a phase displaying a short-range spiral-incommensurate spin order occurs for (i)m = 0,0<a::<Cl, 
and < T ^ Tg* and (ii) m = 0, x = 0, and < T <C Tg*. At m = and temperatures below Tq the system is 
driven both for (i) < a; ^ 1 and < T < Tq and (ii) x = and < T < Tq into a renormalized classical regime 
where the x — and T = long-range antiferromagnetic order is replaced by such a short-range spin order, which is 
a quasi-long-range spin order as that studied in Refj^ for simpler spin systems. 

Let us study the U/At dependence of Ag sa kg Tq by combining the results obtained from the use of the general 
description introduced in this paper with those of the low-temperature approach to the half-filled model of Ref.^-, 
whose investigations focus on temperatures < T ^ Tg*. We use often in the following the magnitudes of the 
energy bandwidth W^^ and energy parameter derived in Refii For the value U/it « 1.525 found in Ref.'^^ to be 
appropriate to the description of the properties of several hole-doped cuprates on combining the description introduced 
in this paper with results on the spin- wave dispersion of the half-filled Hubbard model in a spin-density- wave-broken 
symmetry ground state obtained by summing up an infinite set of ladder diagrams'"^. 

The picture that emerges is that for the Hubbard model on a square lattice at T = there is an overall zero- 
temperature energy order parameter 2|A|, which at a; = has a singular behavior due to a sharp quantum phase 
transition occurring at vanishing hole concentration, such that there is long-range antiferromagnetic order at a; = and, 
as justified below by combining our description with the results of Refst^ '^'^'^'' , a short-range spiral-incommensurate 
spin order for < a; <C 1. Consistently, the energy scale 2|A| reads, 

2|A|=Ai°, x = 0; 2|A| w 2Ao (^1 - , < x < 1 , U/At > uq ~ 1.302 . (82) 

The linear dependence of 2|A| « 2Ag(l — x/x^) on x where x^ ~ 2rs/n and = Ag/4W°^ is justified in Ref.^. Hence 
2|A||2,^Q = is for U/4t > different from 2Ag = limj.^Q 2|A|. Symmetry arguments imply that lim(7/4t_^g 2Aq = 
2|A||a;=o = M*^ so that the two energy scales become the same in that limit. 

The energy scale Ag is found in Refsii^ to play an important role in the square-lattice quantum-liquid physics. 

Upon increasing J7/4t, it interpolates between Ag « ^°/2 « 16^6"''^^ for U/At < 1 and Ag « AW^^ « tt [UY /U 
for U/At ^ 19, going through a maximum magnitude at a U/At value found below to be approximately given by 
U /At ^ uq ^ 1.302. The energy parameter /io/2 is an increasing function of U/At. As given in Eq. (|8T1) . it behaves as 

//g/2 « 16 1 e^'^V^ for U/At < 1 and as /io/2 « [U/2~At] for U/At > 1. In turn, the energy scale AW^^ is a decreasing 
function U/At, which according to Eq. (|38l) is for 13 = 2 given by AW^i = 16t at U/At and decreases approximately 
as AW^i ~ TT [4t]^/[/ for very large U/At. The energy scale Ag interpolates between these two behaviors, vanishing 

both for U/At and U/At — >• oo as Ag « Idte"^^/**^ and Ag « ■jT{At)^/U where it becomes /ig/2 and an energy 
scale 7r(4t)^/[/ associated with the strong < x ^ 1 antiferromagnetic correlations, respectively. This is consistent 
with its maximum magnitude being reached at the intermediate U/At value U/At = uq for which /io/2 ~ AW^i and 
we consider that such a U/At value is that at which the equality /ig(uo)/2 = 4M^]\(uo) holds. 

The ratio = Ao/4T4^f]^ is parametrized in the following as — e^^" where As = | ln(Ao/4T/V°]^)| has an important 
role in such an interpolation behavior. The energy scale Ag can then be expressed as, 

Ag = r, AW°^ = AW^^ e-^= ; A, - | ln(Ag/4W^o ) | , (83) 

where A^ has the limiting behaviors, 

A,, = TT^yAt/U, U/At<t:l; Xs~Atuo/U, uqq <U/At < m; Xs=0, U/At ->-oo, 
uoo ~ {uo/t^Y w 0.171 ; wg w 1.302 ; ui w 1.600 . (84) 

Consistently, the ratio Ts is an increasing function of U /At, which changes continuously from = for U/At — > to 
Ts — 1 ior U/At oo. That for ugg < U /At < ui it is approximately given by Vs ~ g-'^t '^o/u jj^gtead of Vs ~ e^'^'vA*/^ 
for U /At <C 1 is consistent with for large U/At it being such that (1 — Vg) oc At/U rather than (1 — r^) c>c y/At/U. 

The temperature which plays the role of Tg* w Ag/fc^ is plotted in Fig. 3 of Ref<^. Its U/At dependence 
is qualitatively correct, vanishing both in the limits U/At —i' and U/At ^ oo and going through a maximum 
magnitude at an intermediate value 5/4 < U/At < 3/2. Nevertheless, the interpolation function used to produce it, 
provided in Ref. 74 of such a paper, is poor for intermediate values of U/At. However, that does not affect the validity 
of the results of Ref.^^. Indeed, the studies of that reference refer to the temperature range < T ^ T^, for which 
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the accurate dependence of on J7/4t is not needed and the goal of its Fig. 3 is just illustrating qualitatively the 
behavior over the entire coupling range^^. 

Our goal is reaching a more quantitatively accurate U/4t dependence of the energy parameter Aq ~ ksT^ for 
intermediate U/^i values, consistent with the qualitative physical picture of Ref.— . That according to Eqs. ([83l) and 
([M)) one has that Aq = '^W^i e~^' for C//4i > where the parameter is a continuous decreasing function of U /At 
given by As = oo for [//4i — > and As = for U/At oo reveals that the magnitude As = 1 separates two physical 
regimes. Consistently, As = 1 refers to the U/At value U/At = uo at which Ao reaches its maximum magnitude 
max{Ao} = [moM/2] e-i = 4M^,"iK) 

That as found in Ref.^, both the parameter is for approximately mq < U/At < ui given by x'^ = ^Ts/tt = 
2e-4tiio/£^y'7r and at U/At w = 1.525 reads x^ « 0.27 implies that uq sa 1.3. Moreover, in that reference it 
is found that AW^^iu^) = [198.4/295] t w 0.673i at U/4t « = 1.525. Since W^^ = lim^_^ VFsi has the same 
magnitude as the a; = energy bandwidth Wsi, such a result is obtained in that reference from comparison of 
the a; = and m = spin excitation spectra for the high symmetry directions found by use of the description 
introduced in this paper with those estimated by the controlled approximation of Ref."^". On combining the equality 
/io(uo)/2 = 4VFs°i(wo) with the magnitude /io(uo) ~ [2e^/TT]t of Eq. §^ we find AW°^{uo) « [eV7r] t « 1.731 i and 
max{Ao} = [mo(uo)/2] e'^ = ^W^^iuo) e"! « t/n. 

Note that the magnitude Wgi{u^,) — [49.6/295] t « 0.168t obtained in Ref.^ is about 12.25 times smaller than that 
found by use of the limiting expression Airt'^/U « 2.060 f at U/At = 1.525. This reveals that for the Hubbard model 
on the square lattice such a limiting expression is valid for a smaller range of very large U/4t values than for ID. 
Specifically, we find that it is valid for U/U > 19. The value U/At sa 18.70 w 19 is that at which Ant'^/U is given by 
0.168 i, alike W^^ at U/U w 1.525. Therefore, the relation W^^ = J w int'^/U is indeed valid for U/it > 19 and for 
the model on the square lattice the energy scale J ~ Airt'^/U controls the physics for a smaller U/4t range than in 
ID, which corresponds to very large U/At ^ 19 values. Hence the intermediate-C//4i range plays a major role in the 
physics of that model, as confirmed by the related studies of Refs^^ '^^i^^ . 

For the model on the square lattice and approximately the intermediate range uq < U/it < ui the U/it dependence 
of the energy parameter AW^i is of the form AW^^ ~ [g^/tt] t W{U/4t) where W{u) is an unknown function of m = U/At 
such that W{uo) = 1. Since on combining Eqs. (gSl) and ^ one finds Ao = 4M^fi e^"**""/^, fulfillment of the 
maximum condition 9Ao(m)/9u = at u = u/At = uq requires the function W{u) be given by W{u) w (2 — u/uq) for 
< [{u-uo)/iui-uo)] < 1. In turn, the use of the above results 4:W°^{u^) = [198.4/295] t « 0.673 1 and AW°j^{uo) « 
[e^/n] t K, 1.731 1 reveals that the ratio AW°-y{u^,) / AW°^{uq) may be expressed as AW°^{uJ) / AW^-yiu^) « [1 - (m* - wq)] 
so that W{u) « [1 — (u — Wo)] for w « = 1.525. We then use a suitable interpolation function for W{u), which 
has these two limiting behaviors, so that the energy parameter 4^/"°^ is for U/4i intermediate values U/At £ {uq, ui) 
given approximately by. 

Mo < U/At < Ml , 

(85) 

The value mo ~ 1.302 is that obtained from the equation 4M^°i(m*) = [e^ t/n] [1 - (m* - mo)] = [198.4/295] t « 0.673 i. 
On combining the above results we find the following magnitudes for the energy parameter Ao, 

Ao w 16ie~'^V^, L//4t<l, 

= max{Ao} w t/TT , ?7/4t = Mo, 

« e(^-^tuo/u) wiU/it) , Mo < U/U < Ml , 

= 4iy°i w TT [4t] VC/ , U/At > 19 , (86) 

where W{U/4t) is the interpolation function given in Eq. Note that 9Ao(m)/9u = at m = u/At = uq, 

consistently with Ao reaching its maximum magnitude t/TT at that U/At value. The overall U/t dependence of 
Tq « Ao/fcs is similar to that plotted in Fig. 3 of Ref.^^ for with the U/t value at which the maximum magnitude 
is reached shifted from U/t « 5.60 to U/t « 5.21 and that magnitude lessened from maxIfc^Ta;} sa 0.625 1 to 
maxIfc^TQ*} « max{Ao} sa t/vr sa 0.318 1. Otherwise, such a dependence has qualitatively a shape similar to that 
illustrated in the figure. 

Next let us provide strong evidence that for intermediate and large values of U/At and small hole concentrations 
< X <C 1 the short-range spin order of the m = ground state corresponds to that of a spin-singlet incommensurate 
spiral state. In terms of the rotated-electron spins occupancy configurations the ground state is then a spin-singlet 
incommensurate spiral state for U/At > 0, m — 0, and < a; ^ 1. That evidence is found on combining the above 
result that such a ground state is a spin-singlet state with the recent results of Refi^ concerning the spin degrees of 
freedom of a related quantum problem. 



4M^°i «^ —WiU/U); VK(m) sa 1- (M-Mo)e~^^'"^""^ 

TT 

Mo « M* — 1 + e~^7r sa 1.302; M*sa 1.525. 

295 
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As discussed in Section II, for intermediate and large values of U/At the Hubbard model on the square lattice given 
in Eqs. ([1]) and ^ in terms of electron and rotated-electron operators, respectively, can for D = 2 be mapped onto 
an effective t — J model on a square lattice with t, t' = t'{U /At), and t" — t"{U /At) transfer integrals where the role 
of the processes associated with t' = t'{U/At) and t" — t"{U/At) becomes increasingly important upon decreasing the 
U/At value. Rigorous results on the spin degrees of freedom of the t — J model on a square lattice with t, t', and 
t" transfer integrals were recently achieved in Ref.— . The investigations of that paper refer to small values of the 
hole concentration < a; <C 1 and spin density m = and have as starting point a suitable action first introduced 
in Refi^. The use in Refi^ of a staggered CP^ representation for the spin degrees of freedom allows to resolve 
exactly the constraint against double occupancy, which is equivalent to performing the electron - rotated-electron 
unitary transformation. In order to achieve the rigorous result that for small hole concentrations there occurs a 
incommensurate-spiral spin order, the effective action for the spin degrees of freedom is reached after integrating 
out the charge fermionic degrees of freedom and the magnetic fast CP^ modes. Importantly, the dependence on the 
hole-concentration of the coupling constants of the effective field theory is obtained explicitly for small x. 

From the mapping between the above effective t — J model and the Hubbard model of Eqs. ([T]) and @ in the 
subspace with vanishing rotated-electron double-occupancy, the studies of Ref.~ reveal that for intermediate and large 
values of U/At, spin-density m — 0, and small hole concentrations < a; ^ 1 the ground state is an incommensurate 
spiral state. This is a rigorous result, yet the studies of Ref.^ are not conclusive on whether for < x ^ 1 the 
m = incommensurate spiral ground state has short-range or long-range spin order. However, symmetry requires 
that the occurrence of long-range and short-range incommensurate-spiral spin order is associated with a spin-singlet 
ground state and a broken-symmetry ground state without well defined spin, respectively. Therefore, combination 
of our above result that the m — and < x ^ 1 ground state of the Hubbard model ((T|) is a spin-singlet state 
with the result of Refill that such a state has a short-range incommensurate-spiral spin order is consistent with 
for intermediate and large values of U/At the ground state of the model being a spin-singlet state with short-range 
incommensurate-spiral spin order and thus strong antifcrromagnetic correlations. 

Moreover, in terms of rotated electrons the short-range incommensurate-spiral spin order prevails for U/At > 0, 
m = 0, and < x <C 1 and implies that in terms of electrons the system has a short-range spin order. Indeed and as 
discussed above, the ground-state rotated-electron occupancy configurations are for finite values of U /At more ordered 
than those of the electrons for the same state so that a lack of long-range spin order of the rotated-electron spins 
for TO = and < x ^ 1 implies a similar lack of long-range spin order for the spins of the original electrons. The 
occurrence of a short-range incommensurate-spiral spin order for < x <C 1 is fully consistent with the occurrence for 
U/At > of a long-range antifcrromagnetic order and a short-range spin order for x = and < x 1, respectively. 
Above the occurrence of such long-range and short-range spin orders is associated with the number of sites of the 
spin effective lattice obeying the relations = and < N^, respectively. In turn, here the occurrence for 
< X <SC 1 of a short-range spiral-incommensurate spin order is reached by use of completely different arguments. 

The energy order parameter 2|A| of Eq. (15^ is both for x = and x > identified in the studies of Refsjii^^ with 
the maximum pairing energy of the —1/2 and +1/2 spinous of a composite spin-neutral two-spinon si fermion. Its 
magnitude has a singular behavior at x = due to the sharp quantum phase transition from a ground state with 
long-range antifcrromagnetic order at x = to a state with short-range spiral-incommensurate spin order and strong 
antifcrromagnetic correlations for < x ^ 1. The results of Refsiii^^ extend the to = and T = short-range 
spin order to a well-defined range of hole concentrations < x < x*. According to these results, for U/At > uq the 
X dependence 2|A| w 2Ao(l — x/xj) given in Eq. (15^ for < x <C 1 is valid for x e (0, x,) whereas for x > x, 
the energy scale 2|A| vanishes. Here for approximately U/At e {uqjU-^) where u-^ > ui is the U/At value at which 

= 1/2 the critical hole concentration x, equals the [//4t-dependent parameter x'^. For U/At > uq the studies of the 
above references identify it with a critical hole concentration above which there is no short-range spin order so that 
2|A| — ^ as X — >■ X*. Indeed, the energy order parameter 2|A| vanishes in the absence of that order. Consistently, 
the spinon pairing energy vanishes for the Hubbard model on the ID lattice for which 2|A| — for the whole range 
of hole concentrations. The short-range incommensurate-spiral spin order discussed here for < x ^ 1 corresponds 
then to a limiting case of the general short-range spin order that according to the investigations of Refs.— occurs 
for < X < X,. 

The form of the si effective lattice spacing (|79l) is for the square-lattice model directly related to the above spin 
orders. Indeed, that at x = and m = the square si effective lattice has no unoccupied sites and its spacing is 
given by Csi = V^a reveals that then its periodicity has increased relative to that of both the original lattice and 
c effective lattice. This is because at x = the si effective lattice refers to a x reconstruction in which the 
periodicity of the spin-sub-system real-space structure is increased. Such an effect is consistent with the occurrence 
of the long-range antifcrromagnetic order at x = and to = 0. In turn, that for .x > and to = the square si 
effective lattice remains having no unoccupied sites but its spacing reads instead — y^2/{l — x) a is consistent 
with the emergence of the short-range incom mensurate- spiral spin order discussed here. Furthermore, that then the 
si effective-lattice spacing is given by a^i = ■\/2/(l — x) a implies that the si fermion occupancy-configuration states 
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break the translational symmetry. 

The studies of Refsjii^^ reveal that for the square-lattice quantum liquid the si bond particles associated with the si 
fermions are closely related to the spin-singlet bonds of Refj^. The investigations of Ref. — combine the description 
introduced in this paper with some of the results of Refj^ on a related problem to investigate the square-lattice 
quantum-liquid fluctuations and competing orders in the presence of a small 3D anisotropy perturbation associated 
with weak plane coupling. They confirm the validity of the long-range and short-range spin orders found here at 
a; = and a: > 0, respectively, by studying the fluctuations of the phases that control such orders. 

At a; = and m — the spin effective lattice is identical to the original lattice so that the two-spinon si bonds 
of the si fermions refer to the same lattice as in Ref^^^ and the si effective lattice is one of its two sub-lattices. 
Also for < X <C 1 the spin effective lattice is very similar to the original lattice and the si effective lattice spacing 
flsi w [\/2/(l + x/2)] a « \/2a is close to that of the original sub-lattices. As further discussed in Ref.— , at zero 
temperature, x = 0, and m = there are in the the square-lattice quantum liquid strong phase fluctuations whose 
action has a local compact gauge symmetry. That corresponds to the long-range antiferromagnetic phase studied 
above for which monopole-antimonopole pairs of the type considered in Refi^^ are unbound and proliferate. For small 
values < X ^ 1 the motion of the c fermions and the associated kinetic energy play the role of symmetry-breaking 
Higgs terms. The presence of such terms supresses free monopoles and is behind the replacement of the long-range 
antiferromagnetic order at x = by the short-range incommensurate-spiral spin order for < x ^ 1. 

VII. CONCLUDING REMARKS 

In this paper we have profited from the interplay of the global 50(3) x 50(3) x U{1) symmetry found in Ref.^ for the 
Hubbard model on any bipartite lattice with the transformation laws of a well-defined set of operators and quantum 
objects whose occupancy configurations generate a complete set of 5^, 5^, 5^, 5|, 5c, and momentum eigenstates 
under a suitably chosen electron - rotated-electron unitary transformation to introduce a general description for the 
model on a square lattice with ^ 1 sites. Often the problem is also addressed for the model on the ID lattice. 

Within the description introduced in this paper the above complete set of states is generated by occupancy con- 
figurations of c fermions associated with the state representations of the U{1) symmetry, r?-spin- neutral 2z/-ry-spinon 
composite riv bond particles and independent 7y-spinons associated with the state representations of the 7y-spin SU (2) 
symmetry, and spin-neutral 2i'-spinon composite sv bond particles and independent spinous associated with the state 
representations of the spin SU (2) symmetry. The index u = 1,2, ... refers to the number of 77-spinon or spinon pairs in 
each composite object. Evidence is provided in this paper that in addition to being 5,,, 5^, Ss, 5|, 5c, and momentum 
eigenstates, for the Hubbard model on the square lattice in the one- and two-electron subspace defined in this paper 
such states are energy eigenstates. (For the ID model they are energy eigenstates for the whole Hilbert space^.) 

The relation of the composite rji' bond particles (and sv bond particles) to the ?7-spinons (and spinous) has sim- 
ilarities with that of the composite physical particles to the quarks in chromodynamics^. Within the latter theory 
the quarks have color but all quark-composite physical particles are color-neutral. Here the ?7-spinon (and spinous) 
that are not invariant under the electron - rotated-electron unitary transformation have 77-spin 1/2 (and spin 1/2) 
but the 2i'-77-spinon (and 2i/-spinon) composite rjv bond particles (and su bond particles) are 77-spin-neutral (and 
spin-neutral) objects. The components of the microscopic momentum values of the av fermions are eigenvalues of av 
translation generators q^^, of Eq. (I?7|) in the presence of the fictitious magnetic fields Bai/{rj) of Eq. Those are 

associated with extended Jordan- Wigner transformations from which the au fermions emerge from the corresponding 
av bond particles. In the U /'it — >■ 00 limit the av band momentum occupancy configurations that generate the energy 
eigenstates of Hubbard model on a square lattice in subspaces (A) and (B) defined in this paper are described by 
the 2D QHE physics. The structure of the spin configurations associated with a si fermion simplifies in the one- 
and two-electron subspace introduced in this paper where only the charge c fermions and spin- neutral two-spinon si 
fermions play an active role. The physical picture that emerges is that of a two-component quantum liquid of charge 
c fermions and spin-neutral two-spinon si fermions whose momentum values are good quantum numbers. For the 
Hubbard model on a square lattice in the one- and two-electron subspace the composite si fermion consists of two 
spinous in a spin-singlet configuration plus an infinitely thin flux tube attached to it. Thus, each si fermion appears to 
carry a fictitious magnetic solenoid with it as it moves around in the si effective lattice. The square-lattice quantum 
liquid is further investigated in Refs.^'^^. 

Concerning previous studies on the large-C/ Hubbard model and t — J model on a square lattice involving for 
instance the slave particle formalis m^'^'^i'^^ or Jordan- Wigner transformational^, the crucial requirement is to impose 
the single occupancy constraint. Here that constraint is naturally implemented for all values of U/At > 0, since the 
spins associated with the spin SU{2) state representations refer to the rotated electrons of the singly occupied sites. 
Moreover, for the above schemes the spinless fermions arise from individual spins or spinous. In contrast, within the 
extended Jordan- Wigner transformation considered in this paper and Refii the si fermions emerge from spin- neutral 
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two-spinon composite si bond particles. Elsewhere the two-spinon occupancy configurations in the spin effective 
lattice of the si bond particles are studied in detail for the model ([T]) in the one- and two-electron subspaceiS.. This is 
a needed step for the investigations of Refsj^i^^ of the square- lattice quantum liquid of c and si fermions with residual 
interactions. Following the results of these references, it is expected that such a quantum liquid captures the universal 
properties of general many-electron models with short-range interactions on the square lattice. 

The general rotated-electron description introduced here for the Hubbard model on the square lattice is consistent 
with a ground state with long-range antiferromagnetic order for half filling, short-range incommensurate-spiral spin 
order for < x ^ 1, and short-range spin order for finite hole concentrations x below a C//4t-dependent critical value 
X* above which the system is driven into a disordered state without short-range spin order. Strong evidence is given 
in Refs i^^i^^ that for a well-defined hole-concentration range x e (xc,a;*) where Xc ~ 10~^ and S (0.23,0.31) for 
approximately Uq ^ U /At < u^^ a, long-range superconducting order coexists with the short-range spin order in the 
square-lattice quantum liquid perturbed by small 3D anisotropy associated with weak plane coupling. 

Finally, the results of Refsi^ '^^i'^^ suggest that the description introduced in this paper is useful for the further 
understanding of the role plaid by the electronic correlations in the unusual properties of the hole-doped cuprates. 
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Appendix A: Full information about the quantum problem when defined in the LWS-subspace 

In this Appendix we confirm that full information about the quantum problem can be achieved by defining it in 
the LWS-subspace spanned by the energy eigenstates that are both LWSs of the jy-spin and spin algebras. By that 
we mean that all A/'-electron operator matrix elements between energy eigenstates such that at least one of them is 
a non-LWS can be evaluated exactly in terms of a corresponding quantum problem involving another well-defined 
A/'-electron operator acting onto the LWS-subspace. Here N = 1,2, ... refers to one, two, or any other finite number 
of electrons so that the expression of the general A/'-electron operator under consideration reduces to an elementary 
creation or annihilation electronic operator (A/" = 1) or involves the product of two or more such elementary operators 
[M>2). 

In addition, we show here that use of the model global symmetry provides a simple relation between the energy of 
any non-LWS and the corresponding LWS so that all contributions to the physical quantities from non-LWSs can be 
evaluated by considering a related problem defined in the LWS-subspace. The use of such a symmetry reveals that 
for hole concentration a; > and spin density m > the ground state is always a LWS of both the 77-spin and spin 
algebras. For simplicity here we consider matrix elements between the ground state and a non-LWS. Those appear 
in Lehmann representations of zero-temperature spectral functions and correlation functions, but similar results can 
be obtained for matrix elements between any excited energy eigenstates. 

Expressions (|8))- (ITT|) for the c fermion, yy-spinon, and spinon operators refer to the LWS-subspace. A similar 
representation can be used for instance for the HWS subspace, which is spanned by all energy eigenstates with 
Sa = where a = r],s. (There are also two mixed subspaces such that S,^ = ±S^ and Sg = TSg.) The HWS 
representation is suitable for canonical ensembles referring to electronic densities larger than one and negative spin 
densities since then the ground states are HWSs of both the ?7-spin and spin algebras. The studies of this paper refer 
to the LWS-subspace. In this Appendix we provide the expressions for the c fermion, yy-spinon, and spinon operators 
that are suitable for the HWS subspace. 

1. How to map the problem onto the LWS-subspace 

We start by confirming that there is a well-defined A/'-electron operator G)jv such that any matrix element 
(/I OatIV'gs) of a A/'-electron operator Oj\f where |/) is a non-LWS of the ?7-spin algebra and/or spin algebra and \ipGs) 
the ground state of the Hubbard model on the ID or square lattice can be written as (/| Oj\f\ipGs) = {f-LHS\Q^f\iJjGs)- 
Here \ f.LWS) is the LWS that corresponds to the state |/). Within the two SU{2) algebras the latter state can be 
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expressed as, 



n \f-LWH}, (Al) 



a— 7]. s 



where 

Ca =^L„,.i/„0+ n i[La + l-l], (A2) 



1/2 



is a normalization constant, La _i/2 < = 25*0, and the 77-spin flip (a = 77) and spin flip {a — s) operators S*^ 
are the ofl-diagonal generators of the corresponding SU{2) algebras provided in Eq. These operators remain 

invariant under the electron - rotated-electron unitary transformation and thus as given in that equation have the 
same expression in terms of electron and rotated-electron creation and annihilation operators. 

As confirmed in Section IV, for a hole concentration x > and spin density to > the ground state is a LWS of 
both the ry-spin and spin SU{2) algebras and thus has the following property, 



Sa IV^Gs) = ; a^r), s. (A3) 



The operator Q^/ is then such that. 



{f\0^r\^PGs) = {f-LWS\ J] -1/2 o^l^,.^) 

a— 77, s ^ 

- {f.LHS\e^f\i;Gs)■ (A4) 

By the suitable use of Eq. (jA3p . it is straightforward to show that for iv^^^_i/2 > and/or Ls _i/2 > the operator 
Sj^f is given by the fohowing commutator, 



a— ?7, s * CK— s 



and by, 



= Oa^ , (A6) 



for i^._i/2 ~ Lg _if2 — 0. If the commutator on the right-hand side of Eq. (jA5l) vanishes then the matrix element 
(/I OJ^^\^|}Gs) under consideration also vanishes. 

Let Ef denote the energy eigenvalue of the non-LWS |/) and Ef,LWS that of the corresponding LWS \f.LWS). 
To find the relation between Ef and Ef,Lws one adds chemical-potential and magnetic-field operator terms to the 
Hamiltonian ([Ij , what lowers its symmetry. Since our description refers to the LWS-subspace, we use units such that 
the sign of the chemical potential n > is that of the hole concentration x — [N^ — N]/N^ > and the sign of 
the magnetic field H > that of the spin density m = [A^-^ — Ni]/N^ > 0. Since such operator terms commute 
with the Hamiltonian ([T]) , the rotated-electron occupancy configurations of all energy eigenstates correspond to state 
representations of its global symmetry for all densities. Moreover, the use of such commutation relations reveals that 
the energy eigenvalues Ef and Ef,LWS are related as, 

Ef—Ef,Lws+ ^ A*Q-^Q,-i/2, (A7) 

ct—rj,s 

where /x^ = 2/i is twice the chemical potential whose range is 2/i" < 2/i < [/ + ADt for D = 1, 2, 2/1*^ is the half-filling 
Mott- Hubbard gap considered in Subsection IV-C, (Xg = 2/i_B H > 0, and /is is the Bohr magneton. 

It follows from Eq. (|A7I) that Ef > Ef,Lws- For hole concentration x > and spin density to > such an inequality 
can be shown to be consistent with the ground state being always a LWS of both the 7y-spin and spin algebras. On 
the other hand, the model global symmetry requires that Ef = Ef,Lws for the half-filing and zero-magnetization 
absolute ground state. Such a requirement is fulfilled: for half filling one has that /i e (— /x", /i"). The value /i = 
corresponds to the middle of the Mott-Hubbard gap, whereas the magnetic field H vanishes for zero magnetization. 
The absolute ground state corresponds to 5^ = Ss = and Sc — /2 and hence is both a LWS and a HWS of the 
ry-spin and spin algebras. 
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Within a Lehmann representation the A/'-electron spectral functions are expressed as a sum of terms, one for each 
excited energy eigenstate. In spectral-function terms associated with non-LWSs one can then replace the matrix 
element (/| OaaIV'gs) by {f.LHS\&_^f\tpcs) provided that the excited-state energy Ef is expressed as in Eq. (|A7|) . 
Then the original A/'-electron spectral function can be written as a sum of spectral functions of well-defined A/'-electron 
operators of the form given in Eqs. (IA5[) or (|A6[) . all acting onto the LWS-subspace. 

These results confirm that full information about the present quantum problem can be obtained by defining it in 
the LWS-subspace. 



2. Operators of the three elementary objects within the HWS representation 

Within the HWS representation the c fermion operators read, 

4,c - 4„i (1 - "^,t) + ~c,,,in,^A ■ (A8) 

The ?7-spinon and spinon operators are given by Eq. (jH]) where now the operator n^.^c = fp. ^ frj,c of Eq. (jlOp involves 
the c fermion operator (|A8[) and the rotated quasi-spin operators qp, — q^, ±i qf and read, 

Qf, - (4)^; 1p,=f^r,.i-l^ (A9) 

respectively. After inverting such relations one finds that within the HWS representation the local rotated-electron 
operators can be expressed in terms of the local c fermion and rotated quasi-spin operators as follows, 



4„t = 4,iflc-^''-''fr..c). (AlO) 

The commutation and anti-commutation relations (|12l) - (fT6)) are valid for all representations and hence are the same 
for both the LWS and HWS representations. 



Appendix B: Equivalence of the quantum numbers of the present description and those of the exact solution 

for ID 

1. Relation to the quantum numbers of the exact solution for Na ^ 1 

The studies of this paper establish that the spinous (and 77-spinons) that are not invariant under the electron 
rotated-electron transformation are part of spin (and ry-spin) neutral 2z^-spinon (and 2i/-r7-spinon) composite sv (and 
7]!/) bond particles, which have a bounding (and an anti-bounding) character. Moreover, for the ID model some neutral 
2^-spinon (and 2i/-77-spinon) configurations exist that define suitable si^ (and rjiy) bond particles whose occupancies 
generate the energy eigenstates for f//4t > 0. The detailed structure of such configurations remains an open problem. 
It is studied in Ref both for the ID and square lattices for the case of the one- and two-electron subspace as defined 
in this paper, where only the c fermions and si bond particles have an active role. 

By generalizing the procedures used in Ref.^^* for the si bond particles to the remaining av bond-particle branches 
introduced in this paper, one finds that for the model on the ID lattice such objects can be associated with av 
bond-particle operators gl^ of the general form, 

g 

Here ax ^av,g is a superposition of operators belonging to a number of families, which increases for increasing values 
of the number v = 1,2, ... of spinon (a = s) or 7y-spinon {a = rj) pairs, the spatial-coordinate Xj refers to the av 
effective lattice, the index j — 1, Na^^ to the sites of that lattice, and the index g to different types of 2i/-site links 
belonging to the same family as defined in that paper for the si bond particles. The c effective lattice equals the 
original lattice and thus has length L — Na a and lattice constant a so that xj — j a for Na ^ 1- The av effective 



(Bl) 
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lattices have also length L yet the lattice constants read aau = [Na/Na^^] a> a and are such that L = Na^^ Uai, and 
Xj = j a^u where j = 1, Na^^. Furthermore, according to the same results the number of sites of the av effective 
lattice Na^^ is for ID given by expression (|42|) with I? = 1 so that Na^„ = N^i, + N^^ where Na^ = N^^ and Nau is 
the number of composite av bond particles and hence of occupied sites of the av effective lattice. The corresponding 
number of unoccupied sites iV^j^ of that lattice is provided in Eq. (|43p with N^^ = Na^^ ■ Here Ca is the number of 
"occupied sites" of the 77-spin (a = ?/) and spin (a = s) effective lattice, which refers to the subspaces with constant 
values of 5c, Srj, and Ss and is given in Eq. ([31]) where iV^ = Na^^ for ID is the total number of sites of these 
effective lattices. Each subspace with constant values of Sc and hence also of Na = [Na — SS'c] and = "iSc can 
be divided into smaller subspaces with constant values of Sc, S",,, and Ss and hence also of C,, and Cs- Furthermore, 
the latter subspaces can be further divided into even smaller subspaces with constant values for the set of numbers 
{Nrjiy} and {A^s^}, which must obey the sum-rule of Eq. (PT|) . 

Concerning the operators ax^.ai^.g of the general expressions given in Eq. (|B1[) . the simplest case refers to the si 
bond particles whose number of families is two and for instance for the iVg\ = configuration state the operator 
axj,si,g is given byi^. 



E 



6^ 





'1 






'1 




s ) 
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2 " 




^X+Xi 
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^x+xi 


"x-xi j 



(B2) 



and bx,si,i.g = [bl, gi i gj ■ Here the spinon operators are given in Eq. ([9]) and for simplicity we have omitted the 

family index rf, which reads d = I for ID so that xi corresponds to the general two-site link vector fd ; introduced in 
Rcf.18. 

The studies of Ref^iS, lead to si bond-particle operators that are shown in that paper to obey a hard-core algebra. 
For ID it follows straightforwardly that the two-77-spinons 77I bond-particle operators also obey such an algebra. 
One then considers the four-spinon s2 bond-particle and four-77-spinon 772 bond-particle operator problems and finds 
that operators of the general form (jBip can be constructed by procedures similar to those used for the si and rjl 
bond-particle operators so that the hard-core algebra is obeyed. 

We skip here the details of such derivations since the basic procedures are the same as those already used for the 
two-spinon si bond particles in Refjii^ except that the algebraic manipulations are for 1/ > 1 much more cumbersome. 
Moreover, it turns out that for interesting applications within the dynamical theory of Ref^SS, referring to the model on 
the ID lattice, such as the evaluation for finite excitation energy lo of one- and two-electron matrix elements between 
the ground state and excited states and corresponding spectral-weight distributions, only excited states involving c 
and si fermions and a finite number of sv fermions with v > 1 spinon pairs lead to finite spectral weight. The latter 
sv fermions have vanishing energy and momentum and remain invariant under the electron - rotated-electron unitary 
transformation so that the only effect of their creation and annihilation is in the numbers of si fermions and si band 
discrete momentum values. Therefore, for the excited states contributing to the dynamical theory of Refi^ only the 
c fermions and si fermions play an active role. 

The only general result this is needed for the goals of this Appendix is that alike for the two-site one-bond operators, 
composite aiy bond-particle operators of general form given in Eq. (|Bip can for ID be constructed so that when acting 
onto the LWS Hilbert space they anticommute on the same av effective-lattice site. 



{9xj,ai/T 9xj,ai'} — 1; {9xj,auT 9xj,ai/} ~ {9xj,ayT9xj,ai/} — 0, 

and commute on different sites, 

[9xj,avT 9x-, ,au] — [9xj,auT9xj,,au] ~ [9xj,au,9xj,,au] —0: j ^ j ■ 



(B3) 



(B4) 



Furthermore, such operators commute with the c fermion operators and operators corresponding to different aiy 
branches also commute with each other. 

In this Appendix we confirm that the general hard-core algebra of Eqs. (|B3[) and (|B4[) combined with the universal 
number expressions given in Section IV leads to discrete momentum values for the c and fermions that coincide 
with the quantum numbers of the exact solution for the whole LWS Hilbert subspace it refers to. 

It follows from such an algebra that one can perform an extended Jordan- Wigner transformation that transforms 
the ai/ bond particles into fermions with operators /J^ ^n/- Alike in the general expressions provided in Eq. 
such operators are related to the corresponding bond-particle operators as, 



J Xj ^OLU 



9x 



fxj .au 



where 



^j,av — E/ 



fx-/,aufx,,,av 4'j 



Vj' ,j,au 



arctan 



9x 



yy - Vj 



< 



< 27r. 



(B5) 



(B6) 
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However, for ID the coordinate Zj = Xj +i yj is such that yj = and hence reduces to the real-space coordinate Xj 
of the ah' bond particle in its av effective lattice. Therefore, for ID the phase (pji j^av can for all av branches have 
the values (j)j\j,av = and <j)j\j,av — tt only. Indeed, the relative angle between two sites of the av effective lattice in 
a ID chain can only be one of the two values. Then the av phase factor of Eq. (jB6P is such that. 



where 4>au{xj) = (f>j,au- 

The c fermion operators have the anticommuting relations given in Eq. (jl2p . which for ID read, 

{fl.c , = <^„' ; {4,, , flJ = {/.„c , = . (B8) 

Moreover, the av bond-particle operators that emerge from the Jordan- Wigner transformation associated with Eqs. 
(|B5p and (|B6p have similar anticommuting relations given by, 

{fxj.av J fx-,.av} = <^j,j' ] {fxj,au i fx-,,av\ ~ {fxj,ai' , fx-,,au\ = , (B9) 

and the c fermion operators commute with the av fermion operators and av and a' v' fermion operators such that 
av ^ a'v' also commute. 

One can introduce c fermion operators given in Eq. ([8]) associated with discrete momentum values, which for ID 
read. 





/T , = ^ > ^ e+^'i'-^' fl,^, ; /,^,, = ^ > ^ e-'"^-^^' /.,„c ; j = 1, Na ; x, = ja; L = aN, . (BIO) 



Iqj.c— fly- / . Jx^,,ci JQj.c 

For the ID Hubbard model the av fermion operators f^. given in Eq. (j24p labeled by the discrete momentum 
values Qj such that j = 1, Na^^, which are the conjugate variables of the av effective lattice spatial coordinates 
Xj, can be introduced provided that the ratio Na^^/Na involving the number Na^^ of sites of the av effective lattice 
given in Eqs. and (|43p is finite for Na — >■ oo, 

fqj,av ~ / J.J ^ ^ e ''^ J fxjt,au^ fqj,av — / ^ ^ 6 '^^ ' fxj,,au', j = ■■■■> ^a^^ i —jlau, (Bll) 

where L = a^^Na^^- 

In ID the phase factor e"^^ "" does not have any effect when operating before f^. It follows that in ID the 
expression of the Hamiltonian does not involve the phase (j)j^au- Moreover, expression of the ID normal-ordered 
Hamiltonian in terms of the c and av fermion operators reveals that such objects have zero-momentum forward- 
scattering only. This is consistent with the integrability of the model in ID and the existence of an infinite number 
of conservations laws for the limit Na — ^ oo that our description refers toi^. For the ID model the occurrence of 
such conservations laws is behind the set of av fermion numbers {Nau} being good quantum numbers^. This is in 
contrast to the model on the square lattice, for which such numbers are not in general conserved, yet the quantity 

The Jordan- Wigner transformations phases 4'j,ai, have direct effects on the boundary conditions, which determine 
the discrete momentum values qj of both the c and av fermion operators of Eqs. (jBlOp and (jBlip . respectively. In ID 
the periodic boundary conditions of the original electron problem are ensured provided that one takes into account 
the effects of the Jordan- Wigner transformation on the boundary conditions of the c ferniions and av fermions upon 
moving one of such objects around the chain of length L once. 

As discussed in Subsection IV-B, for both the model on the ID and square lattices the rotated-electron occupancies 
of the sites of the original lattice separate into two degrees of freedom only. (In the remaining of this Appendix 
we limit our considerations to the ID problem.) Those of the 2Sc sites of the original lattice singly occupied by 
rotated electrons separate into (i) 2Sc sites of the c effective lattice occupied by c fermions and (ii) 2Sc sites of the 
spin effective lattice occupied by spinous. Those of the [Na — 2Sc] sites of the original lattice doubly occupied and 
unoccupied by rotated electrons separate into (i) [Na — 2Sc] sites of the c effective lattice unoccupied by c fermions 
and (ii) [Na — 2Sc] sites of the ry-spin effective lattice occupied by ry-spinons. 

For instance, the 2v sites of the spin (and ry-spin) effective lattice referring to the occupancy configuration of one 
local 2i^-spinon composite sv fermion (and 2i^-77-spinon composite rjv fermion) correspond to the spin (and 77-spin) 
degrees of freedom of 2v sites of the original lattice whose degrees of freedom associated with the global U{1) symmetry 
found in Ref3 are described by 2v sites of the c effective lattice occupied (and unoccupied) by c fermions. 



51 



An important point is that the independent spinons and independent 7/-spinons are not part of the Jordan- Wigner 
transformations that transform the av bond particles onto av fermions. It follows that when one c fermion moves 
around its effective lattice of length L it feels the effects of the Jordan- Wigner transformations through the sites of 
the spin and ry-spin lattices associated with those of the su and tjv effective lattices occupied by sv and 77^ fermions, 
respectively. Indeed, we recall that the sites of the c effective lattice on the one hand and those of the spin and ?7-spin 
effective lattices on the other hand correspond to the different degrees of freedom of rotated-electron occupancies of 
the same sites of the original lattice. 

Since the c fermions do not emerge from a Jordan- Wigner transformation and each av fermion corresponds to a 
set of 2iy sites of the original lattice different from and independent of those of any other a'v' fermion, when a av 
fermion moves around its av effective lattice of length L it only feels the Jordan- Wigner-transformation phases of its 
own lattice associated with both the av fermions and av fermion holes so that its discrete momentum values obey 
the following periodic or anti-periodioc boundary conditions, 



iqj L 



Here the phase factor reads 1 and —1 for [Na^^ — 1] even and odd, respectively. The term —1 in \Na^^ — 1] can 
be understood as referring to the site occupied by the av fermion moving around its effective lattice and must be 
excluded. For the av fermions the unoccupied sites of their av effective lattice exist in their own right. Indeed, note 
that according to Eq. (jB5[) both the creation and annihilation operators of such objects involve the Jordan- Wigner- 
transformation phase ipj^av- As a result such a phase affects both the av fermions and av fermion holes. That justifies 
why the phase factor e*'^'^"""^^! of Eq. (jB12p involves all the Na^„ — [Nai, + Nj^^] sites of the av effective lattice 
except that occupied by the moving av fermion. Hence it involves both the [Nau — 1] sites occupied by the remaining 
fermions of the same av branch and the corresponding N^^^ av fermion holes. 

In contrast, the c fermions are affected by the sites occupied by av fermions only. Indeed, only the sets of 2v sites 
of the spin (and 77-spin) effective lattice associated with each occupied site of the sv (and rjv) v = 1,2, 3, ... effective 
lattices and the sites of the spin (and 77-spin) effective lattice occupied by independent spinons (and independent 
ry-spinons) correspond to sites of the original lattice whose degrees of freedom associated with the c fermion U{1) 
symmetry are described by the occupancy configurations of the c effective lattice. However, the independent spinons 
(and independent 77-spinons) do not undergo any Jordan- Wigner transformation so that due to the Jordan- Wigner- 
transformation phase (pj' ^ai' of each of the Nav av fermions at sites j' = 1, Na^^ of their av effective lattice the c 
fermion discrete momentum values are determined by the following periodic or anti-periodioc boundary condition. 



Jl e^(0,'+i,„.-0,',„J ^ e-E^^A^o. . (B13) 

ay j' — l 

Again the phase factor on the right-hand side of Eq. (|B13[) reads 1 and —1 for even and odd, respectively. 

The above results imply that the discrete momentum values qj of both c and av fermions have the usual momentum 
spacing qj+i — qj = 2'k jL and read, 

27r 27r 
qj^—If; j = l,...,iV,_; qj = j;Ij^ j = l,...,Na. (B14) 

However, following the boundary conditions (IB12[) [and (|B13[) ] the numbers (and I^) where j = 1,2, Na^^ 
(and j = 1,2, Na) appearing in this equation are not always integers. They are integers and half-odd integers 
for [Na^^ — 1] (and X^cny ) even and odd, respectively. Furthermore, as a result of the periodic or anti-periodic 
character of such boundary conditions these numbers obey the inequality j/""! < [Na^^ — l]/2 for both [iV^^^ — 1] odd 
and even (and the inequality < [Na - l]/2 for J2au ^au even and -[Na ~ 2]/2 < < Na/2 for Y.au '^^^ odd). 
For the one- and two-electron subspace one can separate the numbers /j*^ and of Eq. (jB14p in two terms 

corresponding to an integer number and a small deviation as /J^ = [Mj^ + ^TITtI ^"^^ ^ [-^/^ Irl^]' respectively. 
The corresponding c and si fermion discrete momentum values then read. 



qj = ^Af!; + qVNc\ A/7=j-^=0,±l,±2,...; 



= ^A/^i+g^i/A^.i; ^ 3 ~ ^ ^ Q,±l,±2, j = 1, iV,^, , (B15) 

Here (and (7°^) is given either by qJ? = or gj? = ■n[Nc/L\ (and = or (7°^ = 7r[iVsi/L]) for all j — 1, Na (and 
j = 1, Na^-^) discrete momentum values of the c (and si) band whose momentum occupancy describes a given state 
and the si effective lattice length L — Na^^ asi where a^i = L/Na^^ — [Na/Na^-^] a is the si effective lattice constant. 
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Importantly, the c fermion and ap femiion discrete momentum values obtained from our Na ^ 1 operational 
description of the quantum problem correspond to the Bethe-ansatz quantum numbers of the exact solution. Indeed, 
the discrete momentum values of Eq. (|B14p can be expressed as, 

Qj^j^Ij-, j = l,:;Na; = — J/; j = 1, 7V,^„ ; ^ — J- ■ j = 1, 7V^,,^ , (B16) 

where Ij = /|, J^" = fj'^ , and = /j"^ are the exact-solution integers or half integers quantum numbers involved 
in Eqs. (2.12a)-(2.12c) of Ref.— and defined in the unnumbered equations provided below these equations (in the 
notation of that reference v — n and j = a in J^'^ and JJ). Moreover, the numbers on the right-hand side of the 
two inequalities given just above Eq. (2.13a) of that reference correspond to Na^^j2 and respectively, so 

that these inequalities read \Jj^\ < and jJj'l < Na^^/2. That is fully consistent with the above inequality 

< [-^qq^ ^ 1]/2 where a = ri,s and Na^,, is given in Eq. (H^ . A careful comparison of the notations and 
definitions used in Ref and here confirms that there is also full consistency between the even or odd character of 
the integer numbers [A^a,„ — 1], [-^a„ — 1], and J2au -^av considered here and those that determine the integer of 
half- integer character of the quantum numbers Jj"^ = Ij'^ , Jj = /j"^, and Ij = Ij, respectively, in that reference. 

We emphasize that for the aiy fermions the discrete momentum values qj of Eq. (jB16|) are the eigenvalues of the 
translation generator in the presence of the fictitious magnetic field of Eq. (P5|) . which for ID reads Ba,y{xj) = 
'^j'=ij ^x-,,av ^i^j' ~ Xj) 6x3. Hence the corresponding exact-solution quantum numbers are the eigenvalues of such a 
translator operator in units of 2t:/ L. 

Since for ID the numbers {Na^} of ai> fermions are conserved, that the discrete momentum values qj of the 
c and av fermions are good quantum numbers is consistent with the momentum operator commuting with the 
unitary operator V'' as defined in this paper. That operator generates exact U /At > energy and momentum 

eigenstates lws^u /At) = V"^ |*lws;oo) of general form given in Eq. ^ where in ID lws^u /At) = LWSvU/it) 
with ?7/4t-dependent energy eigenvalues but [//4t-independent momentum eigenvalues from the U/4i — 00 energy 
and momentum eigenstates \^ lwS\oo) of Eq. (|58p . Hence the momentum eigenvalues are fully determined by the 
U/At > 1 physics. 

The use of the exact solution of the ID problem confirms that the momentum eigenvalues have the general form 
given in Eq. (I63p . For ID they may be written as, 

^ = + E <if NsAqf) + J2 Y.^"" - ir] N^iir) + ^ M^.-i/2 ■ (B17) 

j=l y j' = \ T)v j' = \ 

Here M^ _xii is the total number of 7]-spin-projection —1/2 ?7-spinons and the distributions Nc{qj) and Nau{qj') are 
the eigenvalues of the operators Nc{<lj) — fqj.c fqj,c and Nau{<lj) — flj^au fqj,au, respectively, which have values 1 and 
for occupied and unoccupied momentum values, respectively. One may obtain expression (|B17p from analysis of the 
ID problem for U/At ^ 1 without the use of Bethe ansatz, whose starting point is that for U/t 00 the electrons 
that singly occupy sites do not feel the on-site repulsion. Consistently, expression (jB17p is that also provided by the 
exact solution after the rapidities are replaced by the quantum numbers Ij = Ij, Jj" = fj^ ^ and = If related to 
the discrete momentum values qj by Eq. (IB16P . 

That the physical momentum (jB17P is for U/At > additive in the c and av fermion discrete momentum values 
follows from the latter being good quantum numbers whose occupancy configurations generate the energy eigenstates. 
It follows from the direct relation to the thermodynamic Bethe ansatz equations of Ref.^^ that the c fermions obtained 
here from the rotated electrons through Eq. ([5]) and the aiy fermions that emerge from the Jordan- Wigner transfor- 
mations of Eq. (jB5P are for the ID lattice the c pseudoparticles and av pseudoparticles, respectively, associated in 
Ref.'^'' with the Bethe-ansatz quantum numbers. The momentum quantum numbers of Eq. (jB16p are precisely those 
given in Eq. (A.l) of Rcf.'^^ and the numbers Na^^ of Eq. (|^^ equal the numbers N*^ defined by its Eqs. (B.6), 
(B.7), and (B.ll) with the index c replaced by rj. Furthermore, the ?7-spinons and spinous considered here are for 
ID the holons and spinous of that reference, respectively. Also the independent 77-spinons and independent spinous 
are for ID the Yang holons and HL spinous, respectively, of Refj^l. (In the notation of that reference HL stands for 
Heilmann and Lieb.) 

2. Relation to the algebraic formulation of the exact solution 

We just confirmed that for ID the discrete momentum values of the c fermion operators f^. ^, which within our 
description emerge from the electron - rotated-electron unitary transformation, and those of the av fermion operators 
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fqj av^ which emerge from that transformation and an extended Jordan- Wigner transformation, equal the quantum 
numbers of the exact Bethe-ansatz solution. Such a result was obtained in the limit Na 3> 1 that the description 
considered in this paper and in Rcfs.— refers to. However, such a connection corresponds to the quantum numbers 
only and the relation of the c and au fermion operators to the exact solution of the ID model remains an open 
problem. 

The relation of the building blocks of our description to the original electrons is uniquely defined yet corresponds 
to a complex problem. Such building blocks are the c fermions, 77-spinons, and spinous. For U /At 3> 1 the rotated 
electrons become electrons and the c fermion creation operator ft ^ becomes the quasicharge annihilation operator 

Cr of Ref.— . Therefore, in that limit the c fermions are the "holes" of the quasicharge particles of that reference 
whereas the spinous and 77-spinons are associated with the local spin and pseudospin operators, respectively, of the 
same reference. The transformation considered in Ref.^^ does not introduce Hilbert-space constraints. It follows that 
suitable occupancy configurations of the objects associated with the local quasicharge, spin, and pseudospin operators 
considered in that reference exist that generate a complete set of states. However, only in the limit U/At ^ 1 suitable 
occupancy configurations of such basic objects generate exact energy eigenstates. 

The point is that rotated electrons as defined in this paper are related to electrons by a unitary transformation. 
And such a transformation is such that for U/4t > rotated-electron occupancy configurations of the same form as 
those that generate energy eigenstates for U /At 1 in terms of electron operators do generate energy eigenstates for 
finite values of U /At. The c fermion, 77-spinon, and spinon operators are related to the rotated-electron operators as 
the quasicharge, spin, and pseudospin operators of Ref.^^ are related to electron operators. 

Importantly, note that the validity of the c fermion, spinon, and 77-spinon operational description constructed in 
this paper and in Refs.— for the Hubbard model on a square and ID lattices is for the ID problem independent of its 
relation to the exact solution. For the LWS subspace that such a solution refers to the validity of our operational de- 
scription follows from the transformations behind it not introducing Hilbert-space constraints. Such a transformations 
correspond to explicit operator expressions in terms of rotated-electron operators: For the c fermion operators such 
an expression is given in Eq. ([S]) and for the spinon and 77-spinon operators in Eqs. And the rotated-electron 

operators are related to the original electron operators by the unitary transformation ct ^ — V^'f ct^V whose unitary 
operator is within our description uniquely defined. For U/4t > it has been constructed to inherently generating 
a complete set of energy eigenstates of the general form LWS;U/it) = \'^lwS;od) with {\'^ LWS;oa)} being a com- 
plete set of suitably chosen U/At » 1 energy eigenstates. For ID such states are such that \^ LWS;U/it) — \^LWS;U/4t} 
where \^LWS;U/it) = \^lwS;oo) are the states of general form given in Eq. ([TO)) . 

In order to clarify the relation of the c and av fermion operators to the exact solution of the ID model rather than the 
so called coordinate Bethe ansatai^iii it is convenient to consider the solution of the problem by an algebraic operator 
formulation where the HWSs or LWSs of the 7;-spin and spin algebras are built up in terms of linear combination of 
products of several types of creation fields acting onto the hole or electronic vacuum, respectivel y^^'^^ . The model 
energy eigenstates that are HWSs or LWSs of these algebras are often called Bethe states. Here we briefly discuss how 
in ID and for Na 3> 1 the c and av fermion operators emerge from the creation fields of the algebraic formulation of 
the Bethe states. 

That the general description introduced in this paper for the Hubbard model on the square and ID lattices is 
consistent with the exact solution of the ID problem at the operator level as well confirms its validity for ID and this 
is the only motivation and aim of this Appendix. However, since that refers to a side problem of the square-lattice 
quantum liquid of c and si fermions studied in this paper, in the following discussion we skip most technical details 
that are unnecessary for its general goals. Nevertheless provided that our analysis of the problem is complemented 
with the detailed information provided in RefsJ^iii the resulting message clarifies the main issues under consideration. 

The algebraic formulation of the Bethe states refers to the transfer matrix of the classical coupled spin model, 
which is the "covering" ID Hubbard model"'''*. Indeed, within the inverse scattering method*^ the central object to 
be diagonalized is the quantum transfer matrix rather than the underlying ID Hubbard model. The transfer-matrix 
eigenvalues provide the spectrum of a set of [Na — 1] conserved charges. The creation and annihilation fields are labeled 
by the Bethe-ansatz rapidities A, which may be generally complex and are not the ultimate quantum numbers of the 
model. Many quantities are functions of such rapidities. For instance, the weights a(A) and 6(A) considered in the 
derivation of RefJ^ satisfy the free-fermion condition a(A)^ + 6(A)^ = 1. (A possible and often used parametrization 
is a(A) = cos(A) and 6(A) = sin(A).) The reparametrization A = [a(A)/6(A)] e2''(^) - [6(A)/a(A)] e-^'^t^) - U/2 where 
the constraint h{X) is defined by the relation sinh[2/i(A)] — [U /2] a(A)6(A), plays an important role in the derivation in 
the context of the quantum inverse scattering method of the non-trivial Boltzmann weights of the isotropic six-vertex 
model given in Eq. (33) of Refpi^. 

The diagonalization of the charge degrees of freedom involves a transfer matrix of the form provided in Eq. (21) 
of that reference, whose off-diagonal entries are some of the above mentioned creation and annihilation fields. The 
commutation relations of such important operators play a major role in the theory and are given in Eqs. (25), (40)- 
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(42), (B.1)-(B.3), (B.7)-(B.ll), and (B.19)-(B.22) of the same reference. The solution of the spin degrees of freedom 
involves the diagonalization of the auxiliary transfer matrix associated with the monodromy matrix provided in Eq. 
(95) of Ref.— . Again, the off-diagonal entries of that matrix play the role of creation and annihilation operators, whose 
commutation relations are given in Eq. (98) of that reference. The latter commutation relations correspond to the 
usual Faddeev-Zamolodchikov algebra associated with the traditional ABCD form of the elements of the monodromy 
matrix. In turn, the above relations associated with the charge monodromy matrix refer to a different algebra. The 
corresponding form of that matrix was called ABCDF by the authors of Ref.^^. 

The main reason why the solution of the problem by the algebraic inverse scattering methodic was achieved only 
thirty years after that of the coordinate Bethe ansatai^iii is that it was expected that the charge and spin monodromy 
matrices had the same traditional ABCD formal, consistently with the occurrence of a spin SU (2) symmetry and 
a charge (and 77-spin) SU{2) symmetry known long ago^, associated with a global 5*0(4) — [SU{2) x SU{2)]/Z2 
symmetryi*^. Fortunately, the studies of Ref.^® used an appropriate representation of the charge and spin monodromy 
matrices whose structure is able to distinguish creation and annihilation fields as well as possible hidden symmetries, 
as discussed by the authors of that reference. 

A hidden symmetry beyond S'0(4) was indeed identified recently: It is the charge global U{1) symmetry found 
in Refi^. The studies of that reference reveal that for U/it > the model charge and spin degrees of freedom are 
associated with U{2) = SU{2) x U{1) and SU{2) symmetries, rather than with two SU{2) symmetries, respectively. 
The occurrence of such charge U{2) = SU{2) x U{1) symmetry and spin SU{2) symmetry is fully consistent with 
the different ABCDF and ABCD forms of the charge and spin monodromy matrices of Eqs. (21) and (95) of Ref J^, 
respectively. Indeed, the former matrix is larger than the latter and involves more fields than expected from the global 
50(4) = [SU{2) X SU{2)]/Z2 symmetry alone, consistently with the global 50(3) x 50(3) x U{1) = [50(4) x U{l)]/Z2 
symmetry of the model on the ID and any other bipartite lattice^. 

Our general description of the model on a square and ID lattices takes such an extended global symmetry into 
account. Furthermore, for ID it was also implicitly taken into account by the appropriate representation of the charge 
and spin monodromy matrices used in Refj^. Such a consistency is a necessary condition for the c and av fermion 
operators emerging from the fields of the monodromy matrices upon diagonalization of the ID problem for Na ^ 1. 

Initially the expressions obtained by the algebraic inverse scattering method for the Bethe states include both 
wanted terms and several types of unwanted terms. The latter terms are eliminated by imposing suitable restrictions 
to the rapidities. Such constrains lead to the Bethe-ansatz equations and ultimately to the real integer and half- 
integer quantum numbers, which in units of 27r/L are the discrete momentum values that label the c and av fermion 
operators. 

After solving the Bethe ansatz for both the charge and spin degrees of freedom, one reaches the charge rapidities \j 
with [Na — 25,,] values such that j = 1, [Na — 2S^\ and the spin rapidities with [Na/2 — 5,, — Ss] values such that 
j = 1, [Na/2 - S^- S,]. Here 5,, = 5^ and Ss = SI for a HWS and 5",, = -5,^ and Ss = -5| for a LWS of both 
the ry-spin and spin algebras yet the rapidities can be extended to non-Bethc tower states provided that the number 
of their values are expressed in terms of Sjj and Ss rather than of 5^ and S^ , respectively. Equations of the same form 
as those obtained by the coordinate Bethe anstaz are reached by the algebraic operator formulation provided that one 
introduces the charge momentum rapidities kj and spin rapidities \j given by^^ z_(Aj) = [a{\j)/b{Xj)] e^'*'^^^ = e^'^^'^ 
where j = 1, [Na - 25",,] and Xj = -iXj/2 where j = 1, [Na/2 - 5,, - 5^1/2. 

Before discussing the structure of the rapidities for Na ^ 1, let us confirm that the numbers [Na — 25,,] and 
[Na/2 — Srj — Ss] of discrete charge momentum and spin rapidity values, respectively, provided by the Bethe-ansatz 
solution are closely related to the occupancy configurations of the rotated electrons that are not invariant under the 
electron - rotated-electron unitary transformation. Indeed, such numbers can be rewritten as, 

[Na - 2S^] = [2Sc + 20„] ; [Na/2 - 5„ - Ss] = [C^ + C's] , (B18) 

where the numbers O,, and Cs are those of Eqs. (pij) . (|47p . and ([?T|) . 

Here 25c is the number of elementary rotated-electron charges associated with the singly occupied sites and 20,, = 
J2i, N^-qf the number of such charges associated with the rotated-electron doubly occupied sites whose occupancy 
configurations are not invariant under the above unitary transformation. Moreover, [Cri + Cs] = X^ai/ 
the number = v N^i^ of down spins (and up spins) of such rotated-electron doubly occupied sites plus the 
number Cs ^'Yliv v N sv of down spins (and up spins) of the rotated-electron singly occupied sites whose occupancy 
configurations are not invariant under that transformation. 

The structure of the charge momentum and spin rapidities simplifies in the limit Na ^ 1^^- Then the [Na — 25,,] — 
[25c + 2Cr/\ charge momentum rapidity values separate into two classes corresponding to 25c and 20,, of these values, 
respectively. The set of charge momentum rapidity values kj such that j = 1, 2S'c are real and are related to the 
integer or half- integer quantum numbers of Eq. (|B14p such that = [27r/L]/| are the 25c discrete momentum 
values occupied by the c fermions, out of a total number Na of such values. That corresponds to the 25c elementary 
charges of the rotated electrons that singly occupy sites. 
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Moreover, for Na ^ I the [Na/2 — S^ — Sg] = [C^ + Cg] spin rapidity values separate into two classes corresponding 
to Crj and Cg of these values, respectively. The point is that the down spins and C,, up spins of the rotated 
electrons that doubly occupy sites and whose occupancy configurations are not invariant under the electron - rotated- 
electron unitary transformation combine with the 2C,, elementary charges left over by the above separation of the 
[Na — 2Sri] — [2Sc + '2Cri] charge momentum rapidity values. Therefore, C,, spin rapidity values out of [C^ + Cs] combine 
with 2Crj charge momentum rapidity values out of [2Sc + 2Cri], leading to C,, new rapidity values associated with 
the ry-spin singlet configurations. Those describe the 77-spin degrees of freedom of the rotated electrons that doubly 
occupy sites and whose occupancy configurations are not invariant under the electron - rotated-electron unitary 
transformation. 

In turn, the Cg spin rapidity values left over by the separation of the [Na/2 — 5^ — Sg] — [C,, + Cg] spin rapidity 
values describe the spin-singlet configurations associated with the rotated electrons that singly occupy sites and whose 
occupancy configurations are not invariant under the electron - rotated-electron unitary transformation. The charge 
degrees of freedom of these rotated electrons are described by the above set of charge momentum rapidity values kj 
such that j = 1, 25c. 

According to the Na 3> 1 results of Ref.-!!, the Crj ?7-spin rapidity values (and Cg spin rapidity values) further 
separate into rjv rapidities of length v = 1, 2, ... with Nr^v values (and sv rapidities of length v — 1, 2, ... with N g^ 
values). Obviously, the sum-rules C,, = v N^v and Cg = v Ngy are obeyed. The results of that reference reveal 
that sv rapidities of length v = \ are real and the imaginary part of the remaining branches of (xv rapidities where 
a = ri,s has a simple form. 

Therefore, for Na ^ 1 the original Bethe-ansatz equations lead to a system of [1 -|- C,, -I- Cg] (infinite in the limit 
Na — > 00) coupled thermodynamic equations whose solution gives the charge momentum rapidity values kj such that 
j — 1, 2Sc and the values Aj^au of the real part of the av rapidities such that j = 1, N^u as a function of the 
occupancies of the real integer or half-integer quantum numbers /| and 1°"^ of Eq. (IB14[) . Alike for the c fermions, the 
quantum numbers 1°"^ such that qj — [2t:/L] Ij"^ are the Na^, discrete momentum values occupied by av fermions, out 
of a total number Na^^ of such values. (Note that here we used a notation suitable to our rotated-electron description 
of the problem and that in Ref^ii such quantum numbers are denoted as in Eq. (IB16[) with ly replaced by n and the 
values Aj^^i, and Aj^g^ of the real part of the rapidities by and A", respectively, where a plays the role of j.) 

Summarizing the above discussion, in order to reach the c and av fermion operators the algebraic operator for- 
mulation of the diagonalization of the quantum problem starts by building up the Bethe states in terms of linear 
combination of products of the above mentioned several types of creation fields acting onto a suitable vacuum. The 
diagonalization of the charge and spin degrees of freedom involves the transfer matrices given in Eqs. (21) and (95) of 
Ref.^5., respectively, whose off-diagonal entries are some of these creation and annihilation fields. The different form of 
such matrices is consistent with the model global SO{3) x 50(3) x U{1) symmetry, which is taken into account by our 
general description both for the model on the ID and square lattices used in our studies of the square-lattice quantum 
liquid of c and si fermions. The creation and annihilation fields obey the very involved commutation relations given in 
Eqs. (25), (40)-(42), (98), (B.1)-(B.3), (B.7)-(B.ll), and (B.19)-(B.22) of Refill and are labeled by rapidities, which 
may be generally complex and are not the ultimate quantum numbers. 

However, such creation and annihilation fields and their involved algebra generate expressions for the Bethe states 
that include both wanted terms and several types of unwanted and unphysical terms. Indeed, they act onto an 
extended and partially unphysical Hilbert space, larger than that of the model. The unwanted and unphysical terms 
are eliminated by imposing suitable restrictions to the rapidities that change the nature of the fields and for Na ^ 1 
replace them by the c and av fermion operators labeled by the real integer and half-integer quantum numbers of the 
diagonalized model. Hence the Bethe-ansatz equations obtained by imposing suitable restrictions to the rapidities 
describe the relation between the rapidities and the ultimate quantum numbers associated with the c and av fermion 
operators. 

In addition to emerging from the elimination of the unwanted and unphysical terms of the Bethe states generated 
by the initial creation and annihilation fields, the c and av fermion operator algebra refers to well-defined subspaces, 
which are spanned by energy eigenstates whose number of ?7-spinons, spinous, and c fermions is constant and given by 
Na^ = [Na — 25c], Na^ = 25c, and Nc = 25c, respectively. Hence the number 25c of rotated-electron singly occupied 
sites and the numbers Na,^ and Na^ of sites of the ry-spin and spin lattices, respectively, are constant. 

As discussed in this paper, for both the model on ID lattice considered in this Appendix and the model on the 
square lattice the 5c > vacuum of the LWS subspace |0^s) given in Eq. ([^ is invariant under the electron - 
rotated-electron unitary transformation. For D — 1 the Na^ independent -1-1/2 spinous of such a vacuum are the 
spins of Na^ spin-up electrons, the Na^ independent -1-1/2 ry-spinons refer to the Na^ sites unoccupied by electrons, 
and the Nc c fermions describe the charge degrees of freedom of such electrons of the fully polarized state. 

The corresponding LWSs that span the above subspace refer to rotated electrons rather than to electrons and have 
the general form given in Eq. (j60p where the set of numbers {Na^} obey the sum-rule associated with the expressions 
given in Eq. pip , the ?7-spin (a — rj) and spin (a = s) can have values 5ct = 0, Na^/2 such that 5ct = [Na^/2~Ca\, 
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and the set of numbers {-/Va„„} of discrete momentum values of each av band are weU defined and given by Eq. 
and gS]). 

As mentioned previously, each subspace with constant values of Sc and hence also with constant values of iVa,^ = 
[iVa — 25*0] and Na^ = 2Sc, which the vacuum of Eq. ([23| refers to, can be divided into smaller subspaces with constant 
values of Sc, Srj, and Sg and hence also with constant values of and Cs- Furthermore, the latter subspaces can 
be further divided into even smaller subspaces with constant values for the set of numbers {Nrf^} and {Nsu}, which 
obey the sum-rule of Eq. ([5Tt . 

Last but not least, since in contrast to the initial creation and annihilation fields the c and av fermion operators 
generate the Bethe states free from unwanted and unphysical terms, are labeled by the quantum numbers of the 
diagonalized model, and their algebra does refer to the full LWS or HWS subspace but instead to subspaces spanned 
by well-defined types of Bethe states, there is no contradiction whatsoever between the charge ABCDF algebrapi^ 
and spin ABCD traditional Faddeev-Zamolodchikov algebra^i associated with involved commutation relations of the 
initial fields and the anticommutation relations of the c and av fermion operators provided in Eqs. (jB8l) and (jB9l) . 
Indeed, the initial creation and annihilation fields act onto an extended and partially unphysical Hilbert space whereas 
the c and av fermion operators act onto well-defined subspaces of the model physical Hilbert space. 

Upon acting onto such subspaces, the operators f^. ^ of the c fermions and f^.^av of the av fermions have the 
expected simple anticommutation relations associated with those provided in Eqs. (IBSP and (|B9I) for the corresponding 
operators /J^ ^ and /J^ ^j^, respectively. Moreover, c fermion operators commute with the av fermion operators and 
av and a'v' fermion operators belonging to different av 7^ a'v' branches also commute with each other. 

Interestingly, in ID our description may be used in the study of spectral functions involving matrix elements 
between Bethe states belonging to different subspaces with integer and half-integer character, respectively, for some 
of the quantum numbers of Eq. (IB14p . In that case the anticommutation relations of the operators /J .^^ and fl -^^v 
remain as given in Eqs. (jB8|) and (jB9|) yet effective anticommutation relations, which take into account the shifts of 
the discrete momentum values upon the subspace transitions, must be used for the operators f"^. ,, and f\.^av 

Finally, for the model on the ID lattice in the one- and two-electron subspace as defined in this paper one may 
introduce a unitary transformation that under the ground-state - excited-state transitions, shifts the discrete mo- 
mentum values qj — > qj + Qc{qj)/L of the c fermions and qj — > qj + Qsiilj)/ L of the si fermions. Here Qc{qj)/'2 
and Qsi{qj)/2 are well-defined overall phase shifts generated by such transitions. That allows the evaluation for 
finite excitation energy of one- and two-electron matrix elements between the ground state and excited states and 
corresponding spectral-weight distributions for finite excitation energy to^. In the w — >■ limit that method provides 
the well known low-energy expressions predicted long ago by other techniques^ and for finite w values successfully 
describes the unusual finite-energy one-electron spectral-weight distributions observed by photo-emission experiments 
in quasi-one-dimensional organic metals^. 

Fortunately, the ID dynamical theory of such papers refers to the one- and two-electron subspace where only the 
c and si fermions have an active and explicit role. The anti-commuting algebra of the c and si fermion operators 
and the expressions of Eqs. ^ and (|24l) apply to both the model in the square and ID lattices and hence together 
with the consistency with the ID exact solution found in this Appendix confirms the validity of the ID dynamical 
theory. The general c and si fermions considered here for the model on the square and ID lattices correspond to 
the c and si pseudoparticles of the ID problem of Refi^, respectively. In turn, the above shifted momentum values 
qj = qj + Qc{qj)/L and qj = qj + Qauiqj)/ L refer to the canonical momentum values of that reference and in its 
language label the c and si pseudofermions, respectively, of the ID pseudofermion dynamical theory. 



In this Appendix we perform the subspace-dimension summation of Eq. (j52p that runs over Sc, Srf, and Ss integer 
and half-odd- integer values. Here we consider the square lattice so that D = 2 in Eq. ([5^ . yet the derivation proceeds 
in a similar way for D = 1. The subspace dimensions have the form dr ■ Yla=7i s J^{Sa,Ma) given in Eq. where dr 
and J\f{Sa, Ma) are provided in Eq. (HH]) . Recounting the terms of Eq. ([5^ . one may choose Sjj to be the independent 
summation variable what gives. 



Appendix C: Subspace-dimension summation 




S'^=0 S„=i 



^0 



[^-S,] [^~S„] 



(CI) 



E E E 



S,=0 5^=0 Sa=S, 



2 2 
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One can then rewrite the summation (j52p in the form, 

Nll2 [Aff /2-S,] 



1 I i^\2{S^+S,) 



where U{Sri,Ss) denotes the 5,, and 5^ dependent summation over Sc as follows, 



E 



" 1 + (_i)2(s=+Sc) /Ar2 



25, 



111 ^ Q ^ Q I Vili^Q-S-l 



JV2 



25e 



1 



25e 
Sc Ss 1 

1 



[Sc ' Ss)l {Sc + SsY. [Sc -Ss- 1)! (Sc + Ss + 1)! 



1 



(iV2/2 -Sc- 5,)! (iV2/2 - 5e + 5^)! - Sc - S„ - 1))! - 5, + 5, + 1)! 



(C2) 



(C3) 



In order to evaluate S{Sri, Sg) it is useful to replace the variable 5, by fc = 5c — 5s. To simplify the notation we 
then introduce. 



S = S^ + Ss^ 5(5^, Ss); 5„ - 5, = 2?(5^, Ss) . 

Due to the parity factor, in the summation over k only the terms with k integer survive so that. 



(C4) 



k=0 



^-5 



1 



1 



k\{s-v + k)i {k-i)\{s-v + k + iy. 
1 



1 



E ^al< 



(7V2/2 - 5 - fc)! (iV2/2 + V-ky. {Nl/2 - S - k - 1)1 + V - k + 1)! 

1 1 



A:=0 



kl{S-V + fc)! {N^/2 -S-k)\ (7V2/2 + 2? - fc)! 
1 1 



fc! (5 - D + fc)! (iV2/2 - 5 - fc - 1)! {N^/2 + V-k+l)l 
1 1 



(fc - 1)! (5 - 2? + fc + 1)! (7V2/2 - 5 - fc)! (A^2/2 + 2? - fc)! 
1 1 



(fc - 1)! (5 - 2? + fc + 1)! (7V2/2 - 5 - fc - 1)! (Ar2/2 + 2? - fc + 1)! 

where now the fc summation runs over integers only. 

In order to perform the summation (|C5[) we rearrange the terms as follows. 
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1 
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1 
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Next, by using the identity, 
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N 

E 

fc=0 
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k \N -k 
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we carry out separately the summations in expression (jC7p . what gives, 

k }\Nll2-S-k) \NII2~S, 



E 
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Nll2^V 
k-\ 



NII2-V 
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E 
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Introducing these results in expression (jC7[) for S leads to, 

Nl \ W Nl \ f Nl 
,Nl/2 v) [ [nI/2 - 5 j + [nI/2 -S-2 

Nl \ r/ N^ 

,Nl/2-S-l 
S(5,I?). 
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Expression (jC2l) for Ntot can now be rewritten as, 

Mot 



E E -'^M (25, + l)(25. + l)x 



,Ni/2-{S^-Ss) 
Nl 

- Ss) - 1 



Nl/2 - (5, + Ss)) + \NII2 - (S, + S,)-2 

Nl \ ( Nl 

NII2 - (S, - Ss) + UaV2 - (S, -Ss)-l 



(C14) 



where the summations run over both integers and half-odd integers. The use of the notation (IC4[) then allows rewriting 
(IC14I) in compact form, 



A/---E E 



25 



-{S + V +1){S -V +1) xS(5,X') 



(C15) 



S„=0 Ss=0 



where the summations run again over both integers and half-odd integers. 

We can perform the summations of Eq. (IC15I) in the integers S and V instead of in 5,, and Ss- Indeed, the first 
factor cancels all the terms with S and V non-integer so that. 



lN^/2-Sr,] 



• • -(5^ and Ss both cither integers or half odd integers) 



S„=0 Ss=0 



Thus we find, 



+S 

E] ■ ■ ■ {S and V integers) . 
5=0 v=-s 



■^tot=J2 E ((5 + 1)^-2?') xS(5,I?). 



5=0 X'=-5 



The use of the result (jC13|) then leads to, 

5 

Mo*= E E ((5+1)^-2?^ 

5=0 15=^5 



N^/2~S-1 



Nl \ ( Nl 
Nll2-S)^\Nll2-S-2 

Nl \ ( Nl 
N^/2-V + l) + [n^/2~V-1J\ 



(C16) 



Replacing the variable 5 by 5' = 5 + 1 we reach a more tractable expression for Aftc 

s'-l 

^tot^ E E ns',D), 

5' = 1 I5=-5' + l 



(C17) 



where 



T{S',V) = {{S')^ -V^) x-S{S'-l,V) 



= is" - v 



NII2 - V 



iV2 



NII2-S'^\) \NII2-S' -\ 



iV2 



(C18) 



Nf. 



Nl/2-S'JWNl/2-V + lJ \Nl/2-V-l 
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is completely symmetric in the summation variables. 

Since T{S' , T) = ±S') ~ 0, we can extend the summation over T) of Eq. (jC17|) to 2? = ±5'. We then formally extend 
the summation over S' to 5' = because the corresponding term vanishes: T{S' = 0,1? = 0) =0. Futhcrmore, 
T{±S' ,V) = T{S', ±1?) = T{S' .V), and due to the symmetry S' ^ V we can write, 



s'-l 



NI/ 2+1 



S' = 1V=-S' + 1 S' ,V=-{Nl/2+l) 

Let us introduce the numbers p and q such that, 

S' + Nl/2 + l=p ^ S' =p- {Nl/2 + 1) 

V + Nll2 + l = q ^ V = q- {Nll2 + 1) . 

The use of l|C19p then allows rewriting (IC17|) as, 

Nl+2 



\ E HNl + 2-q)-p{Nl + 2-p)] 



p,q=0 



9-1 



p 



,P-2. 



p-1 



^ * a 

q 



q-2 



This expression can be simplified noticing that. 



N 
x-2 



Replacing in Eq. (jC20|) one then finds, 

N^+2 



Mtot = \Y. + 2-q)-p{Nl + 2-p) 

P:l] = . 

Nl+2 , 
p,q=0 ^ 

=¥{ e^(^.^+2-,)(Me 

L g=0 ^ p=0 

Finally, the use of the identities. 



= -2 



iV2 
q-l 



N 
x-l 



N + 2 

X 



9-1 



N^ + 2 
P 

nI+2 



p-i 



N 



N+2 

E 

fe=0 

N 



N 
k - 1 



N 

E 

fe=0 

N+l 

E 

fc=i 



iV2 + 2 

nI+2 



Ar2 

p-i 



+ 2 



(q^p) 



q=0 



N„+2 



p=0 



7V2 
* a 

p-i 



= 2 



*:=0 



E 

k=l 



N 
k-1 

m 



E ( y 



fc'=0 



N-2 



{k-l)\{N -k - 1)! 



7V(7V-1) E 



iV-2 
fc- 1 



= 7V(iV- 1)2 



N-2 



and 



/c — — 1 k ' 1 — 



N 



Y,{k' + l){N-k' + l)[l] = E [fc'(A^-fc') + (A^ + l)](r, 



TV 



fe'=0 ^ ' k'=0 

-,N-2 , fAT , T\oN_nN-2\ 



N 



= N{N - 1)2'^-^ + {N + 1)2^ = 2^-^ [N{N - 1) + A{N + 1)] 



(C19) 



(C20) 



(C21) 
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leads to, 

AA,„, ^ i I [N^ + 3iVf + 4] 2^a-2 X 2^'+2 - (iV^ + 2){NI + 1)2^' x 2^' } = 122^- x 2 = 4^' , (C22) 
which is the desired result. 
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